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Abstract 

Within  a  defined  domain,  geophysical  phenomena  often  are  characterized  by  smooth  continuous 
power  spectral  densities  having  a  negative  power  law  slope  dependence.  Frequently,  Fourier 
transform  analysis  has  been  employed  to  generate  synthetic  scenes  from  pseudorandom  arrays 
by  passing  the  stochastic  data  through  a  Fourier  filter  having  a  desired  correlation  structure  and 
power  spectral  dependency.  This  report  examines  the  possibility  of  producing  synthetic 
structure  by  invoking  autoregression  analysis  in  conjunction  with  the  Fourier  method.  Since 
computations  that  apply  multidimensional  fast  Fourier  transforms  to  large  data  arrays  consume 
enormous  resources  and  time,  the  goal  of  this  study  is  to  seek  an  alternative  method  to  reduce 
the  computational  burden.  Future  editions  of  the  Phillips  Laboratory  Atmospheric  Radiance 
Codes  will  feature  an  ability  to  calculate  structured  radiance.  The  methods  explored  herein 
provide  a  process  that  can  complement  or  in  some  cases  supplement  methods  presently  being 
used.  Line  of  sight  computations  using  the  three-dimensional  database  will  be  presented. 
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Figures 


1 .  Title  Page 

ji 

2.  Stepping  Slones  -  References  from  prfvious  work 


3.  Linear  Difference  Equation  -  An  autoregressive  moving  average  (ARMA)  model  for  a 
discrete  spatial  series.  v(n) ,  that  approximates  deterministic  and  stochastic  processes  can  be 
represented  by  the  filter  linear  difference  equation: 


4.  Power  Spectra  Model  -  A  z  transform  analysis  of  the  difference  equation  shows  that  the 
AR\IA  power  spectral  density  is  contained  in  the  equations  in  the  slide. 


5.  Yule-Walker  Equations  -  The  relationship  between  the  autocorrelation  sequence  and  the 
pure  autoregressive  model  is  contained  in  the  equations  in  the  slide. 


6.  Levinson  Recursion  Algorithm  -  A  standard  "'Levinson”  algorithm  that  takes  advantage  of 
the  Hermitian-Toeplitz  matrix  equation  was  employed  to  solve  for  the  AR  parameters.  These 
are  presented  in  the  equations  in  the  slide. 


7.  PSD  of  Simulated  Series  -  the  simulated  PSD  is  obtained  from  the  formula  given  in  the 
slide.  For  a  given  altitude,  the  forward  and  backward  method  solves  the  squares  problem 
given  in  the  slide. 


8.  PSD  and  ACF  Models  -  In  practice,  a  given  one-dimensional  single-sided  PSD  was  con¬ 
structed  from  the  PSD  and  ACF  power  law  models  as  shown  in  the  slide. 


9.  Coherence  Length  Parameter  -  the  equivalent  width  or  correlation  length  as  a  function  of  the 
PSD  may  be  written  as  presented  in  the  slide. 
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10.  PSD  -  The  model  autocorrelation  function  used  in  this  report  leads  to  an  isotropic  two- 
dimensional  power  spectral  density  function  given  by  the  equation  in  the  slide.  The 
simulated  data.  S(t,  ?),  for  a  layer  is  described  by  the  second  equation  in  the  slide. 


11.  Transverse  Correction  -  In  pursuit  of  avoiding  the  time-consuming  convolution  integral,  an 
approximation  was  developed  that  modified  the  PSOik^.k/,)  at  =0.  This  is  given  by  the 
equadon  in  the  slide.  The  horizontal  two-dimensional  discrete  Fourier  Transform  used  in  the 
syntheses  of  the  database  results  in  a  well-known  rectangular  window  sampling  effect. 
Assuming  negligible  truncation  of  the  autocorrelation  funcdon  in  the  line-of-sight  direction, 
truncation  in  the  transverse  direction  leads  to  side-lobe  leakage  of  the  frequency  response 
funcdon  of  the  digitized  rectangular  window.  In  the  transverse  direction,  the  effects  of 
leakage  can  be  computed  as  the  convolution  of  the  “theoretical”  PSD[k(..k()  with  the 
square  of  the  magnitude  of  the  frequency  response  of  a  digital  sine  function. 


12.  Line-of  Sight  Correction  -  Since  the  constmeted  database  has  a  finite  and  relatively  large 
line-of-sight  data  spacing,  a  natural  truncadon  of  the  line-of-sight  power  spectral  density 
funcdon  must  occur  that  would  cause  significant  aliasing  of  the  line-of-sight  PSD's  for  any 
particular  realization.  Computational  “aliasing”  of  the  line-of-sight  PSD's  was  added  to  the 
“theoretical”  PSD's  by  the  equations  in  the  slide. 


13.  Conclusions  Slide.  Self-contained 


04/18/94 


Memo 


Backup  Table 


It  is  anticipated  that  the  following  Table  will  be  presented  in  answer  to  questions. 

I .  The  chosen  geometry  is  a  segment  of  a  sphere  that  consists  of  a  volume  70  km  in  height  by 
approximately  3000  kra  in  the  sensor  LOS  direction  by  approximately  5 1 .2  km  in  the 
direction  transverse  to  the  line-of-sight.  A  spatial  resolution  element  has  an  incremental 
altitude  of  0. 1  km  by  an  incremental  transverse  distance  of  approximately  0.1  km  by  an 
incremental  LOS  distance  of  approximately  15.7  km.  The  actual  data  spacing  in  the 
horizontal  dimension  increased  slightly  with  altitude  so  that  points  in  the  vertical  dimension 
stayed  on  radial  lines.  The  altitude  range  accounts  for  dominant  structure  effects  (50-120 
km).  Also  its  vertical  and  transverse  dimensions  are  sufficient  to  fill  the  field  of  a  limb¬ 
viewing  two-dimensional  sensor  array.  The  vertical  and  transverse  spatial  resolution  is 
governed  by  the  vertical  and  horizontal  correlation  lengths  of  the  temperature  variations  but 
also  is  chosen  to  provide  high  field-of-view  spatial  resolution.  The  geometrical  LOS 
dimensions  are  chosen  to  provide  sufficient  length  to  describe  the  integrated  path  radiance. 
Table  I  summarizes  the  database  specifications. 
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Backup  Figures 


It  is  anticipated  that  the  following  Figures  will  be  presented  in  answers  to  questions. 

1.  Database  Geometry.  Figure  la  illustrates  the  volume  of  the  spherical  segment  defined  by  the 
database.  The  radial  altitude  range  starts  at  50  km  and  ends  at  120  km.  At  50  km  altitude  the 
curved  length  in  the  LOS  direction  measures  3014  km  while  the  curved  length  in  the 
transverse  direction  measures  51.2  km.  These  distances  increase  slightly  with  altitude  due  to 
the  curved  geometry.  Figure  lb  illustrates  the  resolution  volume  at  50  km.  The  radial 
altitude  resolution  is  0. 1  km  while  the  resolution  of  the  curved  lengths  in  the  LOS  and  trans¬ 
verse  directions  is  15.7  km  and  0.1  km  respectively.  Again  the  horizontal  resolution  lengths 
increase  slightly  with  altitude. 

2.  Computational  scheme.  Correladons  along  the  vertical  lines  are  calculated  (one  layer  at  a 
time)  and  then  the  horizontal  correlation  calculation  is  applied  to  each  layer. 

3.  Horizontal  correlation  lengths  plotted  as  a  function  of  altitude.  A  logistic  dose  response 
equation  is  fitted  to  the  data.  The  outside  curves  represent  the  95%  confidence  limits. 

4.  Vertical  correlation  lengths  plotted  as  a  function  of  altitude.  A  logistic  dose  response 
equation  is  fitted  to  the  data.  The  outside  curves  represent  the  95%  confidence  limits. 

5.  Root  mean  square  values  of  the  relative  temperature  fluctuations  plotted  as  a  function  of 
altitude.  A  rational  polynomial  is  fitted  to  the  data.  The  outside  curves  represent  the  95% 
confidence  limits. 

6.  Graphs  of  Vertical,  Transverse,  and  Line-of-Sight  one-dimensional  Power  Spectral  Densities 

plotted  against  spatial  frequency.  The  left  panel  plots  the  vertical  PSD’s  for  the  theoretical 
model  (light  gray  solid  curve),  autoregressive  “generator”  model  (dashed  curve),  and  the 
autoregressive  vertical  PSD  from  the  simulated  database  (dotted  curve).  The  middle  panel 
plots  the  transverse  PSD's  for  the  theoretical  model  (light  gray  solid  curve),  autoregressive 
PSD  from  the  simulated  database  (dotted  curve),  and  a  periodogram  from  the  simulated 
database  (solid  curve).  The  right  panel  plots  the  line-of-sight  PSD's  for  the  theoretical  model 
(light  gray  solid  curve),  autoregressive  PSD  from  the  simulated  database  (dotted  curve),  and 
a  periodogram  from  the  simulated  database  (solid  curve).  The  plots  are  for  50.6  km  altitude 
and  one-dimensional  high  frequency  spectral  slope  (5)  =  -5/3,  where  o2  =  8.34xl0'4,  L^v  = 
1.72  km,  L(;h  -  32.9  km.  The  vertical  spatial  spacing  =  0.1  km,  transverse  spacing  =  0.1  km, 
and  line-of  sight  spacing  =  15.7  km.  The  curves  show  log-log  plots  of  PSD's  measured  in 
{ST  /  Temperature^'  ,  ^  , 

- — - ; - and  wavenumber  measured  in  km'L  The  vertical  autoregressive 

Wavenumber  ® 

PSD  for  the  simulation  was  computed  by  calculating  a  6  coefficient  forward  and  backward 
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19.  PSD  estimates  of  the  horizontally  scanned  simulated  atmosphere  from  Autoregressive 
analysis. 

20.  Two-Dimensional  image  temperature  structure  scene  from  3-D  simulation.  Artificial 
temperature  sensor  located  at  85  km. 
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autoregressive  model  of  the  vertically  simulated  data  from  layers  extending  from  50.0  km  to 
50.6  km  in  the  forward  direction  and  layers  extending  from  5 1.2  km  to  50.6  km  in  the  back¬ 
ward  direction.  Note  that  the  “transverse”  direction  is  defined  as  the  horizontal 
perpendicular  to  the  “line-of-sight”  direction.  The  transverse  and  line-of-sight  PSD’s 
represent  the  1-D  component  of  the  2-D  horizontal  PSD  having  a  high  frequency  horizontal 
spectral  slope  =  -8/3.  Their  autoregressive  PSD’s  result  from  computing  a  10  coefficient 
autoregressive  model  of  the  transverse  and  line-of-sight  simulated  data  in  the  layer  at  50.6 

km.  The  transverse  Nyquist  frequency  is  approximately  l/(2x0.1)  =  5  km'^  and  the  line-of- 
sight  Nyquist  frequency  is  approximately  l/(2xl5.7)  =  0.032  km’^. 

7.  One  dimensional  periodograms  of  the  first  four  samples  of  the  line-of-sight  simulated  data 
in  the  layer  at  50.6  km. 

8.  Same  as  figure  6  except  altitude  =  1 19.3  km,  =  0.0277,  =  2.42  km,  and  L^h  =  49. 1 

km.  The  transverse  spacing  is  0.101  km  and  LOS  spacing  is  15.9  km. 

9.  Histogram  representing  the  probability  distribution  function  of  the  simulated  data  at  50.6 
km.  Solid  line  is  the  theoretical  Gaussian  pdf.  The  variance  of  the  data  is  8.35xl0"4. 

10.  Cumulative  distribution  function  of  the  simulated  data  at  50.6  km  compared  to  the 
theoretical  curve. 

11.  Sample  page  from  the  database  explanation  file. 

12.  Alias  branches  of  estimated  “theoretical”  LOS  PSD. 

13.  ‘Theoretical”  transverse  PSD  versus  a  512  point  PSD  “corrected”  for  truncation  of  the 
autocorrelation  function.  The  plot  has  the  same  input  parameters  previously  given  for  50.6 
km  altitude.  The  one-dimensional  high  frequency  spectral  slope  is  (S)  =  -5/3,  <7^  =:  8.34x10"*, 
L^i,=32.9  km,  and  transverse  spacing  =  0.1  km. 

14.  Approximation  correction  for  rectangular  window  discrete  Fourier  Transform  digital  sine 
function. 

15.  Viewing  geometry  of  an  imaginary  temperature  sensor  platform  located  at  an  altitude  of  85 
km. 

16.  Integrated  temperature  structure  along  selected  pixel  rows  comprising  zenith  look  angles 
from  the  minimum  of  94.60  to  the  maximum  of  960 

17.  Integrated  temperature  structure  along  selected  pixel  columns  comprising  horizontal  look 
angles  firom  -.01870  to  +.01870. 

18.  Fourier  periodograms  of  the  horizontally  scanned  simulated  atmosphere  for  the  same  look 
angles  of  Figure  16. 
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Applicable  to  Radiance  Structure  Models 
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Model  for  Smooth  Isotropic  Two-Dimensional 
Power  Spectra,  PL-TR-93-2224 


L-inG9r  PifforGncG  Ec|U9tion 
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Power  Spectra  Model 
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Yule  Walker  Equations 
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Levinson  Recursion  Algorithm 
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PSD  of  Simulated  Series 
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PSD  and  ACF  Models 
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ACF(x) 


Coherence  Length  Parameter 
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2-D  P! 
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PSD  is  adjusted  to  account  for  rectangular  windowing  of  decretized 
samples.  This  is  the  "Sine"  function  effect. 

For  low  Nyquist  frequency  along  the  LOS  direction,  the  PSD  is  adjusted  to 
account  for  an  aliasing  effect. 
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LOS  correction  for  small  Nyquist 

frequency 


Autoregressive  synthesis  (6  coef.)  accurately 
models  1-D  vertical  geophysical  PSD's 
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Generation  of  entire  3-D  database  in  2.5  Hrs. 


Database  Specifications 
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LOS  192  15.7  3014.  .00033  .0318  (LOW) 
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Figure  2.  Calculational  scheme.  First,  correlations  along  the 
vertical  lines  are  calculated  (one  layer  at  a  time)  and  then  the 
horizontal  correlation  calculation  is  applied  to  each  layer. 
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Horizontal  Correlation  Length  (km) 

Eqn  8013  y=a+b/(1-!-(x/c)*^)  [LogisticDoseRspj 
r2=0.g9 1459845  DF  Adj  r2=0.989182471  FitStdErr=2.07219665  Fstat=31 9. 167443 

a=32.764354  b=52.968001 
0=135.07913  d=-6.4666086 


50  100  .  150  200  250  300 

Altitude  (km) 


Rank  18  Eqn  8013  y=a+b/(1+(x/c)‘^)  [LogisticDoseRsp] 


Coef  Det 
0.9914598452 


DF  Adj 

0.9891824705 


Fit  Std  Err 
2.0721966466 


F-value 

619.16744310 


Farm  Value 
a  32.76435426 

b  52.96800100 

c  135.0791315 

d  -6.46660858 


Std  Error 
1.295628218 
2.144297606 
2.234742679 
0.680074818 


t-value 

25.28839199 

24.70179552 

60.44504933 

-9.50867230 


95%  Confidence  Limits 
30.01210311  35.51660541 

48.41295545  57.52304655 
130.3319571  139.8263059 

-7.91126428  -5.02195287 
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v 


Lev  (km) 


Vertical  Correlation  Length  (km) 

Rank  4  Eqn3013  y=a+b/(1+(x/c)‘^)  [LogisticDoseRsp] 
r2=0.977486661  OF  Adj  r2=0.9721 89405  FitStdErr=0.546523677  Fstat=260.508669 

3=1.7292654  b=3.3478628 
c=  161. 52008  d=-7.9024983 


Rank  4  Eqn8013  y=a+b/(1+(x/c)'^)  [LogisticDoseRsp] 


Coef  Det 
0.9774866610 


DFAdj  r2 
0.9721894048 


Fit  Std  Err 
0.5465236769 


F-value 

260.50866885 


Farm  Value  Std  Error  t-value  95%  Confidence  Limits 

a  1.729265412  0.210585414  8.211705540  1.286422835  2.172107988 

b  8.347862811  0.536310568  15.56535207  7.220048927  9.475676694 

c  161.5200848  3.853308016  41.91725244  153.4169175  169.6232522 

d  -7.90249825  1.339177504  -5.90100881  -10.7186706  -5.08632592 


fl(r  V 
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/3  P 


'^6T/T 


(5T/T)Rf^s  Values 
[LgstcDoseRsp]  y=a+b/(1+(x/c)'') 

r2=0.99459057  DF  Adj  r2=0.99026303  FitStdErr=0.00972201 58  Fstat=367.72483 
3=0.0087695291  b=0.31 03953 
c=11 8.04806  d=-3. 1509968 


Rank  931  Eqn8013  [LgstcDoseRsp]  y=a+b/(1+(x/c)‘^) 


Coef  Det 
0.9945905715 


DF  Adj  r2 
0.9902630287 


Fit  Std  Err 
0.0097220158 


F-value 

367.72482563 


Farm  Value 
a  0.008769529 

b  0.310395297 

c  118.0480550 

d  -3.15099683 


Std  Error 
0.009675816 
0.024196036 
5.745293881 
0.473460616 


t-value 

0.906334825 

12.82835344 

20.54691326 

-6.65524591 


95%  Confidence  Limits 


-0.01500493 

0.250943187 

103.9312859 

-4.31433753 


0.032543990 

0.369847408 

132.1648241 

-1.98765613 


Date  Time  File  Source 

Jun  1,1994  10:59:59  AM  d:\tcwin2\rms1.pm 
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Vertical  Transverse  LOS 


{^.uq/_(x/jL5)}  aSd 


{^.uq/_(jyX5)}  QSd 
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OAc  tCvP  Wavenumber  (krn’ ' )  Wavenumber  (km'  ^ )  Wavenumber  (km 

IP  /  ^  Altitude  =  120  km,  1-D  Slope  =  -5/3 
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Horizontal  Transverse  PSD 


Transverse  PSD 
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Plane  Wave  Scintillation  in  an  Onion  Skin  Model 
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Definition  of  Onion  Skin  Geometry 
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Onion  Skin  Model  and  the  Log  Amplitude  Variance 

The  Rytov  solution  for  a  plane  wave  originating  at  L  and  observed  at  s=0  (point  detector) 
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0.56fe^  {sec(v|/  )}'^‘ (z)  (z  -  z,  )^rfz 


Onion  Skin  Model 
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(Re  +  ^  (Re  +  ^I  )^  Sin^(\|/  ) 


Define  the  altitude  weighting  function, 


Wp{z,\\f)  =  0.56it/Hsec(\}/)}  /^(z  -  zj) 


Plane  wave  altitude  weighting  function  for  zenith  angle  of  60 
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Z(km) 


Plot  of  plane  wave  altitude  weighting  function  for  zenith  angle  of  89 


Properties  of  W(z,\|/) 
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Finite  Aperture  Effects 
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Path  weighting  runction  for  plane  wave  scintillation.  The  solid  curve  is  the  usual  5/6 
power  weighting  function. 
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receiver  at  an  altitude  of  12.5  kni. 


Onion  Skin  Model  Predictions  for  Airborne  Measurements 
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results  than  an  actual  profile  that  had  a  strong  layer  in  the  first  kilometer 
of  the  path. 


Inner  Scale  Effects 
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The  effecl  of  inner  scale  is  shown  for  a  1  min  receiver.  The  calculation  uses  the  Clear  1 

profile  and  a  receiver  at  12.5  km  altitude. 
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Comparison  to  Data 
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Profiles  used 


530 


14031 

LAUNCH:  05-06-06  01:19:22  UT 


531 


TEMP  ("C)  LOG  (C^) 


Predictions  of  with  Rytov/onion  skin  model 
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5  cm  Aperture,  Clear  I 
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Comparison  of  Burke's  astronomical  data  (Burke, JJ.,  JOSA,60,  1262-1263,1960)  with  predictions 

using  CLEAR  1  and  an  observer  at  ground  (1.216  km  MSL) 


Summary/Conclusions 

Onion  skin  formalism  provides  a  model  for  horizontal  structure  of  turbulence 
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COMPARISON  OF  A  MODEL  DESCRIBING  PROPAGATION  THROUGH  OPTICAL 
TURBULENCE  (PROTURB)  WITH  FIELD  DATA 


Robert  W.  Smith 

U.S.  Army  Test  and  Evaluation  Command 
Ft.  Belvoir  Meteorological  Team 
Ft.  Belvoir,  Virginia  22060 

Jennifer  C.  Ricklin 
U.S.  Army  Research  Laboratory 
Battlefield  Environment  Directorate 
White  Sands  Missile  Range,  New  Mexico  88002 

Kevin  E.  Cranston  and  John  P.  Cruncleton 
Physical  Science  Laboratory 
New  Mexico  State  University 
Las  Cruces,  New  Mexico  88003-0002 


ABSTRACT 

The  PROTURB  model,  developed  by  the  Battlefield  Environ¬ 
ment  Directorate  of  Army  Research  Laboratory,  calculates  an 
estimate  of  optical  turbulence  strength  and  its  effects  on  vis¬ 
ible  and  infrared  imaging  and  laser  system  performance.  The 
PROTURB  model  was  earlier  compared  to  data  obtained  from  the 
REBAL  92  field  test  held  in  Bushland,  Texas,  and  data  supplied  by 
the  iUrmy  Test  and  Evaluation  Command  Ft.  Belvoir  Meteorological 
Team^.  Recent  PROTURB  updates  included  modifying  the  radiative 
flux  portion  of  the  energy  balance  model.  Here  we  present  a  new 
comparison  using  data  from  three  locations  that  differ  greatly  in 
meteorological  conditions,  including  hot  and  moist  conditions, 
hot  and  dry  conditions,  and  moderate  climatic  conditions.  There 
also  exists  a  considerable  difference  in  site  characteristics 
such  as  ground  conditions  and  vegetation.  This  latter  data  set 
was  collected  during  the  summer  of  1993  at  the  following  sites: 
Eglin  AFB,  Florida,  Ft.  Hunter-Liggett ,  California,  and  Ft.  A.P. 
Hill,  Virginia.  We  discuss  the  data  collection  and  reduction 
effort,  and  show  how  this  data  compares  with  PROTURB  predictions. 


2 . INTRODUCTION 

The  PROTURB  model  consists  of  five  distinct  modules;  the 
PROTURB  interface,  surface  boundary  module,  planetary  boundary 
layer  module,  imaging  statistics  module  and  the  laser  propagation 
module.  The  PROTURB  user-friendly  interface  controls  the  input 
of  data,  the  calculation  procedure  and  the  display  of  the  output 
data  in  both  numerical  and  graphical  form.  The  surface  boundary 
layer  module  is  based  on  a  radiative-energy  flux  balance  model^ 
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that  in  turn  drives  a  Monin-0 bukhov  similarity  theory  model^  for 
unstable  atmospheric  conditions.  This  module  is  responsible  for 
the  optical  turbulence  profile  that  extends  from  10  times  the 
roughness  length  up  to  the  Obukhov  height.  The  surface  boundary 
layer  module  takes  into  consideration  the  effects  of  moisture  and 
gives  approximations  during  stable  atmospheric  conditions  and  in 
and  around  the  neutral  events.  The  planetary  boundary  layer 
module  is  also  based  on  Monin-Obukhov  similarity  theory  and  the 
Ekman  spiral'^.  This  module  is  responsible  for  the  continuation 
of  the  optical  turbulence  profile  from  the  Obukhov  height  to  the 
top  of  the  planetary  boundary  layer.  The  imaging  statistics 
module  is  primarily  based  upon  the  work  of  D.  L.  Fried^.  This 
module  determines  the  image  statistics  that  correspond  to  an 
image  that  has  propagated  over  the  designated  path.  The  current 
laser  propagation  module  is  based  on  Rytov  theory  developed  at 
ARL  West^.  This  module  is  responsible  for  determining  laser 
propagation  statistics  such  as  irradiance  variance  and  beam  size. 

3^  FIELD  DATA  COLLECTION 

All  field  metorological  data  were  collected  by  the  TECOM  Ft 
Bel voir  Meteorological  Team.  The  same  type  of  instruments  were 
used  at  each  location  and  proper  calibration  procedures  were 
followed.  Figure  1  locates  the  three  locations.  The  optical 


FT.  A.P.  HILL,VA 


EGLIN  AFB,FL 


FIGXniE  1  LOCATION  OF  DATA  SITES 

scintillometers  used  were  the  Lockheed  model  CN2.  All  data  col¬ 
lection  was  automated  and  included  temperature,  relative  humidi¬ 
ty,  pressure,  rain  rate  and  amount,  visibility , soil  temperature, 
long  and  short  wave  radiation,  cloud  amount,  cloud  height, 
and  spectral  atmospheric  transmission. All  data  except  clouds  and 
transmission  were  collected  in  a  single  data  set.  An  example  is 
shown  in  table  1. Clouds  and  transmission  were  reported  in  seper- 
ate  files  which  are  not  shown  here. It  must  be  pointed  out  that  in 
each  case,  the  collection  system  was  designed  for  the  test  sup¬ 
port  and  not  this  comparison  effort.  The  standard  measurements 
were  taken  at  a  location  very  near  the  scintillometer  receiver 
while  the  cloud  data  was  taken  from  our  base  site  which  could  be 
several  kilometers  away  from  the  other  measurments.  Data  from  Ft. 
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Day  Tfme  Prasa  UDO  USP  PW  T««p  RM  Soil  Vsby  PSP  PSP  RMRT  SCINT  Pl«  Pl« 

hhna  *6  dg  M/S  M/S  dg/c  X  dg/c  lui  VER  (JOT  «iB/h  VOLT  VER  OUT 

196  15  1012.8  143  0.35  0.60  23.2  94  14.9  10.54  0.00  0.00  0.00  0.016  391.5  U2.4 

196  30  1012.6  20  0.04  0.40  23.1  95  14.8  9.15  0.00  0.00  0.00  0.013  390.8  U1.7 

196  45  1012.6  52  0.09  0.40  23.1  95  14.6  10.86  0.00  0.00  0.00  0.005  391.1  U1.5 

196  100  1012.4  104  0.21  0.60  23.1  95  14.6  11.51  0.00  0.00  0.00  0.017  390.2  440.3 

196  115  1012.3  60  0.31  0.90  23.0  95  14.5  9.81  0.00  0.00  0.00  0.012  390.8  U0.4 

TABIjE  1.  FORMAT  OF  FIELD  DATA 

AP  Hill  and  Eglin  AFB  were  averaged  for  one  minute  while  10 
minute  data  was  collected  at  Ft.  Hunter  Liggett.  Data  was  col¬ 
lected  over  the  entire  24  hour  period  and  examples  of  tempera¬ 
ture,  humidity,  solar  radiation  and  are  shown  in  figure  2  to 

illustrate  the  differences  between  the  locations.  Since  much 
more  data  was  collected  than  could  possibly  be  used  in  this 
effort,  a  few  representative  days  were  chosen  from  each  location. 
Complete  details  are  available  from  the  Ft.  Belvoir  team. 


TEMPERATURE  SOLAR  RADIATION 


FIGURE  2 


GRAPHICAL  EXAMPLES  OF  SELECTED  DATA  EXAMPLES 


4 . PROTURB  METHODOLOGY 

The  field  data  was  reviewed  and  10  to  15  minute  averages 
were  determined.  The  geographical  locations  of  the  tests  were 
determined  to  within  1/4  of  a  degree.  The  data  point  times  were 
converted  from  local  time  to  decimal  GMT  time.  The  original 
values  for  the  relative  humidity,  wind  speed,  cloud  properties 
and  reference  temperature  were  then  converted  to  units  acceptable 
to  PROTDRB.  Table  2  and  3  are  samples  of  the  results  of  that 
effort.  Table  2  is  the  time  dependent  inputs  and  Table  3  is  an 
example  of  the  PROTURB  inputs  with  an  identification  of  the 
various  input. 


Time 

Temperature 

Wind  Speed  Relative 

Cloud 

Cloud 

Cloud 

Ceiling 

GMT  degrees  Kelvin 

m/sec 

Humidity 

Cover 

Type 

Density 

Height 

2000 

289.1500 

6 

0.94 

99 

0 

4 

1180 

2025 

288.9500 

.7 

0.94 

99 

0 

4 

900 

2050 

288.7500 

1.1 

0.94 

99 

0 

4 

700 

2075 

288.7500 

.3 

0.95 

99 

0 

4 

650 

2100 

288.7500 

.3 

0.95 

99 

0 

4 

200 

TABLE 

2  AN  EXAMPLE  OF  METEOROLOGICAL 

DATA  PREPARED 

1  FOR  PROTURB 

1993  9 

8 

35.8  -121.2 

2100.0 

20 

1 

600  285.00 

310.25 

973,0  2.0 

0.416 

2.0 

0.05 

0.26  50 

3.0 

10.0  1000.0 

1.06 

2 

4 

0.92  3 

2 

3  0 

284.65 

10.0 

4.0 

0.0  32.0 

0 

0 

LINE  1 

LIME  2 

LINE  3 

LINE  4 

LINE  5 

YEAR 

CLOUD  HGT 

ALBEDO 

GROUND 

CUT  OFF  FREQ 

EMISSIVITY 

MONTH 

REFERENCE 

INTEGRATION 

SOIL 

FREQ  CALC 

TEMPERATURE 

SUB-INTERVAL 

MOISTURE 

INTERVAL 

DAY 

EST.MAX 

RECEIVER 

SOIL 

DEBUGGING 

TEMPERATURE 

HEIGHT 

COMPOSITION 

FLAG 

LATITUDE 

REFERENCE 

SOURCE 

VEGETATION 

PLANETARY 

PRESSURE 

HEIGHT 

FACTOR 

EXTENTION  FLAG 

LONGITUDE 

WIND  SPEED 

PATH 

SNOW 

LENGTH 

TYPE 

GMT 

RELITIVE 

WAVE 

EST.  MIN 

HUMIDITY 

LENGTH 

TEMPERATURE 

%  CLOUD 

REFERENCE 

CLOUD 

FOCAL 

COVER 

HEIGHT 

TYPE 

LENGTH 

GROUND 

ROUGHNESS 

CLOUD 

ENTRANCE 

COVER 

ELEMENT 

DENSITY 

PUPIL  DIA. 

TABLE  3  AN  EXAMPLE  OF  A  PROTURB  INPUT  DATA  SET 

The  unstable  atmosphere  diurnal  output  from  the  PROTURB  model 
was  then  compared  to  the  corresponding  measured  values  for  Cn2 
and  the  incoming  solar  radiation.  Error  analysis  was  performed 
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on  the  unstable  atmospheric  condition  portion  of  the  numerical 
results. 

Incomplete  knowledge  of  the  test  site  forced  the  assumption  of 
some  of  the  values  for  the  air  pressure,  the  soil  characteris¬ 
tics,  and  the  vegetation  characteristics.  These  values  were 
estimated  through  the  use  of  site  photographs,  topographic  maps 
and  data  from  surrounding  sites.  The  most  critical  of  these 
inputs  is  the  estimation  of  soil  wetness  which  was  estimated 
through  the  use  of  rainfall  data  and  on-site  observations.  The 
beam  path  was  considered  to  be  a  constant  height  of  2.0  meters 
above  the  ground.  The  remaining  parameters  are  presented  above 
in  the  sample  input  cards. 


5.  RESULTS 


The  results  of  the  comparisons  are  presented  in  Figures  3 
through  5  which  show  both  the  and  the  incident  shortwave 
radiation.  Figure  3  shows  comparisons  from  Ft.  AP  Hill.  Agreement 
in  solar  radiation  is  very  good.  The  C_^  does  show  some  large 
differences  but  are  still  considered  goc^  when  one  considers  the 
terrain  of  the  path.  The  comparisons  do  show  a  vast  improvement 
over  the  predecessor  model  ”IMTDRB''  as  discussed  by  Smith  and 
Corbin^.  Figxire  4  shows  Ft.  Hunter-Liggett  data  which  also  looks 
quite  good.  The  deviation  shown  on  6  September  is  probably  due  to 
differing  cloud  conditions  between  the  base  site  and  the  observa¬ 
tion  site  3  kilometers  away.  Figure  5  shows  Eglin  AFB  data.  Here 
the  cloud  observation  was  6  km  from  the  data  site.  Again  the 
agreement  is  considered  very  encouraging.  If  perfectly  represen¬ 
tative  data  were  available,  better  agreement  would  be  expected. 

The  daylight  portion  of  the  error  analysis  is  presented 
in  Tables  4  and  5.  These  results  include  only  times  when  the  at¬ 
mosphere  was  unstable.  The  percent  variation  was  determined 
by  the  relation  %Variation  =  100.0*( (Measured-PROTURB) / 
Measured) .  The  mean,  standard  variation,  standard  error  and 
skewness  of  the  percent  variations  were  then  determined.  The 
radiation  difference  is  presented  in  units  of  watts/meter^  and  is 
a  straight  difference  between  the  measured  and  calculated  values. 

The  values  for  percent  variation  of  C_^  at  Ft.  A.  P.  Hill 
were  quite  large  and  negative.  The  negative  values  are  an  indi¬ 
cation  that  the  PROTURB  model  predicted  larger  C-^  values  than 
were  measured.  The  negative  values  for  the  skewness  also  indi¬ 
cates  the  presence  of  a  systematic  error.  It  is  believed  that 
the  main  source  of  error  is  in  the  estimation  of  a  2.0  meter 
height  for  the  scintillometer  positions.  The  actual  site  pos¬ 
sessed  several  shallow  canyons  along  the  propagation  path.  This 
would  mean  that  the  actual  average  path  height  was  greater  than 
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TABLE  4 .  Variation 

statistics.  Data  points 
that  indicate  stable  or 
neutral  conditions  have 
been  excluded  from  the 
statistical  analysis. 
Also,  experimental  data 
that  indicate  equipment 
malfunction  have  been 
excluded. 
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1,37096 

TABLE  5.  Short  wave 
radiation  difference 
statistics. 


2.0  meters.  An  increase  in  the  assumed  path  height  would  reduce 
the  Cn2  values  generated  by  PROTURB  thus  reducing  the  error.  The 
occurrence  of  fog  and  rain  along  the  path  should  also  be  consid¬ 
ered  a  source  of  error  since  PROTURB  does  not  include  the  actions 
of  such  atmospheric  additives.  In  general,  the  PROTURB  values 
tend  to  possess  the  same  trends  as  the  measured  values  and  though 
the  error  is  quite  high  the  results  are  comparable  with  that  of 
models  with  similar  input  simplicity. 

The  values  for  percent  variation  of  at  Ft.  Hunter- 

Liggett  were  very  reasonable.  There  appears  to  a  sea  breeze 
effect  present  in  the  measured  data  in  both  the  morning  and 
afternoon.  The  PROTURB  model  assumes  that  the  atmospheric  condi¬ 
tion  is  generated  by  the  surface-air  interface  and  not  influenced 
by  outside  forces  such  as  weather  fronts  and  nearby  surface 
inhomogeneities  such  as  a  land-sea  interface.  Thus  the  prevailing 
deviations  were  seen  in  the  morning  and  early  evening  were  less 
pronounced  possibly  due  to  the  effect  of  the  nearby  ocean.  The 
surface  at  Ft.  Hunter-Liggett  was  relatively  flat  leading  to 
better  overall  results. 

The  values  for  percent  variation  of  at  Eglin  AFB, 

were  moderately  successful.  The  PROTURB  model  was  more  sensitive 
to  wind  speed  and  cloud  coverage  than  the  measured  data.  In 
general,  PROTURB  would  have  predicted  larger  values  for  results 
had  it  not  been  for  this  sensitivity.  The  over  prediction  of  C_ 
can  be  attributed  to  PROTURB 'S  lack  of  sensitivity  in  relation  to 
soil  type  and  soil  moisture  effects.  An  upgraded  soil  model  is 
currently  under  development  and  should  be  installed  by  the  time 
of  public  release. 

The  shortwave  radiation  comparison  at  all  sites  were  very 
good.  The  primary  deviation  is  in  the  shape  of  the  clear-sky 
curve  directly  after  sun  rise  and  before  sunset.  This  greater 
measured  radiation  at  these  points  can  be  assumed  to  be  related 
to  the  value  of  the  diffuse  sky  radiance  at  low  solar  elevation 
angles.  It  appears  that  the  PROTURB  cloud  reduction  factors 
represent  reasonable  approximations  to  the  effects  of  clouds. 

A  more  complete  analysis  is  being  prepared  for  the  use  of 
the  Army  Research  Laboratory  by  Cranston  and  Cruncleton. 


CONCLUSIONS 

Given  that  some  necessary  approximations  were  made  in 
transforming  from  a  nonuniform  terrain  to  a  uniform  terrain  with 
fixed  site  parameters,  the  PROTURB  model  did  a  reasonable  job  in 
the  determination  of  Cj^^.  It  was  noticed  that  the  significant 
decreases  and  increases  in  ,  due  to  clouds  and  wind  effects. 
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are  parallel  in  both  PROTURB  and  measured  values,  but  that  the 
PROTURB  model  appears  to  be  overly  sensitive  to  these  factors. 
The  primary  factors  responsible  for  this  sensitivity  are^  the 
Angus~*Leppan  wetness  and  cloud  parameters,  and  the  thermal  diffu— 
sivity  of  the  soil.  These  factors  play  important  roles  in  the 
energy  partition  within  the  ener^  balance  model  and  should  be 
scrutinized  further.  The  comparisons  presented  here  considered 
only  the  unstable  condition.  Noting  the  collected  turbulence  data 

shown  on  figure  2,  it  is  recommended  that  an  effort  be  made  to 
extend  the  model  into  night  or  stable  conditions. 
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These  turbxiles  modify  the  cloud 


microphysical  particles  to  such  an  extent  that 


this  is  the  primary  controlling  factor. 
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SH ARC-3:  A  Model  for  Infrared 
Atmospheric  Radiance  at  High 
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ILLUSTRATIVE  SHARC  CALCULATION 
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RADIATIVE  EXCITATION  MODULE 
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COMPLETE  ROTATIONAL  LEVEL  REDISTRIBUTION 
NO  LINE  OVERLAP 
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GEOMETRY  MODEL  INSURES  THAT  LOS  TRAJECTORIES 
INTERSECT  AURORA  AS  DESIRED 
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DEC  VAX  8800  :  4.5 

DEC  VAX  780  :  33.4 
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2  System  of  vibrational  energy  levels  and  radiative 
transitions  of  CO2 

Wc  considerBd  (ill  vibrational  levels  and  (til  ro-vibrational  lines  for  tiie 
isotope  wbicli  are  available  from  HITRAN-92i  tbe  excitation  of 
1/2“  and  V2r  vibrations  up  to  5i/2  and  81/3  respectively. 

Radiative  processes:  (a)  absorption  of  direct  solar  radiation  and  the  ra¬ 
diation  coming  from  the  lower  boundary,  (b)  absorption  of  atmospheric 
radiation  and  (c)  spontaneous  and  stimulated  emission. 

Quantitative  characteristics  of  the  model: 

Total  number  of  vibrational  levels:  L„=127 

Toted  number  of  ro- vibrational  levels:  T=7499 

Total  number  of  bands:  324 

Total  number  of  ro- vibrational  lines:  29993 
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2.1  Collisional  processes 


2.1.1  R-T  processes 

•Assumption  of  rotationeil  LTE  (infinite  large  rate  of  the  R-T  exchange) 

2.1.2  Intra-molecular  V-V  processes 

C02in^V3)  4-  iV/(iV2, 02)  C02{n'‘'{v3  -  1))  +  M  ^  AE, 

where  n'  =  n  -f  1,  n  -I-  2,  n  -|-  3,  n  4. 

C02in‘v2,m)  -h  M{N2,  O2, 0)  -*  C02{'n}' vz,Tn')  -H  M  -f  AE, 

where  m  and  m'  denote  arbitrary  levels  from  the  groups  n^vz  and  n^'vz  respectively. 


2.1.3  Inter-molecular  V-V  processes 


C02{n%)  -f  '¥2(0)  ->  C02[n^\vz  -  1))  +  Nzil)  AE 
COzin^vz)  T"  ^^2(0)  — ^  (702((n  “t"  1)^  (us  —  1))  “b  02(1)  +  AE 
C 02{n^Vz)  "b  ^^2(1)  — ^  OOzK'fi  +  2)^  U3)  -r  02(0)  -h  AE 

-h  (00*^0)  (n''(t;3  -  1))  +  (00°1) 

•^¥2(1-)  "b  ^^2(0)  — ^  -¥2(0)  -b  02(1)  “b  AE 


2.1.4  V-T  processes 


C02(n‘vz)  +  iV/(iV2, 02)  COzdn  -  1)^3)  -h  M  +  AE,  |/  -  l'\  =  1 
C02(t1’^Vz)  “b  0  — *■  OOzK^^  —  1)^  V3)  -b  0  -b  AE 
C02(n‘vz)  -b  0  ->  C02(n‘'(vz  -l))  +  0  +  AE 
Nzil)  +  0 NziO)  +  AE 
02(1)  "b  M(N2,  O2, 0)  — *  02(0)  “b  M  "b  AE 


2.2  Reference  atmosphere 


Temperature  profiles,  mixing  ratio  profiles  for  N2,  O2  and  0:  MSISE-90 
reference  atmosphere. 

The  volume  mixing  ratio  for  CO2  was  compiled  on  the  basis  of  experi¬ 
mental  data,  presented  by  Trinks  and  Fricke  (1978).  Offermann  et  al  (1981) 
and  Rincland  et  al  (1992). 


2.Z  Computational  features 

On  each  iteration  step  we  solve  the  radiative  transfer  equation  in  each  of 
29993  ro- vibrational  lines  using  the  Feautrier  technique  for  16  frequencv 
and  4  angle  grid  points. 

In  order  to  reach  the  convergence  criterion  10“"^  for  vibrational  popula¬ 
tions^  6  and  40  accelerated  lambda  iterations  are  required  for  day  and  night 
conditions  respectively. 
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2.4  Results 


2.4.1  Verificacion  of  standard  approximations 

1 .  Effect  of  highly  excited  levels 

2. ” Emission”  and  ’’radiative  temperature”  approximations  for  week  bands 

3.  Spontaneous  emission 

^  ^  '^vj-^vjyy  ~  j 

j  j 

where 

r^r>  9v(j)^^Pi~%^) 

lyLTEf-  /TIN  _  ksT  J 

"  ^  Q4T) 
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-  Vibrational  temperatures  of  CO2  bend-streteh  states 

-  Atmospheric  infrared  cooling  rates 
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lly  variable  emission  and  absorption  lineshapes 

Dependence  on  temperature,  pressure 


Modified  Curtis  Matrix  (MCM)  Model 
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ineshape,  temperature,  density,  pressure 


Correlating  the  Models 
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Treatment  of  lower  boundary 


Model  Atmospheres 


Kinetic  Temperature  (K) 


Temperature  Difference 


Summary  of  Results,  Vibrational  Temperatures 
‘Fast”  case  (realistic  kinetic  parameters) 

-  Below  70  km,  excellent  agreement 
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MCM  consistently  predicts 
lower  vibrational  temperatures  in  the  thermosphere 
higher  vibrational  temperatures  at  the  mesopause 
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Temperature  Difference  (%) 
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-  Result:  MCM  calculates  more  total  radiative  excitation  at  mesopause 

Residual  discrepancies  explained  by  MCM  path-averages 


rate  (deg/day) 


Summary  of  Results,  Cooling 
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Peak  differences  occur  in  narrow  altitude  ranges 


Conclusions 

The  two  radiative  transfer  algorithms  give  consistent 
results  for  populations  and  cooling 
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from  non-LTE  atmospheric  emissions 

Cooling  rates  are  accurate  enough  for  use  in  dynamical 
models  and  thermal-balance  calculations 
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WAVENUMBER  (cm 


Spectral  Fits 
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TVend  of  population  ratio  is  consistent  throughout  data  set 

-  Much  less  variation  seen  than  in  radiance  or  temperature 

-  Possible  effect  of  chemical  production  seen  in  auroral  data 


Nitric  Oxide  Energy  Level 

Diagram 
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Vibrational  Level 


Daytime  and  Nighttime  NO  Spectra 
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WAVENUMBER  (cm 


Spin-Orbit  Altitude  Profile 
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£2-  /j  FRACTION  TEMPERATURE  (K) 

Expectation:  open  symbols  MSIS  Temperature  Profile:  line 

Observation:  filled  symbols  Observation:  symbols 


Origin  of  Subthermal  Spin-Orbit 
Population  Distribution 
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Chemical  production  has  diurnal  variation  not 
observed  in  sublevel  distributions 


Impact  On  Related  Research 
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production  from  collisional  uppumping  in  NO 


Comparison  of  Global  Variations  of  the  Radiance  of 
the  Stratospheric  O3,  CH4  and  HNO3  Spectra  as  Viewed  by 
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CLAES  -  CIRRIS  lA  OVERVIEW  FOR  CORRELATIVE  STUDY 
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CLAES  DESIGN  AND  FUNCTIONAL  DIAGRAM 
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STARE  WODE  VERTICAl.  SCAN  HO  R  1 2.  R  AS  Tf;  ^  STAiRCASE 

MOOE  MODE  SCAN  MCCt 


CIRRIS  lA  viewing  &  illustration  of  gimbal  motion  for  typical  mode 

[Bartschi  B.  et  al,  1987] 


Typical  course  of  the  Shuttle  during  STS-39  flight  for  April  28  and  May  6, 
1991.  Photograph  shows  CIRRIS  lA  orbit  and  typical  look  angles  and 
projected  auroral  intensity  and  locations.  A  portion  of  the  CIRRIS  lA  flight 
was  dedicated  to  stratospheric  gas  measurements.  Courtesy  of  USAF  Phillips 
Laboratory  and  Radex,  INC. 


LATITUDE 


MEASUREMENT  CYCLE 


CLAES  viewing  geometry  [Roche  et  al,  1993] 


CLAES  NORTH-SOUTH  LOOKING  COVERAGE 
DURING  STRATOSPHERIC  GAS  MEASUREMENTS 


-150  -120  -90  -60  -30  "g  30  60  90  120  150  1S0 

LONGITUDE 


CLAES  latitude  versus  longitude  coverage  for  the  month  of  April  1992. 
Points  at  which  calibrations  were  performed  for  each  descending  day  are 
shown.  Violet  indicates  north  looking;  acqua  indicates  south  looking;  blue  is 
the  overlap  regions.  Upper  and  lower  panels  each  contain  about  600  data 
points  from  16  orbits  over  a  24  hour  period  covering  any  latitude  band. 
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Typical  course  of  the  Shuttle  during  STS-39  flight  for  April  28  and  May  6, 
1991.  A  portion  of  the  CIRRIS  lA  flight  was  dedicated  to  stratospheric  gas 
measurements.  Points  indicate  seven  data  sampling  profiles  locations. 


CALCULATED  RADIANCE  AT  THE 
SENSOR  FOR  LIMB  TANGENT 
ALTITUDE  h,  FREQUENCY vi 


'’^LoMieea 


Ri(h,v,)  =  -  J  J  I  B(vi.T)  F(^)  R(v,)  d«I^S 


Vi  s»0 


PLANCK  FUNCTION  AT  v,,  T,  POINT  S 


INSTRUMENT  SPATIAL  RESPONSE  FUNCTION 


INSTRUMENT  SPECTRAL  RESPONSE  FUNCTION  AT 


DERIVATIVE  OF  TRANSMITTANCE  FROM  POINT  S  TO  SENSOR 


WHERE 


Tj  =  exp 


-|  Ml*. 


T)  <1  (Is 


AND 


Kv(l*.T)  = 


:r[(v-Vor  +  a^] 


SPECTRAL  ABSORPTION 
CROSS  SECTION  FOR 
LORENTZ  LINE  SHAPE 


THE  MEASURED  RADIANCE  IS  OBTAINED  AS 

Ni(h,  vi)  =  D|  (VOLTS)  x  CAL|  (RADIANCE/VOLTS) 

BY  COMPARING  ALL  R|  AND  D|  FOR  EACH  TANGENT  ALTITUDE,  SOLUTIONS 
ARE  FOUND  BY  ITERATING  T,  P,  AND  q 


X00907  R  JS  LC 


X90907  R  JS  LC 
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CLAES  ETALON  AND  BLOCKER  FILTER 


NOISSMSNVUL 
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WAVENUMBER  (cm  *) 

Example  of  CLAES  blocker  etalon  transmittances  for  Etalon 
tilt  positions  for  blocker  filter  region  4. 


CIRRIS  lA  INFRARED  AND  ATMOSPHERIC  RADIATION  MODEL 
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The  spectral  distribution  of  infrared  emission  produced  by  minor  species  in  the  stratosphere. 
The  figure  compares  the  AFGL  radiance  model  (FASCODE3P  line-by-line  radiation  code), 
with  data  taken  by  the  Cryogenic  Infrared  Radiance  Instrumentation  for  Shuttle  (CIRRIS  1  A) 
viewing  the  earth  limb  at  altitudes  of  31  and  35  km. 


CLN03^\  :  / 
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R  profiles  used  m  the  analysis  of  mid-latitude  spectra  over  the  southern  hemisphere 


INFRARED  ATMOSPHERIC  RADIATION 
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The  spectral  distribution  of  infrared  emission  produced  by  minor  species  in  the  stratosphere. 
The  figure  shows  the  spectral  emissions  measured  by  the  Cyrogenic  Array  Etalon  Spectrometer 
(CLAES)  viewing  the  earth  limb  at  altitude  of  35  km. 


INFRARED  ATMOSPHERIC  RADIATION 
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CLAES,  CIRRIS  lA  DATA  AND  RADIATION  MODEL  COMPARISONS 


(rao  ,.js  j.mo  SHVAO  HDNViaVH  TVHlDHdS  SMI 
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WAVENUMBER  (cm  ‘) 

Correlation  and  anti-correlation  between  CLAES,  CIRRIS  lA  03  data  compared  with  FASCODE3P  model. 
The  figure  shows  that  the  line  position,  intensity  and  spectral  distributions  measured  by  CLAES  and  CIRRIS 
1 A  agree  with  each  other  and  with  those  of  the  model  at  this  altitude. 
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CLAES  and  CIRRIS  lA  HN03  spectra  compared  with  model  at  23  Km  tangent  height.  The  figure  shows 
the  measured  and  model  spectra  agree  quite  well.  The  model  spectral  resolution  is  0.02  cm'*.  When 
intercomparing  the  residual  spectra,  the  slightly  different  scaling  should  be  regarded. 


CLAES,  CIRRIS  lA  DATA  AND  RADIATION  MODEL  COMPARISONS 


(mo  i.iSj.mo  SIIVAO  HDNViaVH  TVULOHdS  Smi. 
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Correlation  between  CLAES,  CIRRIS  lA  CH^N^O  data  compared  with  line-by-line  model  radiation 
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ic  % 


Q^Z^y  692,5 


692.7  , 

6?2..f  A/a\ 


Fig.  1.  A  small  ^tion  of  solar  central  intensity  spectrum  at  692.5  nm  plotted  at  full  scale  and  at 
10  times  s^e.  The  heapr  lines  are  the  observed  FTS  spectrum  from  James  Brault  at  Kitt  Peak. 
The  resolution  is  522000  and  the  signal-to-noise  is  about  3000  which  is  poor  for  this  work.  The 
continu^  levd  is  imcertain.  The  thin  line  is  the  computed  spectrum.  It  would  be  possible  to 
adjust  the  line  data  for  a  better  match  but  most  of  the  discrepencies  are  from  missing  lines.  There 
are  solar  Imes  of  Ca  I  Ti  I,  Cr  I,  Fe  I,  C  I,  CN,  and  terrestrial  lines  of  H2O  and  The  first 

numbCT  m  ea^  Ime  labd  is  the  last  3  digits  of  the  wavelength  and  the  4th  number  is  the  per  mil 
central  intensity  if  the  line  were  computed  in  isolation. 


Annual  Review  Conference  on  Atmospheric  Transmission  Models 
7-8  June  1994,  Phillips  Laboratory,  Bedford,  MA 
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POAM  II  Additional  Science  Team  Members 
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Polar  Ozone  and  Aerosol  Measurement  (POAM 11)  Science 
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POAM  II  measurements  are  vital  to  understanding  the  Antarctic 
''ozone  hole"  and  possible  development  of  Arctic  ozone 
depletions. 


POAM  II  -  (Polar  Ozone  &  Aerosol  Measurement) 
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OPTICAL  HEAD  ASSEMBLY 
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Figure  1:  Schematic  drawing  showing  the  location  of  the  POAM  II  instrument  components  integrated 
onto  the  SPOT  3  satellite  (left)  and  schematic  of  the  POAM  II  Optical  Head  Assembly  (right). 


SCIENCE  SENSORS  OPTICAL  PATH 
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Figure  2:  POAM  11  science  sensors  optical  path  arrangement. 


Atmospheric  Extinction  vs  Wavelength  at  18  km 
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Altitude  (km)  Altitude  (km) 
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Altllude  (km)  Altitude  (km) 
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Deshler  &  Johnson  McMurdo  profile  Oct  15  compared  with  POAM II  Oct  16 

^ - 1 - 1 - 1 - 1 - 1 - 1 - 1 - \ - 1 - 1 - 1 - 1 - 1 - 1 - 1  T  I  I  I  I  I  I  I  I  I  '  '  I  I  ^ 


(N 


(UD[)  sptupTV 


663 


5.0*10^^  1.0*10"^  1.5*10"^  2.0*10"^  2.5*10"^  3.0*10 

O3  Density  (cm‘^) 


to 

to  Q 

CO 

(ODI)  epnwiV 

to 

CM 

8 

665 

Polar  Vortex  6c  POAM  Events:  931108 
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Large-Scale  Retrieval  of 
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California  Institute  of  Technology 
Pasadena,  California  91109 


674 


uojidjosqv 


Wave  number 


SEASCRAPE 


Q 

1 


9 

•«  ^ 
S  3 
«  6 


•»«< 


k 

•U 

e 

d 

k 

d 


§ 

5  c/) 

§  Kl 

d 
d 

•w 

s 
g 

®  . 
U  Oh 

<£  « 
^•c 

6  2 

•3  -S 

‘C  ^ 

|>  6 

^  -w 

e 

O  O 


d 

3 


I 


d 

Q) 

•»ij 

k 

•w 

Z  « 

H  'U 

hW  ^ 

3 

3 
d<  o 


3 


Q) 

CO 


d 

5 


c3 


cd 


(D 


73 


\,v 

P.  tiH 


CO 


T3 

d 

cd 

CO 

in 

/rH 

HH 

O 

;-( 

P. 

o 

'd 

<u 

a 

CO 

o 


p 

in 

'a 

•rH 

4-J 

3 

6 


cfl  (0 

S  B 

‘d  3 

6 
p 
u 
cd 
P. 


(U 

;h 

8 

d 

a  S 

3 

e  CO 

!»  S 


a 

V 


u 

o 

CO 

CO 

<u 

o 

o 

u 

p. 

73 

3 

Vh 

cd 

Ph 

cd 

d 

o 

CO 

CO 

cd 

§ 

cd 

-f-» 

cd 

nD 

Hi) 

d 


3 

O 

d 

6 


675 


Retrieval  algorithm  integrates  the  forward  model  with  the 
inversion  algorithm  (the  forward  model  can  stand  alone  as  a 
simulation  tool) . 


Features  of  SBASCRAPE 
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Constraints, and  a  priori: 

A.  Hydrostatic  equilibrium. 

B.  Parameter  ranges. 


qfSEASCRAPE  {cont.) 
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Clouds,  aerosols,  and  scattering,  possibility  incorporating  existing 
electromagnetic  scattering  codes. 

Retrieval  of  instrumental  parameters. 
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Utmties 


SEASCRAPE 


updaied  model  data 


JPL 
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=  frequency; 
=  path  length; 


GENLN2  vs.  SEASCRAPE  -  Case  2aa 
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Least  Squares  Process 
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Optimal  Estimate  of 
Atmospheric  Profile 


Execution  Time  vs.  Number  of  Partiais 
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correlatipns  between  the  errors  in  the  retrieved  parameters; 
the  optimal  parameterization  of  the  observational  model. 


JPL  Parallel  Implementation  of  SEASCRAPE 
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-  MIMD  (each  processor. must  be  independent  to  avoid  serious  load  imbalance) 

-  high  performance  (10-30  MFLOPS) 

-  large  DRAM  (16-64  megabytes) 
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ATMOSPHERIC  AND  ENVIRONMENTAL  RESEARCH,  INC. 

Annual  Review  Conference  on  Transmission  Models 
Phillips  Laboratory,  Hanscom  AFB,  MA 
June  7-8, 1994 
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•  DIfIcult/expensIve  to  vary  Instrument  characteristics  and  to  revise  spectroscopic 
parameters  and  atmospheric  representation 

Answer: 

•  Need  to  sacrifice  actual  execution  time  for  enhanced  flexibility  and  versatility 
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Block  Diagram  of  Retrieval  Algorithm 
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K(x  j)  =  Jacobian  or  sensitivity  matrix 
r  =  matrix  of  "damping"  coefficients 


Optimal  Retrieval  Approach 
(see  Rodgers,  1 976  and  Eyre,  1 989) 
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r  can  be  manually  tuned  or  evaluated  from  background  error  statistics  when 
available  (maximum  probabiiity  soiution  (see  Rodgers,  1976) 


Radiative  Transfer  Aigorithm 
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Derivative  calculations  for  temperature  (at  23  levels)  and  water  vapor  (18  levels)  - 
continuous  spectral  coverage  from  500  to  1400  cm"^  at  full  spectral  resolution:  30 
min.  with  old  version  on  Convex  at  Phillips  Lab.(Goal  with  improved  version:  <10 
min). 
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RETRIEVAL  OF  HIGH  ALTITUDE  OH  PROFILES  FROM  BALLOON  BORNE 
MICROWAVE  SOUNDER  (Steve  Miller,  PL/OS)  ■-•-'Jui'i  Durint 
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Future  Work 
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Long  term:  treatment  of  clouds  and  aerosol  by  coupling  the  code  with  CHARTS 
(Moncet  and  Clough, 1992) 


Acknowledgements 
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SIMULATION  OF  STELLAR  OCCULTATION 
MEASUREMENTS 


Liisa  Oikarinen 

Erkki  Kyrola,  Elina  Sihvola,  and  Johanna  Tamminen 
Finnish  Meteorological  Institute 
Department  of  Geophysics 


1.  Stellar  occultation  method 

2.  GOMOS  instrument 

3.  Atmospheric  model 

4.  Data  inversion 


GOMOS  =  Global  Ozone  Monitoring  by  Occultation  of  Stars 

-  An  ESA  Developed  Instrument  (EDI),  funded  under 
the  ENVISAT-1  programme 

-  Main  objective  of  GOMOS:  trend  studies  of  ozone 

-  Measurement:  global  mapping  of  stratospheric  O3,  NO2, 
NO3,  aerosols,  H2O,  temperature  (O2) 

-  Launch  1998 


Role  ofFMI  in  GOMOS: 

-  Scientific  help  in  instrument  development 

-  Development  of  data  retrieval  methods 

-  Data  center 
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The  stellar  occultation  measurement  principle 


O 


STELLAR 

SPECTRUM 


ENVISAT-1 

I 

GOMOS 


Tangent  height  15  - 100  km 

Vertical  resolution  2  km  (CCD  integration  time  0.5  s) 
Stars  up  to  rnagnimde  4  are  used  (bright  limb  up  to  2) 
Heliosynchronous  polar  orbit,  altimde  8(X)  km 
Nominal  orbit  period  100  min 
25  -  40  occultations  per  orbit 
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of  GOMOS  measurements 


Longitude 


GOMOS  INSTRUMENT 


UVIS-spectrometer: 

-  Large  spectral  width:  250  -  675  nm, 

both  Chappuis  and  UV-bands  can  be  exploited: 
ozone  can  be  retrieved  at  15  -  90  km 

-  A  good  spectral  resolution  (1  nm),  i.e.,  many  pixels: 

Counteracts  faintness  of  stars; 

Several  trace  gases  can  be  retrieved:  O3,  NO2,  NO3,... 

-  Two  dimensionality  of  CCD: 

Solar  scattered  backgroung  can  be  measured, 
daytime  measurements  feasible; 

Pointing  errors  counteracted; 

Refractive  bending  counteracted 

IR-spectrometer: 

-  O2  and  H2O  measurement  (756  -  773  and  926  -  952  nm, 
0.06  nm  resolution) 

-  From  O2  neutral  density  and  temperature  to  support 
UVIS-analysis 

Photometers: 

-  470  -  520  nm  and  650  -  700  nm,  maximum  sampling 
rate  1  kHz 

-  Measurement  of  scintillations:  The  stellar  signal  can  be 
corrected 

-  Studies  of  turbulence 
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Why  stellar  occultation? 


Calibration  stability 
Global  coverage 
Altitude  resolution 


Solar  occultation  -  SAGE 

+  relative  measurement 

-  only  two  occultations  per  orbit 

-  night/day  terminator  measurements 

+  good  altitude  resolution  (1  km) 


SBUV-TOMS,  TOYS 

-  absolute  measurement 

+  good  coverage 

-  only  dayside  measurements 

-  poor  altitude  resolution  (8  km) 


Stellar  occultation-GOMOS 

+  relative  measurement 
+  global  coverage  (night  and  day) 

+  sufficent  altitude  resolution  (2  km) 
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ATMOSPHERIC  PHYSICS  FOR  GOMOS 


Phenomena  that  should  be  simulated: 


-  Absorption  by  O3,  NO2,  NO3,  OCIO,  SO2, 

BrO,  H2O,  O2 

-  Scattering  out  of  the  beam  by  Rayleigh  and 
aerosol  scattering 

-  Scattering  of  solar  light  by  molecules  and  aerosols 
into  the  FOV  of  the  instrument 

-  Emission  into  the  FOV  (aurora  etc.) 

-  Refraction:  bending  of  the  ray,  refractive  attenuation 

-  Turbulence:  scintillations 
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;  T=133iOK 

111=1.35  ■ 

(5000  —  iTangeiit  heiglii=30  km 


Wavelength  (A) 


SPECIAL  REQUIREMENTS  OF  GOMOS 
ATMOSPHERIC  SIMULATION  MODEL 


-  compactness  (analytical  density  distributions  etc,) 

-  the  code  should  be  easy  to  modify 

“  realistic  variations  in  constituent  distributions 
(ultimate  conditions,  eg.  the  ozone  hole,  required  in 
addition  to  average  atmospheres) 

-  occultation  geometry 


Comparison  of  FINGOM  (a  code  specially  designed 
for  GOMOS  by  EMI)  and  LOWTRAN 


model 

FINGOM 

LOWTRAN 

absorbers  in  the  vis. 

O3,  NO2,  NO3,  SO2, 
BrO,  OCIO 

O3 

constituent  profiles 

3-D,  seasonal,  diurnal, 
latitudinal  variations 

local  geographical- 
seasonal  atmospheres 

atmospheric  model 
description 

analytic 

tables 

refractive  attenuation 

included 

not  included 

scattered  solar 
background 

single  scattering 

single-  and  multiple 
scattering 

aerosol  model 

simple 

models  for  various 
aerosol  conditions 

•  toial  iransmission 

•  Rayleigh  scarieHng 
-  aerosof  scaaering 

■  O3  absorption 

’  NO2  absorption  ' 

■  refractive  attenuation  i 


500 

Wavelength  (nm) 


Atmospheric  transmission  for  a  ray  path  of  tangent 
height  20  km.  Total  transmission  is  divided  into  the 
main  extinctive  components  Ttot  -  To3TRaylTaerTN02Trefr 


LaliUicIo  • 


Dobson  units 


Total  ozone  distribution  model  as  a  function  of 
latitude  and  season. 
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Transmission  . 
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INVERSION 


T(X,l)  =  lkiX,s)ds  =  jXPiis)cF.iX,Tis))ds  ,  . 

...  ■  .....  /  .  .  i  i .  ./■.■■■ 

Discretizafion: 

0  l-  x{X,l)  =  ^N,(l)a,(X,TM 

The  same  temperature  T{1)  is  used  for  the  whole  ray 
.........  .  .2®  Niil)  =  X PiUj.  )^j  ’  =hne  density  of  constituent  i 

t(a,/)=xsp^  :' 

.  i  jel 

.  .  1°  ■z(X,l)  =  I,k(X,Sj)Asj 

:  .,  2°  k(X,Sj)  =  'Zp^(Sj)G,{Xj^^^  ;  7/ 
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80 


40  ;  60 

Tangent  height  (km) 


ated  line  density  for  ozone,  N02,  and  N03 
inction  of  the  tangent  heigth  of  the  path. 


Altitude,  (km)  , 

The  standard  deviation  in  vertical 
inversion  (upper  line).  The  lower  line 
is  the  std  in  the  data. 
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CONCLUSIONS 


An  atmospheric  model  with  analytical  density  distributions  has 
been  constructed  at  FML 

Realistic  simulations  of  occultation  measurements  are  possible 
with  the  model. 


Problems: 


-3-D  time  dependent  models  are  difficult  to  construct 
from  the  existing  data 
(exception:  O3  and  temperature) 

-  Absorber  cross-sections  are  required  at  160-280  K, 
not  available  at  the  moment 


Open  points: 

-  simulation  of  scintillations 

-  simulation  of  auroral  emissions 
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Overview/Objectives 


Line-by-line  Model 

•  LBLRTM,  (Line  By  Line  Radiative  Transfer  Model) 

•  Vectorizable,  10-20x  faster  (based  on  FASCODE) 

•  Extensively  validated  against  observed  spectral  radiance  data  (HIS) 

Water  Vapor,  Carbon  Dioxide,  and  Ozone: 

•  Mid-latitude  summer,  Mid-latitude  winter,  Tropical  atmospheres 

-  Calculate  line-by-line  upwelling,  downwelling,  and  net  (up-down) 
clear-sky  fluxes  at  5  cm-1  resolution  for  the  full  IR  spectrum 
(10-3000  cm-1) 

-  Cooling  rates  computed  from  the  divergence  of  the  net  flux 

-  Current  US  standard  profiles  used  for  H2O  and  O3;  CO2  =  355ppm 

•  Doubled  carbon  dioxide  considered  for  CO2  =  710  ppm 

CH4,  N2O,  and  the  Halocarbons 

•  Halocarbons  considered  include  CFC-11,  CFC-12,  CFC-22,  and’CCL4 

•  Profiles  for  Methane,  Nitrous  Oxide,  and  Halocarbons  derived  from 

a  2D  chemistry  model 

1990-2000  Changes 

•  Ten-year  profile  changes  in  all  trace  species  considered  for  their 

impact  on  fluxes  and  cooling  rates 

•  Changes  in  water  vapor  have  not  been  considered 

Rapid  Radiative  Transfer  Model 

•  Apply  LBL  flux  and  cooling  rate  results  to  the  selection  of 

appropriate  spectral  bands  for  use  in  k-distribution  model  and 
the  verification  of  the  rapid  model  output. 
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Plates  9  &  1 0 


Conclusions 


Water  Vapor: 

•  The  dominant  greenhouse  gas  on  a  global  basis  for  the  troposphere 

-  For  atmospheres  with  significant  amounts  of  water  vapor, 
increases  in  concentration  of  a  uniformly  mixed  gas  result  in  a 
reduction  in  the  atmospheric  cooling  in  the  lower  troposphere 

•  Cooling  in  the  upper  troposphere  is  dominated  by  water  vapor  in 

the  pure  rotation  region  (0-400  cm'^), 

Carbon  Dioxide  and  Ozone: 

•  Strong  absorbers,  but  only  over  narrow  spectral  regions 

-  COg  eliminates  cooling  below  SOOmb  due  to  water  vapor 

-  from  640-690  cm’^ 

•  Doubling  of  carbon  dioxide  causes  decreased  cooling  in  the  lower 

troposphere  of  about  5% 

•  Ozone  exchange  heating  occurs  from  500  to  20  mb;  peaks  at 

0.3  K/day  at  50  mb  level 

Halocarbons,  CH4,  and  NgO 

•  These  gases  diminish  HgO  cooling  in  lower  troposphere  in  moist 

atmospheres,  but  have  little  affect  on  cooling  in  dry  atmospheres 

•  Effects  on  outgoing  flux  and  surface  downwelling  flux  are  non-negligible 

1 990-2000  Changes 

•  Radiative  effects  of  ozone  are  highly  dependant  on  depletion  profile 

•  Stratospheric  heating  effect  of  ozone  depletion  offsets  about  16% 

of  the  total  cooling  effect  of  carbon  dioxide  at  this  level 

•  700  cm-1  ozone  band  is  significant;  contributes  a  maximum  of  25%  of 

the  total  ozone  cooling  response  to  ozone  profile  perturbations 
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PRESENTATION  OVERVIEW 
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FOR  H2O,  CO2,  AND  0 


SUB-LAYER  INTEGRATION 
APPROACH* 
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.  Bernstein,  Applied  Optics,  Accepted,  May  1994. 


SUB-LAYER  INTEGRATION 
ERROR  ESTIMATION 
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LAYER  OPTICAL  DEPTH  (tJ 


LINE-OVERLAP  CORRECTION 
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VALUE  OF  <5ft/  =  0.01  cm  '  WORKS  WELL  FOR  FLUX 
AND  COOLING  RATE  CALCULATIONS. 
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FLUX  AND  COOLING  RATE 
COMPUTATIONAL  OVERVIEW 
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ICRCCM  DATA  BASE 
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-  SCHWARZKOPF  AND  PELS,  JGR,  96,  9075,  1991. 

-  RES  =  0.005  cm'^ 


Table  1.  Summary  of  Line-by-Line  and  MODTRAN  Flux  Calculations. 


Atmospheric 

Flux 

Model 

Surface  (W/m^ 

Tropopause^*^  (W/m^ 

Top  (W/m^ 

Case 

UP 

DOWN 

UP 

DOWN 

UP 

Tropical 

H2O  with  Continuum 

MODTRAN 

459.12 

389.21 

331.43 

3.22 

331.13 

AER 

458.92 

385.88 

332.39 

3.19 

332.59 

Mid-lattitude  Summer 

H2O  with  Continuum 

MODTRAN 

423.51 

335.84 

321.09 

7.56 

320.80 

AER 

423.53 

333.92 

321.17 

7.32 

320.96 

Subartic  Winter 

H2O  with  Contiauum 

MODTRAN 

247.72 

138.05 

222.49 

16.59 

221.35 

AER 

247.61 

138.23 

222.55 

16.78 

221.52 

MODTRAN 

423.51 

269.58 

334.98 

7.15 

334.77 

Mid-lattitude  Summer 

MODTRAN*^' 

423.51 

269.35 

335.34 

6.98 

335.12 

H2O  without  Continuum 

AER 

423.32 

269.02 

335.77 

6.89 

335.66 

GLA 

423.48 

267.28 

336.59 

6.25 

336.50 

GFDL 

423.56 

266.77 

336.61 

6.71 

336.32 

MODTRAN 

423.49 

76.38 

382.60 

12.48 

383.09 

Mid-lattitude  Summer 

GLA 

423.48 

76.48 

383.08 

12.14 

383.65 

CO2  (300  ppm) 

GISS^'* 

422.85 

76.35 

382.01 

12.22 

382.68 

GFDL 

423.56 

76.28 

383.26 

12.15 

383.72 

MODTRAN 

423.49 

84.40 

377.67 

14.23 

378.99 

MODTRAN<"> 

423.49 

76.22 

382.75 

12.38 

383.22 

Mid-lattitude  Summer 

CO2  (600  ppm) 

GLA 

423.48 

84.66 

377.95 

13.91 

379.37 

GLA® 

421.29 

83.70 

376.67 

13.89 

378.03 

GFDL® 

421.38 

83.60 

376.83 

13.92 

378.07 

Subarctic  Winter 

CO2  (300  ppm) 

MODTRAN 

247.71 

44.17 

231.50 

14.52 

229.83 

GLA 

247.71 

44.22 

231.64 

14.15 

230.17 

GFDL® 

247.36 

43.90 

231.44 

14.16 

229.94 

MODTRAN 

433.49 

6.33 

419.34 

3.26 

411.64 

MODTRAN® 

423.49 

6.16 

419.40 

3.06 

412.11 

Mid-lattitude  Summer 

0, 

GLA 

423.48 

6.00 

419.60 

3.10 

412.02 

GISS® 

422.85 

6.40 

418.61 

3.31 

410.84 

GFDL 

423.56 

6.18 

419.43 

3.04 

412.06 

(a)  99.7  mbar  (TROP),  179  mbar  (MLS),  282.9  mbar  (SAW). 

(b)  3a)  =  0.01  cm-'  (So)  =  0.0  cm’'  for  all  other  MODTRAN  runs). 

(c)  0  -  2,500  cm“^  spectral  interval. 

(d)  0  -  2,200  cm'*  spectral  interval. 
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Greenhouse  Gas  Concentration  Profiles  Retrieved  from  CIRRIS-IA 
Measurements  in  the  7-13  ^m  window  during  STS-39 

G  E  Bingham,  D  K  Zhou,  and  A  J  Steed 

Space  Dynamics  Laboratory,  Utah  State  University,  Logan,  UT  84341 
G  P  Anderson  and  R  M  Nadile 

Geophysics  Directorate,Phillips  Laboratory,  Hanscom  Air  Force  Base,  MA  01731 


Abstract 

The  spectral  and  spatial  distributions  of  the  greenhouse  gas  infrared  emissions  are 
measured  by  the  CIRRIS-IA  (Cryogenic  Infrared  Radiance  Instrumentation  for  Shuttle) 
interferometer  sensor  during  STS-39  flight  (28  Apr  -  6  May,  1991).  The  spectra  in  the 
important  7-11  fim  infrared  window  were  recorded.  High  resolution  spectral  scans  (~  1.5  cm'* 
instrument  resolution)  and  FASCOD-3  with  HITRAN92  database  and  NMC  temperature  and 
pressure  profiles  are  used  in  an  onion  peel  routine  to  retrieve  gas  concentrations  from  absolutely 
calibrated  radiance  data.  Scan  to  scan  and  profile  to  profile  variability  are  compared  to  show 
the  measurement  capability  of  the  cryogenic  infrared  interferometer.  The  local  aerosol  is 
estimated  by  comparing  the  modelled  and  measured  CO2  hot  band  (10.5  /xm)  spectra.  The 
results  of  concentration  profiles  retrieved  from  the  spectra  (30°  south  latitude  9°  East  longitude, 
local  time  of  1.2  hr.)  are  shown  in  the  figure.  Profiles  of  CFC-11,  CFC-12,  HNO5,  and  O3 
developed  in  this  fashion  and  at  different  locations  are  compared  with  ATMOS  (Atmospheric 
Trace  Molecule  Spectroscopy),  LIMS  (Limb  Infrared  Monitor  of  the  Stratosphere),  and  situ 
mean  data  for  these  gases. 


Conclusion 

The  infrared  emissions  from  greenhouse  gases  were  successfully  measured  as  part  of  the 
STS-39  mission  and  the  calibrated  spectra  used  to  obtain  the  spatial  distribution  of  some  of  these 
gases.  As  an  example,  the  vertical  concentration  profiles  of  O3,  HNO3,  CFC-12,  CFC-11,  and 
aerosols  were  retrieved  from  one  of  the  CIRRIS-IA  measured  spectral  profiles  collected  in  the 
southern  hemisphere  (30°  S,  9°  E).  Specific  errors  in  the  retrieval  routine  are  shown.  Other 
data  profiles  are  in  preparation  to  be  published  separately. 

The  CIRRIS-IA  data  set  demonstrates  the  capability  of  a  moderate  resolution 
interferometer  to  monitor  global  distributions  of  the  most  prominent  greenhouse  gases  in  the  7-13 
/xm  window.  Since  these  gases  are  also  active  in  the  ozone  reduction  cycle,  these  data  are 
doubly  important.  Use  of  the  COj  laser  band  at  940  cm'*  provides  an  effective  way  to  estimate 
aerosol  emissions  in  the  region  where  these  emissions  segregate  offset  trace  gas  retrieval.  A 
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major  limitation  of  the  CIRRIS-IA  instrument  complement,  as  it  was  flown,  was  the  specific 
instrumentation  with  the  lack  of  capability  to  measure  gas  ternperamre  and  pressure  while 
simultaneously  collecting  spectral  signatures.  One  of  the  large  error  sources  in  the  retrieval  may 
be  associated  with  the  temperature  and  pressure  fields  of  the  NMC  model.  The  accuracy  of  the 
tangent  height  determinations  associated  with  this  type  of  measurement  was  found  to  be 
extremely  critical  in  the  accuracy  of  the  analysis,  but  was  controlled  in  the  CIRRIS-IA 
measurements. 

Scientifically,  the  CIRRIS-IA  data  set  fills  an  important  gap  in  the  measurements  between 
the  ATMOS  and  UARS  data  sets.  The  CIRRIS-IA  data  set  does  not  provide  a  comprehensive 
global  sampling;  it  is  significant  because  it  was  collected  only  a  few  weeks  before  the  Pinatubo 
eruption.  The  density  profiles  retrieved  from  other  trace  gases  and  other  locations  measured  by 
CIRRIS-IA  have  also  been  performed  and  a  detailed  comparison  between  CIRRIS-IA  profiles 
and  comparisons  with  other  measurements  are  in  preparation. 


Figure  1.  An  example  of  greenhouse  gas  emission  spectrum  collected  at  22.4  hours  local  time 
on  April  29,  1991,  and  tangent  height  of  14.08  km,  latimde  of  42°  S  longitude  of  32°  W,  by 
interferometer  detector  5  filter  3  (770-2100  cm'^.  The  primary  emitters  are  shown  in  this 
region,  with  indication  of  analysis  and  retrieval  regions  used. 

Figure  2.  The  spectral  distribution,  collected  in  one  stair-step  mode  in  the  wavenumber  region 
of  800-1100  cm*  and  spatial  region  of  10-55  km.  The  data  were  taken  at  latitude  30°  S  and 
longitude  9°  E,  at  local  time  1.2  hours  on  April  29,  1991.  This  data  set  was  used  to  present 
the  VMR  profiles  in  Figure  9. 

Figure  3.  CFC-1 1 ,  CFC-12,  and  HNOj  spectral  overlapping  are  shown  in  the  region  of  830-960 
cm  ';  the  HNO,  spectral  emission  is  taken  off  of  CFC-1 1  and  CFC-12  spectrum  by  using  spectra 
of  different  tangent  heights.  The  major  radiance  peaks  are  identified.  The  data  were  collected 
at  22.4  hours  local  time  on  April  29,  1991,  latitude  of  42°  S  longitude  of  32°  W. 

Figure  4.  CFC-12  Ug  Q  branch  at  1161  cm  '  is  obtained  by  subtracting  two  different  tangent 
height  scans;  NjO  spectrum  is  also  shown  in  the  same  manner.  The  data  were  collected  at  22.4 
hours  local  time  on  April  29,  1991,  latitude  of  42°  S  longitude  of  32°  W. 

Figure  5.  An  example  of  comparison  of  CIRRIS-IA  and  modelled  spectra  in  the  CO2  laser  band 
region  at  a  tangent  height  of  17.93  km  which  is  used  to  retrieve  the  aerosol  extinction. 

Figure  6.  The  temperature  and  pressure  data  (from  NMC)  plotted  in  (a)  and  (b)  at  same  time 
and  location  of  the  CIRRIS-IA  data  collection  (30°  S,  9°  E).  The  temperature  and  pressure 
profiles  show  the  vertical  error  bar  of  altitude,  and  horizontal  error  bar  of  temperature.  The 
data  points  were  the  chosen  boundaries  of  the  onion-peel  layers,  (c)  the  retrieved  aerosol 
extinction  at  0.55  /xm  is  plotted  as  open  cycles  which  are  actual  measured  points.  Some  other 
measurements  are  reproduced  for  comparison. 
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Figure  7.  Examples  of  comparison  of  CIRRIS-IA  and  modelled  spectra  in  the  (a)  O3,  (b) 
HNO3,  (c)  CFC-12,  and  (d)  CFC-11  band  regions  which  are  used  to  retrieve  their  VMR  profiles. 

Figure  8.  The  comparison  between  measured  and  modelled  spectra  in  the  large  spectral  region 
of  800-1 100  cm  *  by  using  retrieved  VMR  profiles  of  O3,  HNO3,  CFC-12,  CFC-1 1 ,  and  aerosol. 

Figure  9.  The  VMR  profiles  of  (a)  O3,  (b)  HNO3,  (c)  CFC-12,  and  (d)  CFC-11  retrieved  from 
30°  S  and  9°  E,  Midnight  April  29,  1991  data  collection  (plotted  in  Figures  2  and  8).  The  open 
triangle  data  points  were  the  chosen  boundaries  of  the  onion-peel  layers,  and  open  cycles  which 
are  actual  measured  points  with  VMR  error  bar.  Some  other  measurements  are  reproduced  for 
comparison. 
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Abstract 

Methane  has  indirect  effects  on  climate  due  to  chemical  interactions  as  well  as 
direct  radiative  forcing  effects  as  a  greenhouse  gas.  We  have  calculated  the  indirect,  time- 
varying  tropospheric  radiative  forcing  and  GWP  of  O3  and  stratospheric  H2O  due  to  an 
impulse  of  CH4.  This  impulse,  applied  to  the  lowest  layer  of  the  atmosphere,  is  the 
increase  of  the  atmospheric  mass  of  CH4  resulting  from  a  25  percent  steady  state  increase 
in  the  current  emissions  as  a  function  of  latitude.  The  direct  CH4  radiative  forcing  and 
GWP  are  also  calculated.  The  LLNL  2-D  radiative-chemistry-transport  model  is  used  to 
evaluate  the  resulting  changes  in  the  O3,  H2O  and  CH4  atmospheric  profiles  as  a  function 
of  time.  A  correlated  k-distribution  radiative  transfer  model  is  used  to  calculate  the 
radiative  forcing  at  the  tropopause  of  the  globally-averaged  atmosphere  profiles.  The  O3 
indirect  GWPs  vary  from  ~27  after  a  20  yr  integration  to  ~4  after  500  years,  agreeing 
with  the  previous  estimates  to  within  about  10  percent.  The  H2O  indirect  GWPs  vary 
from  ~2  after  a  20  yr  integration  to  -0.3  after  500  years,  and  are  in  close  agreement  with 
other  estimates.  The  CH4  GWPs  vary  from  -53  at  20  yrs  to  -7  at  500  yrs.  The  20  year 
CH4  GWP  is  -20%  larger  than  previous  estimates  of  the  direct  CH4  GWP  due  to  a  CH4 
response  time  (-17  yrs)  that  is  much  longer  than  the  overall  lifetime  (10  yrs).  The 
increased  CH4  response  time  results  from  changes  in  the  OH  abundances  caused  by  the 
CH4  impulse.  The  CH4  radiative  forcing  results  are  consistent  with  IPCC  values. 
Estimates  are  made  of  latitude  effects  in  the  radiative  forcing  calculations,  and  UV  effects 
on  the  03  radiative  forcing  calculations  (10%). 


776 


I.  Introduction 


The  current  concerns  about  global  warming  and  ozone  depletion  have  necessitated 
detailed  studies  of  the  potential  greenhouse  effects  of  various  trace  gases  that  have  been 
introduced  into  the  atmosphere.  The  mechanism  usually  used  for  comparison  of  the 
greenhouse  potential  of  trace  gases  is  the  global  warming  potential  or  GWP  as  defined  by 
EPCC  (1990).  The  GWP  is  the  ratio  of  the  radiative  flux  change  at  the  tropopause  caused 
by  the  introduction  of  a  unit  mass  impulse  of  a  trace  gas  into  the  atmosphere  to  the 
radiative  flux  change  at  the  tropopause  caused  by  the  introduction  of  a  unit  mass  impulse 
of  CC)2-  An  essential  part  of  the  GWP  determination  is  the  calculation  of  the  radiative 
forcing,  which  is  defined  as  the  radiative  flux  change  at  the  tropopause  produced  by  a  unit 
change  in  the  number  of  molecules  of  a  particular  gas  with  all  other  abundances  held 
constant.  Usually  the  amount  of  abundance  change  used  in  a  radiative  forcing  calculation 
is  chosen  such  that  the  change  is  just  enough  to  produce  a  numerically  significant  flux 
change  value  at  the  tropopause.  Parameterized  experiments  for  radiative  forcing  have 
been  published,  for  example,  in  IPCC  (1990,  1992),  and  by  Ramanathan  et  al.  (1987).  A 
model  for  determining  GWPs  has  been  developed  by  Wuebbles  et  al.  (1992).  This  model 
approximates  GWP’s  with  an  accuracy  of  a  few  percent  and  has  been  published  by 
Wuebbles  and  Grossman  (1992)  for  CF4,  and  Grossman  and  Wuebbles  (1993)  for  SFe.  In 
addition  to  the  direct  GWP  and  radiative  forcing  discussed  above,  there  is  the  indirect 
GWP  and  radiative  forcing.  The  indirect  effects  are  caused  when  the  atmospheric 
chemical  processes  produce  a  change  in  the  abundance  of  species  2  after  the  introduction 
of  an  impulse  of  species  1.  A  GWP  and  radiative  forcing  effect  results  from  the  indirect 
abundance  changes  of  species  2  as  well  as  the  direct  GWP  and  forcing  effect  of  the 
change  in  abundance  of  species  1.  Indirect  GWP  effects  of  ozone  and  stratospheric  water 
vapor  have  been  discussed  in  IPCC  (1990,1992).  IPCC  (1990)  gives  an  indirect  ozone 
GWP,  due  to  methane  introduction,  of  24  at  a  time  of  20  years.  IPCC  (1992)  gives  an 
indirect  water  GWP  of  between  5  and  30  percent  of  the  direct  methane  GWP.  IPCC 
(1992)  recommended  that  the  indirect  GWP’s  should  be  withdrawn  due  to  incomplete 
understanding  of  the  atmospheric  chemical  processes  which  cause  the  ozone  changes. 
Yet,  according  to  EPCC  (1990),  the  indirect  ozone  forcing  due  to  methane  introduction 
could  represent  one  of  the  main  contributions  to  indirect  GWP  effects. 

The  main  purpose  of  this  paper  is  to  reexamine  the  indirect  effect  on  O3  and 
stratospheric  H2O  due  to  a  tropospheric  impulse  of  CH4  and  to  determine  the  resulting 
GWPs.  The  LLNL  2-D  radiative-chemistry-transport  model  is  used  to  generate  a  set  of 
annually  and  globally  averaged  model  atmospheres  which  contain  the  perturbed  €>3,  H2O, 
and  CH4  abundances  as  a  function  of  time.  A  correlated  k-distribution  model  for  the 
atmospheric  transmission  is  used  to  calculate  the  tropospheric  radiative  forcing  for  CH4, 
O3,  H2O,  and  CO2.  These  forcing  values  are  then  used  to  calculate  the  direct  GWP  for 
CH4,  and  the  indirect  GWP  for  O3  and  H2O.  Estimates  are  made  for  the  UV  contributions 
to  the  indirect  O3  radiative  forcing  values  and  latitude  effects  on  the  radiative  forcing 
calculations  for  O3,  H2O,  and  CH4.  Lastly,  we  make  suggestions  for  future,  more 
advanced,  indirect  O3  and  H2O,  GWP  calculations. 
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n.  Model  Atmospheres 

A.  Model  Structure 

The  version  of  the  LLNL  zonally-averaged  two-dimensional  chemical-radiative- 
transport  model  (Kinnison  et  al.,  1994,  Patten  et  al.,  1994,  Wuebbles  et  al.,  1991)  used  for 
these  calculations  determines  the  atmospheric  distributions  of  44  chemically  active 
atmospheric  trace  constituents  in  the  troposphere  and  stratosphere.  The  model  domain 
extends  from  pole  to  pole,  and  from  the  ground  to  60  km.  The  sine  of  latitude  is  used  as 
the  horizontal  coordinate,  with  uneven  increments  corresponding  to  approximately  10 
degrees  in  latitude.  The  vertical  coordinate  corresponds  to  the  natural  logarithm  of 
pressure;  the  scale  height  is  7.2  km  and  the  surface  pressure  is  1013  mb.  The  vertical 
resolution  is  1.5  km  in  the  troposphere  and  3  km  in  the  stratosphere. 

B.  Model  Transport 

In  the  dynamics  portion  of  the  two-dimensional  code,  the  circulation  field  is 
obtained  following  the  approach  of  Garcia  and  Solomon  (1983).  Since  this  approach  is 
only  partially  implemented  in  the  model  at  the  present  time,  we  use  climatological 
temperature  and  mean  zonal  wind  fields.  Both  the  temperature  and  mean  zonal  wind  field 
vary  continuously  over  the  annual  cycle.  The  temperature  field  is  based  on  the  reference 
model  of  Barnett  and  Comey  (1985),  while  the  mean  zonal  wind  is  based  on  data  from 
Homing  et  al.  (1988).  The  method  used  to  obtain  the  circulation  is  based  on  combining 
the  zonal  mean  momentum  equation  and  the  thermodynamic  equation  into  a  form  that, 
along  with  the  thermal  wind  equation,  yields  a  second  order  diagnostic  equation  for  the 
residual  mean  meridional  stream  function.  The  net  heating  rate  is  calculated  knowing  the 
temperature  and  chemical  species  distribution  and  includes  latent  heating.  The  present 
code  does  not  include  a  separate  planetary  or  gravity  wave  calculation,  but  instead 
calculates  the  Eliassen-Palm  flux  directly  from  the  zonal  mean  momentum  equation.  This 
is  done  by  solving  the  zonal  mean  momentum  equation. 

The  transport  of  chemical  species  is  accomplished  through  both  advection  and 
turbulent  eddy  transport.  Advection  terms  are  treated  using  the  second  order,  two 
dimensional  transport  algorithm  of  Smolarkiewicz  (1984).  The  transport  caused  by  eddy 
motion  is  parameterized  through  the  diffusion  coefficients,  Kyy  and  Kzz.  Values  of  Kzz 
do  not  vary  with  season,  but  are  different  between  the  troposphere  and  stratosphere.  The 
value  of  Kzz  in  the  troposphere  is  4  m^  sec-^  while  Kzz  in  the  stratosphere  has  values  of 
0.1  m2  sec”l  in  the  lower  stratosphere  increasing  to  0.25  m^  sec“l  near  the  middle 
stratosphere  and  returning  to  0.1  m2  sec“l  at  the  strato pause.  Values  of  Kyy  in  the 
stratosphere  are  calculated  using  a  similar  method  to  that  of  Newman  (1988)  by  dividing 
the  zonal  mean  momentum  equation  (i.e.,  the  Eliassen-Palm  flux)  by  the  horizontal 
gradient  of  the  quasi-gepstrophic  potential  vorticity.  In  the  stratosphere,  a  minimum  value 
of  Kyy  has  been  established  as  1  x  10^  m2  sec”^.  The  values  of  Kyy  in  the  troposphere  are 
assigned  a  value  of  2  x  10^  m2  sec-^  The  tropopause  is  determined  from  the  latitude 
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dependent  annual  mean  temperature  field.  The  height  of  the  tropopause  varies  from  8  km 
near  the  poles  to  14  km  at  the  equator. 

Two  physical  removal  processes  are  included  in  the  model.  Surface  deposition 
occurs  in  the  first  layer  of  the  model.  Most  species  are  assigned  a  st^d^d 
velocity  of  1  cm/sec.  Nitric  oxide,  nitrogen  dioxide,  ozone,  molecular  ^ 

formaldehyde  are  assigned  unique  deposition  velocities.  Also, 
surface  (e.g.  methane)  are  assigned  a  deposition  '^®l®®'^y 
parameterized  in  the  model  as  a  first-order  loss  process  that  occurs  m  the  first  6-7  altitude 

layers. 

C.  Solar  Radiative  Transfer  Model 

To  capture  the  spectral  detaU  needed  for  photodissociation  calculations,  the  two- 
stream  multiple-layer  UV-visible  model  uses  126  wavelength  bins  between  175  nm  Md 
^sHm  The  twlstream  approach  was  chosen  because  of  the  requirements  for 
computational  efficiency  placed  on  radiative  transfer  models  designed  for  ^ 

atmospheric  chemistry  models.  In  this  approach,  solar  radiation  is 

direct  solar  radiation,  downward  diffuse  radiation,  and  upward  diffuse  radiation.  The 
scattering  of  energy  from  the  direct  solar  beam  within  each  individual  layer  is  treated 
using  the  delta-Eddhigton  algorithm  (Joseph  et  al.  197^,  which 

S  scattering  and  absorption  on  the  solar  zenith  angle.  The  scattenng  of  diffu^  radiation 
(i.e.,  previously  scattered  radiation)  from  each  individual  layer  is 
dmnlL  Saean-Pollack  (1967)  algorithm.  Both  algorithms  allow  inclusion  of  the  bulk 
Z^JpK^of  cUils  mid  irosols.  Finally,  dm  adding  n.ed.od  is  used  Ui  calculam 
irradiances  throughout  the  vertically  mhomogeneous  atmosphere. 

D.  Infrared  Radiative  Transfer  Model 

The  infrared  model  is  a  version  of  the  model  described  by  Harshvardhan  et  al. 
fl987).  For  our  usage,  it  has  been  modified  to  improve  the  accur^y  in  t  e  upper 
stratosphere.  It  includes  absorption  and  emission  by  CO2,  O3,  and  H2O.  It  is  based  on 
wide-band  parameterizations  fit  to  line-by-line  calculations.  Inhomogeneous  a  sorp  on 
paths  are  included  by  pressure  and  temperature  weighted  scaling  of  trace  gas  absorber 
Liounts  The  model  provides  for  specification  of  fractional  cloud  cover  within  each 
vertical  model  layer.  Separate  fractions  can  be  specified  for  convective  (deep, 
overlapping)  and  randomly  overlapped  clouds. 

E.  Model  Chemistry 

The  mechanism  we  are  currently  using  in  the  LLNL  two-dimensional  model  has  the 
following  characteristics: 

.  A  total  of  44  species  are  allowed  to  vary  with  time  in  the  troposphere  and 
stratosphere.  Four  additional  species  (O,  0(>D),  N,  H)  ate  treated  as  steady- 
State. 
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There  are  103  thermal  (two-  and  three-body)  reactions  and  48  photolytic 
reactions  in  the  mechanism. 


•  A  complete  representation  of  known  stratospheric  chemical  processes. 

The  mechanism  is  listed  in  its  entirety  in  Table  1.  The  reactions  are  classified  by  type: 
two-body,  three-body,  and  photolytic.  The  kinetic  parameters  for  two-body  reactions  are 
included  as  well. 

The  basic  photochemistry  represents  the  tropospheric  and  stratospheric  interactions 
of  actinic  solar  flux  with  the  individual  species  in  the  families  Ox,  NOy,  ClOy,  HOy, 
BrOy,  as  well  as  with  CH4  and  its  oxidation  products.  This  part  of  the  mechanism 
incorporates  44  transported  species  and  4  species  for  which  abundance  is  determined 
through  the  assumption  of  instantaneous  equilibrium.  The  thermal  reactions  considered 
number  103,  the  photolytic  reactions  number  48.  Source  gases  used  include  NOx,  N2O, 
CH4,  CO2,  CO,  the  chlorine  compounds  CFC-11,  CFC-12,  CFC-113,  CFC-114,  CFC- 
115,  HCFC-22,  CCI4,  CH3CCI3,  CH3CI,  and  the  bromine  compounds  CH3Br,  CF2ClBr, 
and  CF3Br.  Most  of  the  thermal  reaction  rates  were  taken  from  the  NASA  Panel 
recommendations  provided  in  DeMore  et  al.,  (1992).  The  rate  constant  for  CH4  +  OH 
calculated  by  Vaghjiani  and  Ravishankara  (1991)  was  used  because  it  represents  the  most 
recent  work.  Absorption  cross  section  information  has  been  assembled  primarily  from 
DeMore  et  al.  (1992)  and  Gillotay  and  Simon  (1991a,b).  Water  vapor  is  dealt  with  in  a 
imique  way  compared  to  the  other  chemical  species.  The  level  of  water  vapor  is  assigned 
its  climatological  value,  following  the  specific  humidity  of  the  Oort  climatology, 
everywhere  below  the  hydropause.  Above  the  hydropause,  water  vapor  treated  as  a 
transported  species;  methane  oxidation  is  its  major  source  in  this  region  of  the 
atmosphere. 

The  photolytic  loss  rate  constants  are  calculated  by  integrating  the  product  of 
absorption  coefficient,  quantum  yield,  and  solar  flux  over  wavelength  (i.e.  from  175  nm 
to  760  nm).  The  exoatmospheric  solar  flux  is  from  WMO,  1985.  The  solar  flux  is  then 
calculated  as  a  function  of  altitude,  latitude,  and  season  at  each  time  step,  including  the 
effects  of  absorption  by  O2  and  O3  and  multiple  molecular  (Rayleigh)  scattering.  The 
absorption  cross  sections  and  quantum  yields  include  temperature  and  pressure 
dependence  where  appropriate  and  available.  The  Schumann-Runge  band  region  of  the 
02  absorption  is  treated  using  the  technique  of  Allen  and  Frederick  (1982),  modified  to 
match  the  more  recent  lower  cross  sections  of  the  O2  Herzberg  continuum  region.  The 
photolysis  of  NO  is  also  treated  through  the  parameterization  technique  of  Allen  and 
Frederick  (1982). 

The  nonlinearity  of  the  photochemistry  with  respect  to  diurnal  averaging  is 
accounted  for  through  the  calculation  of  individual  altitude,  latitude,  and  seasonally 
varying  factors  for  each  photochemical  process.  Full  diurnal  variability  of  each  species  is 
calculated  off-line  for  four  seasons.  The  factors,  relating  individual  process  rates  from  the 
full  diurnal  calculation  to  the  diumally-averaged  values,  are  then  spline-interpolated  for 
each  time  step. 
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Members  of  the  ClOy  family  currently  thought  to  play  a  role  in  the  Antarctic  spring, 
such  as  the  CIO  dimer,  are  included  in  this  version  of  the  model.  Polar  stratospheric 
aerosol  surface  has  not  been  included.  A  representation  of  background  stratospheric 
sulfuric  acid  aerosol  surface  is  incorporated  into  the  model,  following  the 
recommendations  from  the  WMO,  1989  report  The  reactions  of  N2O5  and  C10NC)2  with 
H2O  on  the  aerosol  surface  are  parameterized  as  first-order  loss  processes  with  rate 
constants  determined  by  specified  surface  area  density,  collision  frequency,  and  reaction 
probability. 

F.  Emission  and  Boundary  Conditions 

In  our  modeling  of  current  and  possible  future  atmospheres,  we  need  to  represent 
the  surface  (and  atmospheric  emissions)  and/or  surface  concentrations  of  a  number  of 
constituents  that  affect  tropospheric  and  stratospheric  chemistry.  A  great  deal  of  evidence 
indicates  that  surface  emissions  and  concentrations  of  globally  important  trace  gases  are 
increasing  (see  WMO,  1991  for  an  overview).  During  the  last  ten  years,  scientists  have 
been  determining  the  atmospheric  concentrations  and  emission  rates  of  a  variety  of  trace 
species  (e.g.  nitrogen  oxides,  carbon  monoxide,  and  hydrocarbons)  on  a  regional  as  well 
as  global  scale.  These  trace  gases  can  affect  the  tropospheric  oxidizing  capacity  due  to 
changes  in  the  levels  of  ozone  and  hydroxyl.  To  quantify  these  effects,  researchers  use 
global  atmospheric  models  to  investigate  the  effects  of  trace  species.  Emission 
inventories  and  measurements  of  the  source  gases  are  required  for  the  initialization  as 
well  as  validation  of  these  models.  The  source  gas  fluxes  must  be  estimated  correctly  in 
order  for  the  models  to  calculate  accurate  species  concentrations.  Discrepancies  between 
measured  and  calculated  species  concentrations  can  be  used  to  check  the  consistency  of 
emissions  inventories,  assuming  all  other  model  parameters  are  correct.  The  emphasis  up 
to  now  has  been  to  use  global  models  to  study  these  effects.  One  goal  of  these  studies  is 
to  investigate  the  impact  on  the  troposphere — ozone,  hydroxyl,  and  oxidizing  capacity. 

The  surface  boundary  conditions  for  the  ambient  version  of  the  two-dimensional 
model  consist  of  fixed  concentration  for  nitrous  oxide,  molecular  hydrogen,  hydrochloric 
acid,  and  the  halogenated  organics  and  fluxes  for  methane,  nitrogen  dioxide,  and  carbon 
monoxide.  The  emission  rates  are  specified  for  each  latitude  zone  in  Table  2. 

Methane  emissions  are  distributed  latitudinally  according  to  the  recommendations 
of  Fung  et  al.,  (1991).  Carbon  monoxide  emissions  are  distributed  latitudinally  according 
to  the  emission  categories  defined  by  Hough  (1991).  The  methane  and  carbon  monoxide 
emission  rates  used  for  these  calculations  are  obtained  by  determining  the  steady  state 
chemical  loss  rate  from  a  previous  ambient  model  run.  However,  seasonal  variations  in 
the  source  gas  emissions  are  not  included  in  the  results  presented  here.  This  can  lead  to 
discrepancies  between  the  magnitudes  of  the  calculated  and  observed  species 
distributions.  Another  potential  source  of  discrepancy  is  the  representation  of  advective 
and  diffusive  transport  in  the  troposphere.  This  would  adversely  affect  short-lived  species 
such  as  nitrogen  oxides.  In  the  current  model,  the  diffusion  coefficients  are  constant 
across  all  latitudes  and  seasons  in  the  troposphere.  We  currently  use  diffusive  transport  to 
represent  fast  vertical  transport  processes  such  as  convection.  If  the  diffusivities  are  too 
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low,  then  species  will  undergo  very  slow  vertical  mixing.  Fast  vertical  transport  is  often 
implicated  in  explaining  unusually  high  levels  of  short-lived  species  in  the  middle  and 
upper  troposphere. 

G.  Ambient  Atmosphere 

The  2-dimensional  model-derived  ambient  atmosphere  used  for  these  calculations  is 
based  on  1990  emissions  of  the  source  gases.  The  two-dimensional  (i.e.  latitude,  altitude) 
distributions  of  the  radiatively  important  trace  gases — water  vapor,  carbon  dioxide, 
ozone,  nitrous  oxide,  and  methane — are  averaged  over  the  annual  cycle  and  from  pole  to 
pole.  The  resulting  globally-  and  annually-averaged  altitude  profiles  for  water  vapor, 
ozone,  and  methane  are  depicted  in  Figure  1.  Carbon  dioxide  is  constant  at  350  ppmv 
throughout  the  model  atmosphere.  Nitrous  oxide  decreases  monotonically  from  308  ppbv 
near  the  surface  to  1.2  ppbv  at  59  km.  The  pressure-temperature  profile  for  the  model 
atmosphere  is  shown  in  Figure  2.  The  tropopause  in  the  globally-averaged  atmosphere,  is 
specified  as  the  altitude  at  which  the  temperature  gradient  in  the  troposphere  decreases  to 
2K/km.  As  indicated  in  Figure  2,  this  occurs  at  166  mb  (-13.2  km). 

H.  Time  Dependent  Atmosphere  Calculations 

The  two-dimensional  model  is  used  to  determine  the  changes  in  ozone  following  the 
pulse  increase  in  the  abundance  of  methane.  Due  to  its  long  lifetime  in  the  troposphere, 
the  assumed  latitude  and  season  of  the  methane  pulse  should  not  impact  significantly  the 
derived  changes  in  ozone.  Therefore,  the  pulse  is  distributed  latitudinally  at  the  start  of 
the  run  as  a  percentage  increase  of  the  ambient  methane  concentrations.  The  pulse  needs 
to  be  sufficiently  large  for  the  chemical  model  to  determine  an  ozone  change  providing  a 
significant  change  in  radiative  forcing,  yet  sufficiently  small  for  the  nonlinear  effects  on 
atmospheric  chemistry  to  be  minimal.  Based  on  the  analysis  of  ozone  radiative  forcing  in 
Lacis  et  al.  (1990),  we  estimate  that  a  sufficiently  large  change  in  tropospheric  average 
ozone  concentration  should  be  approximately  4%.  This  corresponds  to  a  25%  increase  of 
current  methane  emissions  in  our  model.  For  the  calculations  in  this  paper  we  estimate 
the  equilibrium  change  in  the  total  atmospheric  mass  of  methane  for  a  25%  emission 
increase  based  on  model  results  for  several  increased  flux  scenarios.  The  difference 
between  the  perturbed  and  ambient  burdens,  1702.9  teragrams,  is  added  to  the  surface 
layer  in  the  model.  We  then  run  the  model  for  thirty  years  using  ambient  CH4  emission 
fluxes  and  initial  concentrations  for  the  other  species. 

Figure  3  shows  the  integrated  change  in  the  mass  of  methane  calculated  in  the 
model  as  a  function  of  time  before  and  after  the  initial  pulse  as  well  as  the  ambient  value 
of  the  methane  mass.  The  methane  pulse  decays  approximately,  exponentially  in  time 
with  a  lifetime  of  the  order  of  17  years.  The  decay  function  with  time  indicates  that 
significant  nonlinear  feedbacks  with  hydroxyl  have  occurred.  The  methane  depletion  rate 
is  proportional  to  both  the  concentration  of  methane  as  well  as  the  concentration  of 
hydroxyl.  Since  the  hydroxyl  lifetime  is  of  the  order  of  five  minutes,  it  can  be  assumed 
that  the  hydroxyl  concentration  is  in  equilibrium.  The  equilibrium  hydroxyl  concentration 
is  an  inverse  function  of  the  methane  concentration. 
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(1) 


^CH^\ldt  =  -ki[CH4\l(k2-^ki[CH^'^+SQ  , 

where  So  is  the  background  methane  flux  input  into  the  atmosphere,  and  [CH4]  is  the 
concentration  of  methane.  The  value  of  [CH4]  is  the  sum  of  the  background  atmosphere 
methane  plus  the  extra  methane  due  to  the  pulse.  Eventually  the  methane  added  in  the 
pulse  will  decay  back  to  the  background  atmosphere  methane  abundance.  Equation  1 
implies  that  the  methane  pulse  decay  will  be  a  complex  function  of  the  instantaneous 
methane  abundance,  but  it  wiU  closely  approximate  an  exponential  decay  with  a  lifetime 
of  the  order  of  17  years.  While  the  instantaneous  methane  lifetime  in  the  ambient 
atmosphere  is  10-12  years,  the  conjecture  made  in  IPCC  (1990)  that  a  methane  pulse  will 
decay  exponentionally  with  a  10  year  lifetime  will  not  be  the  case  for  real  methane 
pulses.  Other  analyses  (Fisher,  1993,  Prather,  1993)  have  determined  that  an  equivalent 
lifetime  of  the  order  of  17  years  appropriately  describes  the  lifetime  of  a  methane  pulse  in 
the  atmosphere.  We  emphasize  the  point  here  that  we  do  not  use  an  exponential  decay 
function  in  this  paper,  we  use  a  direct  integration  of  the  actual  radiative  forcing  which 
contains  the  actual  value  of  the  methane  abundance  at  any  particular  time.  Since  the  IPCC 
(1990)  GWP  numbers  were  based  on  an  exponential  methane  decay  function  with  a  10 
year  lifetime,  we  expect  differences  in  the  GWP  values.  This  will  be  discussed  in  Section 
V.  below. 

Figure  4  depicts  the  annually-averaged  change,  with  time,  of  the  total  ozone  column 
calculated  after  the  methane  pulse.  Figure  5  shows  the  calculated  change  with  time  of  the 
annually-averaged  tropospheric  ozone  column.  Ambient  values  of  the  total  and 
tropospheric  values  of  the  ozone  column  are  shown  in  Figures  4  and  5.  Total  global 
ozone  takes  about  six  years  before  reaching  the  maximum  effect  while  the  troposphere 
responds  within  two  years.  The  change  in  the  ozone  column  abundance  is  a  result  of 
methane  oxidation.  In  tropospheric  regions  characterized  by  high  nitric  oxide 
concentrations,  the  net  methane  oxidation  reaction  is  (cf.  Wuebbles  and  Tamaresis, 
1993): 

CH^  + 12  -  7.5O2  ->  +  0-4  -  0.50//  -H  3.6  -  3.80,  +  \2-  IAH2O + 0.5  -  O.6H2  . 

We  see  that  complete  methane  oxidation  results  in  net  production  of  ozone.  The 
difference  in  response  times  between  the  troposphere  and  stratosphere  is  due  to  transport, 
particularly  across  the  tropopause.  Since  the  pulse  is  emitted  in  the  surface  zone, 
dynamical  processes  will  transport  the  methane  upwards.  Typically,  it  takes  a  long-lived 
gas  emitted  at  the  surface  about  four  years  to  mix  into  the  stratosphere.  Furthermore,  we 
observe  that  the  tropospheric  ozone  column  peaks  two  years  after  the  pulse.  Thus  if  the 
time  constant  for  chemistry  (two  years)  is  added  to  that  for  transport  (four  years),  we 
estimate  that  it  should  take  six  years  for  the  stratospheric  ozone  column  to  peak. 

Figure  6  shows  the  time  variation  of  the  increase  in  the  annually  and  globally 
averaged  stratospheric  water  vapor  burden  (kg).  The  peak  value  of  the  H2O  change 
occurs  at  a  time  of  10-1 1  years  after  the  methane  pulse.  According  to  Lelieveld  et  al. 
(1993)  the  H2O  levels  in  the  lower  stratosphere  are  largely  dependent  on  transport  of  CH4 
and  O3  from  the  troposphere.  Since  the  stratospheric  ozone  peaks  approximately  6  years 
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after  the  methane  we  expect  the  stratospheric  water  vapor  to  take  at  least  this  amount  of 
time  and  possibly  more  to  reach  its  peak  value. 

It  is  also  important  to  consider  the  uncertainties  in  these  calculations  that  may  result 
both  from  the  using  a  two-dimensional  model  in  the  analysis  of  tropospheric  (and 
stratospheric)  chemistry  effects  and  from  uncertainties  in  the  chemistry  itself.  We  are  not 
able  to  define  precisely  the  limits  on  the  uncertainties  at  this  time;  however,  we  provide  a 
range  as  a  representation  of  overall  uncertainties. 

One  major  uncertainty  is  in  the  amount  of  production  of  ozone,  in  the  troposphere, 
per  unit  oxidation  of  methane.  This  uncertainty  depends  on  the  reaction  of  CH3O2  +  NO 
to  produce  NO2  which  then  photolyzes  (in  the  troposphere)  and  produces  ozone.  The 
primary  uncertainty  in  the  rate  of  ozone  formation  will  then  largely  depend  on  the  amount 
of  NOx  in  the  troposphere  (we  argue  that  the  well-recognized  nonlinear  relationship 
between  NOx  and  ozone  is  unimportant  to  the  effect  of  methane  oxidation  to  produce 
ozone,  as  the  production  of  ozone  depends  only  on  the  NO  conversion  reaction  given 
above).  Unfortunately,  tropospheric  NOx  concentrations  are  poorly  characterized 
globally.  The  two-dimensional  model  gives  a  range  in  the  tropospheric  concentrations  of 
NOx  of  1.8-416  pptv.  By  comparison,  the  MLOPEX  study  (Mauna  Loa  Observatory 
Photochemistry  Experiment)  over  the  Pacific  Ocean  found  a  range  of  2.4-358  pptv  * 
(Carroll  et  al.,  1992).  Measurements  over  polluted  areas  can  be  much  larger,  giving  NOx 
levels  >1  ppbv.  Therefore,  the  model  appears  to  represent  the  expected  values  in  the  clean 
troposphere  reasonably  well,  but  may  greatly  underestimate  the  values  in  polluted  areas. 
The  impact  of  these  polluted  areas  on  the  derived  effects  of  methane  will  not  be  known 
until  we  are  able  to  better  characterize  the  global  amounts  of  NOx-  Most  of  the 
uncertainty  in  the  model  representation  of  NOx  should  be  in  underestimating  the  effects 
of  polluted  areas,  but  this  will  also  depend  on  how  well  the  model  represents  the 
variability  of  NOx  in  less  polluted  areas  of  the  troposphere.  At  this  time,  we  believe  that 
we  can  safely  assume  that  the  total  uncertainty  is  no  more  than  a  factor  of  two. 

in.  Radiative  Forcing  Models 

A.  Infrared  Radiative  Transfer  Model 

We  have  used  an  infrared  radiative  transfer  model  based  on  the  correlated  k- 
distribution  algorithm  (Grant  et.  al.,  1992,  Grossman  and  Grant,  1992)  to  calculate  the 
radiative  forcing  at  the  tropopause.  Calculations  were  done  for  a  globally  averaged 
atmosphere  for  CH4  and  O3  at  times  between  1  and  30  years  after  a  CH4  impulse  has 
been  applied  at  the  bottom  of  the  atmosphere.  Calculations  have  been  made  of  the 
radiative  forcing  for  CO2,  in  the  globally-averaged  atmosphere  and  at  five  different 
latitudes,  2  years  after  the  pulse.  This  was  done  in  order  to  determine  the  scaling  factor 
with  respect  to  the  IPCC  (1990)  expression  for  CO2  radiative  forcing  at  the  tropopause. 
Calculations  have  also  been  made  at  five  different  latitudes  for  CH4  and  O3,  at  a  time  of  2 
years,  in  order  to  determine  the  differences  between  the  radiative  forcing  in  the  globally- 
averaged  atmosphere  and  latitude  averaged  radiative  forcing.  The  tropospheric  radiative 
forcing  of  H2O,  caused  by  the  methane  pulse,  has  been  calculated  at  five  different 
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latitudes  for  times  of  5,  10,  15,  and  20  years  after  pulse.  An  appropriate  latitudinal 
average  is  calculated  to  determine  the  globally-averaged  tropospheric  radiative  forcing. 
The  parameters  of  the  calculations  are; 

O3,  CH4,  wave  number  range  of  900—1650  cm  25  cm  ^  subintervals, 

CO2,  wave  number  range  of  550-1200, 2000-2525  cm-l,  25  cm-1  subintervals, 

H2O,  wave  number  range  of  25—3000  cm“^,  25  cm~^  subintervals, 

Tropopause  of  globally  averaged  atmosphere  at  166  mb. 

Ground  Temperature  of  29 IK  (Peixoto  and  Oort,  1992). 

The  radiative  forcing  is  defined  as  the  difference  between  the  net  infrared  flux  at  the 
tropopause  in  the  ambient,  globally  averaged  atmosphere,  and  the  net  flux  at  the 
tropopause  in  the  atmosphere  with  changed  CH4,  H2O,  and  O3  profiles  calculated  at  the 
various  times  after  the  CH4  impulse. 

A  test  of  the  ability  of  correlated  k-distribution  model  to  calculate  O3  tropospheric 
forcing  was  performed  using  an  LLNL  generated  model  for  an  annually-averaged,  rnid 
latitude,  clear  sky  atmosphere.  A  perturbation  of  15%  of  the  tropospheric  O3  w^  added 
and  the  tropospheric  forcing  calculated.  The  result  of  0.266  W/m^  agreed  to  witto  about 
10  percent  of  the  Ramanathan  et  al.  (1987)  value  for  a  similar  test  case,  assutriing  a  20 
percent  correction  factor  for  clouds.  A  test  calculation  for  CH4  radiative  forcing  in  a  mid¬ 
latitude  summer,  clear  sky  approximation,  was  made  by  Grossman  et  al.  (1993)  and  the 
result  agreed  with  the  Ramanathan  et  al.  (1987)  value  to  within  about  4  percent.  A 
calculation  of  the  tropospheric  radiative  forcing  for  CH4  in  a  globally  and  annually 
averaged  atmosphere  both  with  and  without  a  representative  cloud  distribution  was  made 
by  Grossman  and  Grant  (1993).  The  results  with  clouds  agreed  to  within  1  percent  of  the 
DPCC  (1990)  value.  These  tests  indicate  that  our  correlated  k-distribution  model  has  the 
ability  to  provide  consistent  tropospheric  radiative  forcing  results.  Table  3  shows  the 
radiative  forcing  values  in  W/m^  at  representative  times  between  1  and  30  years  after  the 
CH4  impulse  for  both  the  direct  CH4  and  the  indirect  O3  radiative  forcing.  Also  shown 
are  the  results  of  direct  CO2  tropospheric  radiative  calculations  for  year  2  (maximum  O3 
radiative  forcing)  using  CO2  concentrations  of  1.29  and  2  times  normal  CO2  abundanc^ 
The  CO2  calculations  were  made  to  test  the  systematic  differences  between  the  IPCC 
(1990)  formula  for  CO2  tropospheric  radiative  forcing  and  the  results  obtained  using  the 
LLNL  model  atmosphere.  For  both  cases  shown  in  Table  3  the  results  of  our  calculation 
produce  a  radiative  forcing  which  is  approximately  25  percent  higher  than  that  calculated 
using  the  IPCC  (1990)  formula.  IPCC  (1990)  notes  that  calculations  which  do  not  use  a 
cloud  cover  model,  which  is  the  case  for  our  calculations,  will  result  in  radiative  forcing 
values  which  are  of  the  order  of  20  percent  too  high.  Thus  the  agreement  of  our  CO2 
radiative  forcing  values  with  the  IPCC  (1990)  values  would  be  of  the  order  of  5  percent  if 
we  had  included  a  cloud  cover  model  in  our  calculations.  It  is  important  for  the  GWP 
calculations,  which  will  be  discussed  later,  that  the  CO2  forcing  values  derived  from  the 
IPCC  (1990)  formula  contain  a  25  percent  adjustment  in  order  to  be  comparable  with  the 
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03  and  CH4  forcing  values  given  in  Table  3.  The  results  for  the  time  variation  of  the 
indirect  O3  radiative  forcing,  given  in  Table  3  show  that  the  peak  forcing  occurs  at  a  time 
of  2  years  after  the  initial  methane  impulse.  This  is  indicative  of  the  radiative  forcing 
being  correlated  with  the  time  variation  of  the  tropospheric  ozone  column  density  which 
peaks  at  a  time  of  two  years  (Fig.  5)  rather  than  the  total  ozone  column  density  which 
peaks  at  a  time  of  about  6  years  (Fig.  4),  The  reason  for  this  behavior  of  the  O3  radiative 
forcing  is  that  the  upward  IR  flux  change  is  much  larger  than  the  downward  flux 
difference  at  the  tropopause  at  all  times  and  this  drives  the  temporal  behavior  of  the  O3 
forcing  even  though  the  change  in  the  stratospheric  ozone  column  is  much  larger  than  the 
change  in  the  tropospheric  ozone  column. 

Two  interesting  aspects  of  the  use  of  a  2-D  latidudinally  averaged  model 
atmosphere  are  the  (Question  of  the  variation  of  the  tropospheric  radiative  forcing  as  a 
function  of  latitude  and  the  difference  between  the  latitudinally  averaged  value  of  the 
tropospheric  radiative  forcing  and  the  tropospheric  radiative  forcing  calculated  for  the 
globally  averaged  model  atmosphere.  We  have  calculated  the  direct  CH4  and  indirect  O3 
tropospheric  for  five  different  latitudes  between  +71  deg.  and  -71  deg  for  the  time  at 
which  maximum  O3  radiative  forcing  occurs,  year  2  after  the  CH4  impulse.  The  results 
are  shown  in  Table  4.  The  radiative  forcing  for  O3  and  CH4  has  been  shown,  for  a  CH4 
impulse  of  2.83  times  the  concentration  of  the  lowest  atmospheric  level  for  each  latitude 
(25  percent  of  the  total  tropospheric  CH4  mass),  at  the  actual  tropopause  at  that  latitude. 
There  are  significant  differences  between  the  radiative  forcing  values  calculated  in  the 
polar  regions  and  those  calculated  for  the  mid  latitude  and  the  equatorial  regions.  The 
polar  values  are  lower  and  the  mid  latitude  values  are  higher  than  the  average  value  of  the 
forcing.  The  maximum  latitudinal  variations  are  +12,  -9  percent  for  the  CH4  forcing  and 
+29,  -49  percent  for  the  O3  forcing.  These  excursions  from  the  average  may  have 
significant  climatalogical  consequences  which  will  be  discussed  later  from  the  point  of 
view  of  the  GWP  results.  It  is  interesting  to  note  that  the  results  shown  in  Table  4  indicate 
that  the  radiative  forcing  values  calculated  for  the  globally  averaged  atmosphere,  0.2455 
for  CH4  and  0.1180  for  O3,  differ  by  ~5  percent  (CH4)  and  ~9  percent  (O3)  with  the 
average  of  the  forcing  values  for  each  latitude.  We  consider  these  differences  to  be  small 
and  thus  the  use  of  the  globally  averaged  atmosphere  to  determine  the  average  radiative 
forcing  is  a  reasonable  approach  to  the  GWP  calculation. 

Table  5  shows  the  latitude  effect  of  the  tropospheric  CO2  radiative  forcing  for  a 
CO2  impulse  of  twice  the  normal  concentration  (350  ppm).  The  CO2  factor  in  Table  5  is 
the  ratio  of  the  IPCC  (1990)  value  of  the  CO2  forcing  (4.3668  W/m2)  to  the  CO2  forcing 
obtained  in  this  paper  (Column  2,  Table  5).  This  factor  must  be  applied  to  the  CH4  and 
O3  forcing  results  to  make  them  compatible  with  the  IPCC  (1990)  CO2  formulation  for 
the  GWP  calculations  to  be  discussed  later.  In  addition  to  this  factor,  the  results  shown  in 
Table  4  indicate  a  latitude  variation  of  the  direct  CH4  and  indirect  O3  radiative  forcing. 
The  CH4  and  O3  factors  shown  in  Table  5  are  the  product  of  1,  the  ratio  of  the  value  of 
the  CH4  and  O3  tropospheric  radiative  forcings  at  the  various  latitudes  to  the  global 
average  values  (Table  4),  and  2,  the  CO2  factor  of  Table  5.  The  variation  of  these  factors 
represent  a  rough  estimate  of  the  latitude  variation  of  the  radiative  forcing  part  of  the 
GWP  calculation.  The  results  indicate  roughly  a  ±  10  percent  variation  of  the  CH4  GWP 
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with  latitude  and  a  larger  effect  for  the  O3  forcing  of  about  50  percent  at  the  high 
latitudes. 

Table  6  shows  the  latitude  variation  of  the  H2O  radiative  forcing  due  to  changes  in 
the  stratospheric  H2O  as  well  as  the  latitude  averaged  value  of  the  radiative  forcing.  Since 
H2O  is  a  very  strong  greenhouse  gas,  the  tropopause  must  be  carefully  defined  since  the 
H2O  tropospheric  radiative  forcing  will  strongly  depend  on  the  height  of  the  tropopause. 
The  lower  the  tropopause  the  smaller  the  radiative  forcing.  In  order  to  compare  the 
latitude  averaged  radiative  forcing  for  H2O  with  the  radiative  forcing  in  a  globally 
averaged  atmosphere  a  test  calculation  of  the  radiative  forcing  in  the  global  averaged 
atmosphere  at  year  10  was  made.  The  tropospheric  radiative  forcing  for  H2O  in  the 
globally  averaged  atmosphere  (166  mb  tropopause)  was  approximately  30  percent  less 
than  the  globally  averaged  radiative  forcing  for  year  10  given  in  Table  6.  Since  the 
location  of  the  tropopause  is  a  critical  factor  in  the  radiative  forcing  it  is  recommended, 
for  the  case  of  H2O  forcing,  that  the  radiative  forcing  should  be  calculated  as  a  function 
of  latitude  and  then  averaged  as  opposed  to  the  calculation  of  the  radiative  forcing  in  the 
globally  averaged  atmosphere. 

B.  UV  Effects  On  O3  Radiative  Forcing 

Since  O3  has  an  absorption  band  in  the  UV  spectral  region  the  possibility  exists  that 
there  may  be  a  significant  contribution  to  the  total  radiative  forcing  by  the  differential  UV 
absorption.  We  make  an  estimation  of  the  UV  radiation  forcing  contribution  in  the 
following  manner.  The  atmosphere  is  considered  as  a  two  component  system.  The 
stratosphere  is  taken  to  be  a  pure  absorption  region  while  the  troposphere  is  taken  to  be  a 
scattering  and  absorbing  layer.  The  spectral  region  over  which  the  UV  radiative  forcing 
will  be  calculated  is  taken  to  be  250-340  nm  in  5  nm  spectral  sub  intervals.  The  O3 
absorption  cross  sections  and  Rayleigh  scattering  cross  sections  are  given  by  WMO-16 
(1985).  The  radiative  forcing  calculation  is  done  for  the  year  2  model  atmosphere  (the 
maximum  O3  forcing  model).  The  downward  UV  flux  at  the  tropopause  is  calculated 
using  a  simple  Beer’s  Law  exponential  attenuation  of  the  zero  air  mass  solar  flux  (Goody 
and  Yung,  1989).  The  upward  flux  at  the  tropopause  is  calculated  using  a  delta-Eddington 
scattering  model  in  the  troposphere  which  takes  the  troposphere  as  a  single  scattering 
layer.  In  addition  to  the  diffuse,  reflected  flux  component  from  the  troposphere  additional 
upward  flux  components  due  to  the  attenuated  direct  and  downward  scattered  beams  that 
are  reflected  from  the  ground  are  also  considered.  The  ground  albedo  was  taken  to  be  0.2. 
The  UV  radiative  forcing  calculated  in  this  way  amounted  to  about  9  percent  of  the  IR 
radiative  forcing.  Thus  UV  forcing  effects  should  be  a  10  percent  effect  at  most. 

IV.  Carbon  Cycle  Model  and  CO2  Response  Curves. 

Unlike  many  of  the  greenhouse  gases  emitted  to  the  atmosphere,  CO2  does  not  have 
a  chemical  or  photochemical  sink  within  the  troposphere  or  stratosphere.  Removal  of 
CO2  is  therefore  dependent  on  exchanges  with  other  carbon  reservoirs.  Following  the 
emissions  of  a  hypothetical  pulse  of  CC)2  into  the  atmosphere,  about  10%  will  be 
removed  within  2  to  3  years  through  gaseous  diffusion  into  the  mixed  layer  of  the  ocean. 
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Subsequent  removal  is  dependent  on  downward  mixing  into  the  deeper  ocean.  Ultimately, 
the  rest  of  the  initial  pulse  can  be  removed  by  mixing  into  the  ocean,  but  only  over  a  time 
period  of  hundreds  of  years  (Maier-Reimer  and  Hasselmann,  1987).  To  avoid  the  need  to 
use  a  single  lifetime  for  CO2,  we  have  used  a  coupled  climate-carbon  cycle  model  of  Jain 
and  Bach  (1994)  to  calculate  the  decay  rate  of  CO2.  This  model  contains  a  box-diffusion 
carbon  cycle  model  consisting  of  three  reservoirs,  namely  the  atmosphere,  the  mixed 
ocean  layer  (ca.  75  m)  and  the  deep  ocean  (ca.  4000  m).  The  model  ocean  is  treated  as  a 
diffusive  medium  with  constant  vertical  eddy  diffusivity.  The  model  also  takes  into 
account  the  interaction  with  the  biosphere.  The  biospheric  emissions  are  calculated  by  a 
multi-box  globally  aggregated  terrestrial  biosphere  submodel  developed  by  Harvey 
(1989).  The  model  includes  stimulation  of  net  photosynthesis  or  net  primary  productivity 
by  enhanced  atmospheric  CO2  using  the  fertilization  factor  B. 

To  correctly  estimate  the  GWP  values  and  the  past  and  future  €02  concentrations, 
it  is  important  to  use  carbon  cycle  models  that  properly  balance  the  carbon  budget.  To 
balance  the  carbon  cycle  means  to  remove  any  inconsistencies  between  model  output  and 
observations.  When  a  balanced  model  is  run  from  pre-industrial  times  to  the  current  date 
using  the  observed  emissions,  it  will  give  the  correct  history  or  CO2  concentration 
changes.  However,  the  carbon  cycle  model  used  to  calculate  the  GWPs  given  by  IPCC 
(1990, 1992)  does  not  properly  balance  the  carbon  budget  In  a  separate  study  (Wuebbles 
et  al.,  1993)  we  have  assessed  the  effects  of  balancing  the  carbon  budget  on  GWP  values. 
The  CO2  variations  can  be  made  to  fit  observations  by  assuming  the  larger  ocean  flux, 
and/or  accounting  for  feedbacks  in  the  terrestrial  component  of  the  carbon  cycle  and/or 
assuming  the  existence  of  an  additional  CO2  sink.  The  present  model  is  sufficiently 
flexible  that  all  three  methods  may  be  used;  we  chose  to  use  only  the  first  two  methods. 
In  Figure  7,  the  observed  002  concentrations  are  compared  with  the  results  of  the  carbon 
cycle  model.  The  carbon  cycle  model  including  the  biospheric  feedbacks  satisfactorily 
reproduces  the  observed  atmospheric  CO2  concentration  (balanced  case).  Excluding  the 
effects  of  the  biosphere  (i.e.  unbalanced  case)  results  in  an  overestimation  of  CO2 
concentration.  The  time  variation  of  the  fraction  of  CO2  concentration  remaining  after  an 
atmospheric  impulse  is  shown  in  Figure  8.  It  will  be  used  in  the  GWP  calculation 
described  in  the  following  section. 

V.  Global  Warming  Potential  Model 

A.  Basic  Model 

Global  warming  potentials  (GWPs)  are  evaluated  according  to  the  concept  of  IPCC 
(1990).  A  concentration  change  of  trace  gas  X,  arising  from  an  impulse  of  a  unit  mass  of 
the  trace  gas  X  (for  direct  GWPs)  or  from  chemical  processing  due  to  an  impulse  of  gas 
Y  (for  indirect  GWPs),  in  the  atmosphere  produces  a  tropospheric  radiative  forcing 
change  Af^  (t),  which  decreases  with  time  as  X  is  removed  from  the  atmosphere.  This 
radiative  forcing  change  is  considered  as  directly  proportional  to  the  temperature  change 
at  the  surface  via  a  climate  sensitivity  factor.  Commitment  to  radiative  forcing  from  a  unit 
mass  of  X  is  the  integral  of  Afx  (t)  from  the  time  of  an  impulse  to  the  user-selected  GWP 
time  horizon,  T.  The  GWP  of  X  is  the  radiative  forcing  commitment  from  a  unit  mass  of 
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X  divided  by  the  radiative  forcing  commitment  from  a  unit  mass  of  CC)2,  or  a  CO2  like 
molecule^ 

T 

GWPx  =  r - 

0 

The  IPCC  (1990)-derived  GWPs  assume  that  the  background  atmosphere  is  in 
equilibrium  (that  background  concentrations  are  constant  over  the  GWP  integration 
period).  In  our  calculations,  the  chemical  processes  are  allowed  to  operate  on  the 
background  concentrations  for  the  trace  gases  that  are  present  in  the  atmosphere  just 
before  the  impulse.  This  has,  as  will  be  discussed  later,  a  significant  effect  on  the 
equivalent  lifetime  of  the  decay  of  the  impulse  of  the  trace  gas  X  to  the  atmosphere. 

The  definition  of  GWP  (1990)  requires  integration  of  radiative  forcing  (or 
tropospheric  temperature  change)  resulting  from  a  unit  mass  release  of  trace  gas,  usually 
one  kilogram,  and  of  CO2.  A  small  change  in  atmospheric  concentration  resulting  from 
such  a  small  release  can  produce  computational  errors,  so  larger  mass  releases  are  usu^y 
applied  and  then  the  radiative  forcing  is  scaled  to  a  one  kilogram  release.  Radiative 
forcing  scales  linearly  for  many  gases  at  reasonable  concentrations  (IPCC,  1990);  only 
for  releases  of  gas  much  greater  than  one  gigaton  of  CO2,  CH4,  or  N2O  does  nonlinearity 
occur  (Wuebbles  et  al,  1993).  The  radiative  forcing  of  CO2  is  obtained  from  the  IPCC 
(1990)  formula. 


AF, 


CO, 


=  6.31oge 


^0 


(3) 


for  which  Cqo  (0  is  the  concentration  of  CO2  at  time  t  and  is  the  concentration  of 
CO2  which  would  exist  in  the  absence  of  the  pulse.  In  the  equilibrium  case  used  here, 
cco  is  constant  with  respect  to  time,  so  that  the  argument  of  the  log,  is  replaced  by  one 
plus^  the  impulse  response  function  for  CO2,  r'cOj  W* 


AFcOj  =  6.31oge(l+ico,  W); 


ico,(0- 


^CO,(^)~^CQ, 

^co. 


(4) 


The  impulse  response  function  used  in  this  study  was  discussed  in  the  previous 
section  (Figure  8).  Mass-scaled,  integrated  radiative  forcing  for  CO2  is  then  calculated  by 
the  procedure  described  in  Wuebbles  et  al.  (1993)  at  time  horizons  of  20,  50,  100,  200, 
and  500  yr. 
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B.  Extrapolation  and  integration  of  radiative  forcing 

The  radiative  forcing  values  of  Tables  3  and  6  must  be  fitted  before  numerical 
integration  in  Equation  2  due  to  the  uneven  timestep  required  by  computer  limitations  and 
to  time  horizons  exceeding  the  longest  time  after  pulse  for  which  the  atmospheric  and 
line  models  were  run.  IPCC  (1990)  states  that  CH4  radiative  forcings  are  proportional  to 
ccH  while  tropospheric  O3  radiative  forcings  are  linear  in  cq  (t).  In  turn,  the 

concentration  perturbations  of  the  two  trace  species  may  be  taken  as  negative 
exponentials  in  time  L  However,  as  seen  by  Figures  9  and  10,  the  initially  perturbed 
atmosphere  produces  concentrations  inconsistent  with  a  negative  exponential  in  t.  Points 
for  which  this  is  the  case,  the  first  year  of  the  CH4  radiative  forcings  and  years  1-4  after 
the  pulse  for  O3,  are  excluded  from  fitting  for  this  reason.  The  fits  after  weighting  by  CH4 
mass  are: 

A/C//4  =(exp(-58.841±0.028-0.1189±0.0024r)/^^  ,  (5a) 

A/03  =exp(-30.0410±0.0035  -  0.05663±0.000220  ,  (5b) 

LfH20  =  exp(-3.9645  ±  0.040  -  (0.05991 ±  0.0017)r)  .  (5c) 

For  the  integration  of  the  radiative  forcing  of  CH4,  the  value  in  Table  3  is  used  for  the 
first  year  after  the  pulse.  Years  1  to  3  after  the  pulse  also  are  used  in  the  O3  Af 
integration,  and  the  value  at  year  4  is  replaced  by  a  linear  interpolation  of  years  3  and  5. 
Fit  values  are  used  for  all  years  afterward.  Simpson's  Rule  (Press  et  al.,  1992,  equ.  4.1.13) 
is  then  applied  to  integrate  the  radiative  forcing  over  time.  The  radiative  forcing  for  zero 
time  after  the  pulse  is  set  to  zero  in  each  case  for  integration.  Integrated  radiative  forcings 
for  CH4,  H2O,  and  O3  are  each  multiplied  by  a  factor  of  0.8  to  reconcile  the  clear  sky 
approximation  with  the  IPCC  (1990)  cloud  parameterization  as  discussed  previously. 
Table  7  lists  the  integrated  radiative  forcings  per  kilogram  derived  for  each  of  the  gases. 

C.  GWP  Results 

Table  8  shows  the  results  of  the  GWP  calculations.  The  CH4  GWP  values  are 
significantly  larger  than  those  given  by  IPCC  (1990, 1992)  and  Wuebbles  et  al.  (1993).  In 
the  IPCC  calculations  the  CH4  lifetime  was  given  as  10.5  years,  which  implies  that  the 
OH  abundance  was  held  fixed  at  the  ambient  value.  In  the  present  calculation,  the 
chemical  processes  which  affect  CH4,  O3,  and  H2O  abundances  are  allowed  operate  and 
as  a  result  the  OH  abundance  changes  as  the  chemical  equilibrium  is  perturbed  by  the 
CH4  pulse.  Equation  1  predicts  a  lifetime  of  CH4  which  is  dependent  on  the  OH 
abundance.  This  long  apparent  lifetime  of  the  pulse  is  the  result  of  the  slower  destruction 
of  all  CH4  from  the  decrease  in  OH.  The  radiative  forcing  fits  for  CH4  given  by  Equation 
5a  and  shown  in  Figure  9  indicate  an  apparent  lifetime  on  the  order  of  17  years  for  the 
decay  of  the  CH4  pulse.  Wuebbles  et  al.  (1992)  showed  that  when  the  CH4  lifetime  is 
increased,  the  GWP  values  also  increase  towards  the  results  obtained  in  the  present 
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calculation.  For  a  lifetime  of  12.6  years,  GWP’s  at  20,  100,  and  500  years  of  41.7,  15.9, 
and  7.1  were  obtained.  A  simple  extrapolation  of  the  GWP  vs.  CH4  lifetime  results 
Wuebbles  et  al.  (1992)  predicts  a  GWP  of  approximately  52  at  20  years  which  is  the  same 
as  that  obtained  in  this  paper.  The  results  obtained  for  the  O3  indirect  radiative  forcing  are 
on  the  order  of  10  to  20  percent  higher  those  given  in  EPCC  (1990)  and  Wuebbles  et  al. 
(1993).  Given  the  uncertainties  of  the  calculations,  at  least  20  to  30  percent,  this  is  a 
reasonable  agreement.  The  indirect  GWP’s  for  stratospheric  H2O  obtained  in  this  paper 
are  lower  by  factors  of  3  to  5  than  those  given  in  IPCC  (1990).  This  is  probably  due  tc  an 
oveiestimation  of  the  tropospheric  radiative  forcing  of  H2O  in  IPCC  (1990).  The  ratio  of 
the  indirect  H2O  GWP  to  the  direct  CH4  GWP,  for  the  calculations  in  this  paper  is  about 
4  percent  and  agrees  quite  well  with  the  results  obtained  by  Lelieveld  and  Crutzen  (1992) 
who  give  a  value  of  5  percent  for  this  ratio.  This  effect  is  strongly  dependent  on  the 
height  of  the  tropopause  since  H2O  is  such  a  strong  IR  absorber  and  should  have  a  sharp 
mixing  ratio  gradient  at  the  tropopause.  The  overall  uncertainty  in  the  GWP  calculations 
is  due  to  the  uncertainties  in  the  lifetime  of  CH4,  the  amount  of  NOx  present  in  the 
atmosphere,  and  the  radiative  forcing  calculations.  The  best  estimate  of  the  overall 
uncertainty  in  the  calculations  is  on  the  order  of  20  percent. 


VI.  Conclusions 

We  have  calculated  the  direct  CH4  GWP  and  the  indirect  O3  and  stratospheric  H2O 
GWPs  resulting  from  a  CH4  pulse  equal  to  -36  percent  of  the  total  atmospheric  mass  of 
CH4.  The  direct  CH4  GWP  is  significantly  larger  than  previous  estimates  as  a  result  of 
the  effect  on  the  pulse  response  time  due  to  the  variation  of  the  OH  abundance.  The 
agreement  of  the  indirect  O3  G^^P  effect  with  previous  estimates  is  within  the  uncertainty 
bounds  of  the  calculation.  The  indirect  H2O  GWP  is  much  less  than  the  estimate  given  in 
IPCC  (1990,  1992)  but  agrees  very  well  with  more  recent  calculations  Lelieveld  and 
Crutzen  (1992).  The  H2O  effect  is  very  strongly  dependent  on  the  location  of  the 
tropopause  because  of  the  strong  IR  absorption  of  H2O.  Future  calculations  need  to  be 
made  to  determine  the  effect  of  the  magnitude  of  the  CH4  impulse  on  the  impulse  decay 
time  and  the  effect  of  the  presence  of  hydrocarbons  and  NOx  on  the  direct  and  indirect 

GWP’s. 
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Figure  Captions 

Figure  1.  Globally  and  annually  averaged  profiles  of  water  vapor,  ozone,  and  methane  as 
a  function  of  altitude  for  the  ambient  atmosphere. 

Figure  2.  Pressure-temperature  profile  for  the  ambient  atmosphere.  The  temperatures  are 
globally  and  annually  averaged. 

Figure  3.  Total  mass  of  methane  (in  units  of  teragrams)  in  the  model  atmosphere  as  a 
function  of  time  before  and  after  a  pulse  emission  of  methane  equal  to  the 
equilibrium  increase  of  CH4  produced  by  a  constant  25%  increase  in  total 
emissions.  Feedbacks  in  the  lifetime  of  methane  result  in  about  a  36%  peak  increase 
in  total  atmospheric  methane.  The  ambient  value  of  the  methane  mass  is  indicated 
by  the  dashed  line. 

Figure  4.  Annual  average  total  ozone  column  after  the  pulse  increase  of  methane.  The 
ambient  value  of  the  ozone  column  is  indicated  by  the  dashed  line. 

Figure  5.  Annually  averaged  tropospheric  ozone  column  after  the  pulse  increase  of 
methane.  The  ambient  value  of  the  troposperic  ozone  column  is  indicated  by  the 
dashed  line. 

Figure  6.  Increase  in  annually  averaged  stratospheric  water  vapor  burden  (kg)  before  and 
after  the  pulse  emission  of  methane. 

Figure  7.  The  time  variation  of  the  observed  CO2  concentration  (+),  the  carbon  cycle 

model  prediction  including  biospheric  feedbacks  ( _ ),  and  the  carbon  cycle 

model  prediction  excluding  the  biosphere  ( - ). 

Figure  8.  The  response  of  atmospheric  €02  to  a  pulse  of  1  GTC  (gigatons  of  Carbon). 

Figure  9.  Radiative  forcing  of  methane  vs.  time  (years)  after  the  pulse  emission  of 
methane.  The  correlated  k-distribution  radiative  forcing  model  calculations  are 
denoted  by  the  symbols.  The  solid  line  is  the  fit  to  the  model  calculations  given  by 
Equation  4. 

Figure  10.  Radiative  forcing  of  ozone  vs.  time  (years)  after  the  pulse  emission  of 
methane.  The  correlated  k-distribution  radiative  forcing  model  calculations  are 
denoted  by  the  symbols.  The  solid  line  is  the  fit  to  the  model  calculations  given  by 
Equation  4. 
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Table  1.  Chemical  reaction  mechanism  for  LLNL  two-dimensional  model.  The 
reference  used  for  all  photochemical  data  is  Demore  et  al.  (1992). 


(a)  Two-body  reactions  (Arrhenius  rate  constant  expression  k  =  A- exp 


>. 


No. 

Reaction 

A* 

E/R2 

1 

O+O3  ->2*02 

8.00E-12 

2,060 

2 

OiD  +  Nj  0  +  N2 

1.80E-11 

-110 

3 

OjD  +  O2  0  +  Q2 

3.20E-11 

-70 

4 

OiD  +  O3  ->  2*02 

1.20E-10 

0 

5 

OjO  +  03  — >  Qj  •+  2*0 

1.20E-10 

0 

6 

H20  +  0iD-»2*0H 

2.20E-10 

0 

7 

H2  +  OiD^OH  +  H 

l.OOE-10 

0 

8 

N2O  +  OjD  — >  N2  +  O2 

4.90E-11 

0 

9 

N20  +  0iD->2*N0 

6.70E-11 

0 

10 

CH4  +  0iD->CH20  +  H2 

1.40E-11 

0 

11 

CH4  +  OiD  CH3O2  +  OH 

1.40E-10 

0 

12 

CF2a2+OiD->2*a 

1.12E-10 

0 

13 

CF3CCI3  +  OiD  ^  3*C1 

1.60E-10 

0 

14 

CF3CFa2  +  OiD^2*Cl 

128E-10 

0 

15 

CF3CF2CI  +  OiD  ^  Cl 

7.12E-11 

0 

16 

CHF2C1  +  OiD-»C1 

7.00E-11 

0 

17 

H  +  O3  ->0H  +  C)2 

1.40E-10 

470 

18 

H2  +  0H-^H20  +  H 

5.50E-12 

2,000 

19 

OH  +  O3  ^  HO2  +  O2 

1.60E-12 

940 

20 

0H  +  0-»02+H 

2.20E-11 

-120 

21 

H02  +  0^0H  +  02 

3.00E-11 

-200 

22 

HO2  +  O3  ->0H  +  2*02 

l.lOE-14 

500 

23 

H02+0H^H20  +  02 

4.80E-11 

-250 

24 

H202  +  0H^H20  +  H02 

2.90E-12 

160 

25 

N+02-4NO+0 

1.50E-11 

3,600 

26 

N  +  N0-»N2  +  0 

3.40E-11 

0 

27 

NO  +  O3  ^NQ2  +  02 

2.00E-12 

1,400 

28 

N0  +  H02->NQ2  +  6H 

3.70E-12 

-250 

29 

N02+0-^N0  +  02 

6.50E-12 

-120 

30 

N02  +  03  — >  N03  +  02 

1.20E-13 

2,450 

31 

N2O5  +  H2O  ^  2*HN03 

5.00E-22 

0 

32 

HO2NO2  +  OH  ^  H2O  +  NO2  +  O2 

1.30E-12 

-380 

33 

Cl  +  O3  — ^  CIO  +  O2 

2.90E-11 

260 

34 

Cl  +  H02->Ha  +  02 

1.80E-11 

-170 

35 

a  +  H02-^0H  +  C10 

4.10E-11 

450 

36 

CIO  +  0  Cl  +  O2 

3.00E-11 

-70 

37 

C10  +  0H->H02  +  C1 

l.lOE-11 

-120 

38 

C10  +  H02^Q2  +  H0C1 

4.80E-13 

-700 

39 

C10  +  N0-»N02  +  C1 

6.40E-12 

-290 

40 

C10  +  C10-»Cl  +  C100 

4.00E-13 

1,250 

41 

CIO  +  CIO  — >  CI2  +  O2 

4.00E-13 

1,250 

^  Units  are  cm^/faioiecule~sec. 
^Unit  is  degrees  Kelvin. 
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No. 

Reaction 

A3 

E/R^ 

42 

cx:k)  +  ci^2*cio 

3.40E-11 

-160 

43 

0C10  + 0-4  02  +  00 

2.50E-12 

950 

44 

OCIO  +  OH  ^  HOCl  +  O2 

4.50E-13 

-800 

45 

OCIO  +  NO-4NO2  +  CIO 

2.50E-12 

600 

46 

CIOO  +  CI-+CI2  +  O2 

2.30E-10 

0 

47 

ci202+ci-4Ci2+aoo 

l.OOE-10 

0 

48 

Ha  +  0H-+H20  +  CI 

2.60E-12 

350 

49 

H0C1  +  0H^H20  +  C10 

3.00E-12 

500 

50 

CIONO2  +  0  -+  CIO  +  NO  +  O2 

3.19E-13 

800 

51 

CIONO2  +  0  -+  CIO  +  NO2  +  0 

2.58E-12 

800 

52 

CIONO2  +  OH  ^  HOCl  +  NO  +  O2 

1.32E-13 

333 

53 

CIONO2  +  OH  ->  HOCl  +  NO2  +  0 

1.07E-12 

333 

54 

CIONO2  +  Cl  -4  Ch  +  NO  +  O2 

7.48E-13 

-160 

55 

CIONO2  +  Cl  -4  Cfc  +  NO2  +  0 

6.05E-12 

-160 

56 

Br  +  03  -4B1O  +  O2 

1.70E-11 

800 

57 

Br  +  H02^HBr  +  02 

1.50E-11 

600 

58 

Br  +  N03  ^BrO  +  N02 

1.60E-11 

0 

59 

Br  +  0C10-4Bi0  +  C10 

2.60E-11 

1,300 

60 

Br  +  CI2O2  ^  BrCl  +  ClOO 

3.00E-12 

0 

61 

Br  +  CH2O  ->  HBr  +  HO2+  CO 

1.70E-11 

800 

62 

BtO  +  0  -4  Br  +  O2 

3.00E-11 

0 

63 

BiO  +  H02^HOBr  +  Q2 

6.20E-12 

-500 

64 

Br0  +  0H^Br  +  H02 

l.lOE-11 

0 

65 

Br0  +  N0-4Br  +  N02 

8.80E-12 

-260 

66 

BrO  +  NO3  —4  NC^  +  Br  +  O2 

l.OOE-12 

0 

67 

BrO  +  CIO  -+  Br  +  OCIO 

1.60E-12 

^30 

68 

BrO  +  CIO -»Br+ ClOO 

2.90E-12 

-220 

69 

BrO  +  C10-4Bia  +  02 

5.80E-13 

-170 

70 

BrO  +  BtO->2*Br  +  02 

1.40E-12 

-150 

71 

HBr  +  0-4Br  +  OH 

5.80E-12 

1,500 

72 

HBr  +  0H->Br+H20 

l.lOE-11 

0 

73 

CH4  +  OH  -4  CH3O2  +  H2O 

2.90E-12 

1,820 

74 

CH2O  +  OH  -4  H2O  +  HO2  +  CO 

l.OOE-11 

0 

75 

CH302  +  0^CH20  +  H02 

3.(X)E-11 

0 

76 

CH2O  +  0 ->  HO2  +  OH  +  CO 

3.40E-11 

1,600 

77 

CH3O2  +  NO  -4  HO2  +  CH2O  +  NO2 

4.20E-12 

-180 

78 

CI  +  CH4  ^CH302  +  HC1 

l.lOE-11 

1,400 

79 

Cl  +  CH2O  -4  HCl  +  HO2  +  CO 

8.10E-11 

30 

80 

CH3CI  +  OH  ^  Cl  +  H2O  +  HO2 

2.10E-12 

1,150 

81 

CH3CCI3  +  OH  ^  H2O  +  3*C1 

1.80E-12 

1,550 

82 

CHF2C1  +  0H^C1  +  H20 

1.20E-12 

1,650 

83 

CH3CI  +  Cl  -4  HO2  +  CO  +  2*HC1 

3.30E-11 

1,250 

84 

CH3Br  +  0H->Br  +  H20 

3.60E-12 

1,430 

^ Units  are  cm^/tnolecule-sec. 
^Unit  is  degrees  Kelvin, 


(b)  Three-body  reactions  (special  rate  constant  expressions  such  as  Troe  equation) 


Number 

Reaction 

1 

0  +  02  Os 

2 

H  +  02  ->  HO2 

3 

H02  +  H02^H202+02 

4 

HO2  +  HO2  +  H20^  H2O2  +  O2  +  H2O 

5 

N02  +  H02^H02N02 

6 

N03  +  N02^N205 

7 

N2O5  -»  NO2  +  NO3 

8 

NO2  +  OH  -4  HNO3 

9 

HNO3  +  OH  ^  H2O  +  NO  +  O2 

10 

HNO3  +  OH  H2O  +  NO2  +  0 

11 

H02N02->H02  +  N02 

12 

C10  +  N02->C10N02 

13 

C10  +  C10^Cl202 

14 

ClOO  Cl  -1-  02 

15 

CI2O2  ^  2*C10 

16 

BrO  +  NO2  ^  BrON02 

17 

CO  +  OH^H 

18 

N2O5  2*HN03 

19 

CIONO2  -»  HOCl  +  HNO3 
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(c)  Photolysis  reactions 


Number  I _ Reaction 


1. 

02  +  hv->2*0 

2. 

03  +  hv->0  +  02 

3. 

03  +  hv-^0iD  +  02 

4. 

H2O2  +  hv  ->2*OH 

5. 

H20  +  hv^H  +  0H 

6. 

N02  +  hv^N0  +  0 

7. 

N2O  +  hv  -»N2  +  OiD 

8. 

NO  +  hv-^N+0 

9. 

NO3  +  hv  ^N02  +  0 

10. 

NO3  +  hv  ->NO  +  O2 

11. 

N2O5  +  hv  -^N02  +  no  +  O2 

12. 

N2O5  +  hv  ->2*N02  +  0 

13. 

N2O5  +  hv  ^N02  +  no  +  2*0 

14. 

N2O5  +  hv  ^2*NO  +  O2  +  0 

15. 

HNO3  +  hv  ^OH  +  NO2 

16. 

HO2NO2  +  hv  ->OH  +  NO  +  O2 

17. 

HO2NO2  +  hv  OH  +  NO2  +  0 

18. 

HO2NO2  +  hv  — ^H02  +  NO2 

19. 

Cl2  +  hv-»2*Cl 

20. 

CIO  +  hv  — >  Cl  +  0 

21. 

C10  +  hv-4Cl  +  0iD 

22. 

0C10  +  hv->0  +  C10 

23. 

aoo  +  hv  -4  CIO  +  0 

24. 

0202  +  hv  -»Q  +  ClOO 

25. 

HOCl  +  hv  ^  OH  +  Cl 

26. 

CIONO2  +  hv -^a  +  NO  +  O2 

27. 

C10NO2 + hv  -^a  +  NO2  +  0 

28. 

aoN02  +  hv  ->a  +  NO2  +  0 

29. 

Bid  +  hv  — >  Br  +  Cl 

30. 

Br0  +  hv-»Br  +  0 

31. 

HBr  +  hv->H  +  Br 

32. 

HOBr  +  hv  Br  +  OH 

33. 

Br0N02  +  hv  ->Br  +  NO2  +  0 

34. 

Br0N02  +  hv  ->Br  +  NO  +  O2 

35. 

CH2O  +  hv  ->C0  +  H2 

36. 

CH2O  +  hv  ->H02  +  CO  +  H 

37. 

CH3CI  +  hv  ^CH302  +  d 

38. 

CC14  +  hv^4*d 

39. 

CH3CCl3  +  hv^3*Cl 

40. 

CFCI3  +  hv ->3*d 

41. 

CF2CI2  +  hv  ^2*C1 

42. 

CF3CCl3  +  hv->3*Cl 

43. 

CF3CFCI2  +  hv  ^2*C1 

44. 

CF3CF2CI  +  hv  -^Cl 

45. 

CHF2Cl  +  hv-^Cl 

46. 

CH3Br  +  hv  ->Br  +  CH3O2 

47. 

CF3Br  +  hv  — >Br 

48. 

CF9dBr  +  hv^Br  +  Cl 

;  2.  Surface  emission  inventory  for  the  LLNL  two-dimensional  model,  (teragrams 
per  year) 


Model  Zone 

Latitude 

Center 

Methane 

Carbon 

Monoxide 

Nitrogen 

Dioxide 

1 

-80.1 

0.00 

0.60 

0.00 

2 

-70.9 

0.00 

1.00 

0.00 

3 

-62.3 

0.00 

8.45 

0.00 

4 

-53.1 

0.15 

18.64 

1.98 

5 

^2.5 

4.08 

35.62 

7.27 

6 

-30.0 

18.99 

84.83 

10.35 

7 

-17.5 

24.67 

136.09 

15.63 

S 

-5.7 

39.25 

154.50 

11.67 

9 

5.7 

42.43 

145.11 

12.67 

10 

17.5 

54.99 

167.10 

10.69 

11 

30.0 

81.71 

197.88 

20.37 

12 

42.5 

70.41 

209.05 

27.09 

13 

53.1 

57.51 

126.63 

8.88 

14 

62.3 

30.00 

34.40 

4.83 

15 

70.9 

15.92 

5.02 

0.00 

16 

80.1 

12.77 

3.01 

0.00 

TOTALS 

(Tg/yr): 

452.9 

1328 

131.4 
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Table  3.  Radiative  Forcing  W/m^ 


Year 

Direct  CH4 

Indirect  O3 

C02(1.29) 

CO2  (2.0) 

1 

0.2500 

0.0947 

2 

0.2455 

0.1180 

1.956 

5.453 

3 

0.2312 

0.1165 

6 

0.1950 

0.1013 

7 

0.1844 

0.0957 

8 

0.1744 

0.0909 

9 

0.1649 

0.0861 

10 

0.1558 

0.0818 

15 

0.1176 

0.0614 

20 

0.0884 

0.0462 

25 

0.0665 

0.0348 

30 

0.0500 

0.0261 
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Table  4.  Latitude  Variation  Of  Radiative  Forcing  (W/M^)  For  Year  2 


LATITUDE  (DEG)  DIRECT  CH4 

(P.  TROP.) 

INDIRECT  O3 

(P.  TROP.) 

-71 

0.2403  (292mb) 

0.0625  (292mb) 

-30 

0.2755  (139mb) 

0.1518  (139mb) 

+5.7 

0.2628  (123mb) 

0.1516  (123mb) 

+30 

0.2659  (135mb) 

0.1406  (135mb) 

+71 

0.2244  (302mb) 

0.0603  (302mb) 

Avg. 

0.2603 

0.1297 

%Diff. 

+5  (0.2455) 

+9(0.1180) 

Table  5.  Latitude  Effect  Of  CO2  Forcing  For  An  Impulse  Doubling  Of  CO2 


LATITUDE 

FORCING  (W/m2) 

CO2 FACTOR 

CH4 

FACTOR 

O3 

FACTOR 

Average 

5.453 

0.8 

•  0.8 

0.8 

-71 

4.418 

0.988 

0.967 

0.523 

-30 

5.735 

0.76 

0.852 

0.977 

+5.7 

5.976 

0.73 

0.781 

0.939 

+30 

5.773 

0.76 

0.823 

0.906 

+71 

4.624 

0.944 

0.863 

0.482 
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Table  6.  Latitude  Variation  Of  H2O  Radiative  Forcing  (W/M^) 


LATITUDE  (DEG)  H2O  Forcing 

Years  Year  10  Year  20 


1 - 

-60(278mb) 

2.085  (-3) 

5.302  (-3) 

3.761  (-3) 

-30(139mb) 

4.179  (-3) 

1.021  (-3) 

7.270  (-3) 

Eq.(123mb) 

2.995  (-3) 

6.456  (-3) 

4.715  (-3) 

+30(135mb) 

5.295  (-3) 

1.251  (-2) 

8.911  (-3) 

+60  (278mb) 

1.741  (-3) 

4.032  (-3) 

2.932  (-3) 

Avg. 

3.629  (-3) 

8.535  (-3) 

6.099  (-3) 

803 


C3RCD193-80.ab 


Table  7.  Radiative  Forcing  Integrals  per  Kilogram  Mass  of  the  Trace  Gases 


20  yr 

50  yr 

100  yr 

200  yr 

500  yr 

CH4 

LSOSe-l^* 

2.074e-12 

2.183e-12 

b 

b 

03(CH4) 

7.93  le- 13 

1.127e-12 

1.198e-12 

b 

b 

H20(CH4) 

5.624e-14 

1.012e-13 

1.016e-13 

C02 

3.117e-14 

6.303e-14 

1.073e-13 

L832e-13 

3.612e-13 

“Read  "e-”  as  an  exponent,  thus  1 .505e-12  =  1 .505  x  lO  ”  Units  are  W  yr  m'*  kg  '  at  the  tropopause. 
'’Taken  as  identical  to  the  value  for  time  horizon  of  100  yr. 


Table  8.  Global  Warming  Potentials  and  the  Indirect  GWP  Total  Uncertainty  Estimate 


Sources: 

(a)  Wuebbles  et  al.,  1993 

(b)  IPCC,  1990 

(c) IPCC,  1992. 
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Figure  4  Annual  average  total  ozone  column  after  the  pulse  increase  of  methane. 
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Annually  averaged  tropospheric  ozone  column  after  the  pulse  increase  of  methane. 


4.00E+10 


810 


Tim©  aftor  puis©  (y©ars) 


CO2  Concentration 


811 


1.0  s 


o 


6u|uieiuaj  'OD  uojjoBJj 


812 


100  200  300  400  500 

Time  after  Impulse/years 


A)  Methane 


813 


17th  Annual  Renew  Conference  on  Atmospheric  Radiation  Models 


7-8  June  1994 


Name 

Affiliation 

Phone  Number 

Mr.  Leonard  Abreu 

ONTAR  Corporadon 

508-689-9622 

Prabhat  Acharya 

Spectral  Sciences,  Inc. 

617-270-5982 

Steven  Adler-Golden 

Spectral  Sciences,  Inc. 

617-273-4770 

Ronald  E.  Alley 

Jet  Propulsioa  ILaboratoty 

518-354-0751 

Don  Anderson 

JHU/APL 

Gail  P.  Anderson 

PL/GPOS 

617-377-3687 

Peter  Armstrong 

Stewart  Radiance  Laboratory 

617-377-3626 

Robert  Bcland 

PL/GPOS 

617-377-4604 

Frank  Bennett 

Mitre  Corporation 

617-271-8712 

Alexander  Berk 

Spectral  Sciences,  Inc. 

617-273^770 

Larry  Bernstein 

Spectral  Sciences.  Inc. 

617-273-4770 

Mike  Bossert 

Boeing  Computer  Services 

619-939-0667 

Tyler  Brown 

JFL/NASA 

818-354-6357 

James  T.  Bunting 

PL/GPAB 

617-377-4854 

Johnson  Chan 

ONTAR  Corporation 

508-689-9622 

Kelly  Chance 

Smithsonian 

617^95-7389 

Michael  Cheifetz 

SPARTA 

617-863-1060 

Paola  Ciarpallini 

Smithsonian 

617-495-7497 

George  Cipperly 

Martin  Marrietta  Corporation 

610-531-6266 

F.O.  Clark 

PL 

617-377-4139 

John  Conant 

Aerodyne  Research,  Inc. 

508-663-9500 

814 


17th  Annual  RevieTf  Conference  on  Atmospheric  Radiation  Models 


7-8  June  1994 


Name 

Affiliation 

Phone  Number 

■William  Comette 

PRA 

619-455-9741 

Kathy  Crow 

SciTec,  Inc. 

609-921-3892x248 

Gil  Davidson 

Photometries 

617-935-6500 

Richard  E.  Davis 

NASA-Langley  Research  Center 

804-864-1647 

Adam  D.  Devir 

EORD/Technion/Israel 

972^294233 

Edmond  Dewan 

PL/GPOS 

617-377-4401 

Susan  Downer 

Mission  Research  Corporation 

603-891-0070 

Tom  Dugan 

NRL 

Martin  Fair 

DRA  Famborough  UK 

44-252-395344 

Harry  Finklc 

MTT/LL 

617-981-4905 

John  Feeiey 

Boeing/NAWC 

619-939-0605 

Grace  E.  Gagnon 

PI7GPOS 

617-377-2336 

William  Gallery 

AER 

617-547-6207 

Richard  Garner 

Photometries 

617-935-6500 

Stuart  Gathman 

NRnD 

619-553-1418 

Roy  Goeller 

Los  Alamos  National  Laboratory 

505-665-1107 

Donald  Grantham 

PUGPAA 

617-377-2982 

Allen  Grossman 

Los  Alamos  National  Laboratory 

510-923-6371 

J.A.  Hall 

PIVGPIM 

Lanene  Harada 

WJSA 

703-558-7900 

Sabrina  Harvey 

ONTAR  Corporation 

508-689-9622 

7-8  June  1994 


17th  Annual  Review  Conference  on  Atmospiieric  Rarftgrinn  Models 


Name 

Affiliation 

Phone  Number 

Michael  lacono 

AER 

617-349-2208 

Precila  Ip 

Mission  Research  Corporation 

603-891-0070 

Ken  lucks 

Smithsonian  Ast.  Ops 

617-496-7580 

Simon  Kaelin 

Aerosols/Eth  Zurich 

Brook  Koffinc 

Aerospace  Corporation 

310-336-7412 

Robert  Kuiucz 

SAO 

617-495-7429 

Alex  Kutepov 

University  of  Munich 

49-89-92209430 

Tom  Kyle 

PL/XPG 

617-377-7443 

Olga  Lado-Bordowsky 

BNSSAT 

Steven  Lipson 

PL/GPOS 

617-377-3626 

U.  Makhlouf 

SRL 

617-377-4203 

Steve  Mazuk 

Aerospace  Corporation 

310-336-5614 

Andy  McCann 

ONTAR  Corporation 

508-689-9622 

Sue  McKenzie 

Mission  Research  Laboratory 

603-891-0070 

Joe  Michalsky 

SUNY 

518-442-3809 

Ferenc  Miskolczi 

UMD 

301-799-7693 

MikeNeer 

Scilbc,  Inc. 

609-921-3892x206 

Meg  Noah 

Mission  Research  Laboratory 

LQsa  ODcarinen 

Finnish  Meteorological  Institute 

358-0-1929521 

Clifford  A.  Paiva 

Nswcy 

703-663-4980 

Brunei  Pascal 

Meteo  France 

33-96056765 

816 


7-S  June  1994 


17th  Annual  Review  Conference  on  Atmospheric  Radiation  Models 


Name 

Affiliation 

Phone  Number 

Richard  Picard 

PL/GPOS 

617-377-2222 

Russell  Phillbrick 

ARiypso 

814-863-7682 

Angela  Phillips 

Hughes  Aircraft  Company 

310-616-0124 

Behroo2  Rezai 

801-797-3564 

Lt.  CoL  John  Roadcap 

PITWE 

505-846-4722 

David  Robertson 

Spectral  Sciences,  Inc 

617-273-4770 

Ronald  Rodney 

WPAFB-WrightLab. 

513-255-1978 

Larry  Rothman 

PL/GPOS 

617-377-2336 

John  Schroedcr 

ONTAR  Corporation 

508-689-9622 

Mark  Schmitt 

Los  Alamos  National  Laboratory 

505-667-0655 

John  E.A.  Selby 

Northrop-Grununan 

516-575-6608 

Dima  Seiiverstov 

Mitre  Corporation 

617-377-5314 

Eric  Shettle 

NRL 

202-404-8152 

Richard  Shirkey 

ARLT5ED 

505-678-5470 

Lewis  L.  Smith 

Northrop  Grumman 

516-575-2196 

Roberta  Smith 

US  Armylbcom 

703-704-1762 

Ned  Snell 

AER 

617-547-6207 

Pierre  Simoneau 

ONERA 

33-46734850 

Knut  Stamnes 

University  of  Alaska 

907-474-7368 

Dm  Stevens 

ARL/PSU 

814-863-0851 

■WiUiam  Stoner 

SAIC 

508-667-6365 

7-8  June  1594 


17tli  Annual  Review  Coofermce  on  Atmospheric  Radiation  Models 


Name 

Affiliation 

Phone  Number 

James  Telford 

University  of  Nevada 

702-677-3201 

David  Tofsted 

ARL 

505-678-3039 

Martine  Voltaire 

ONTAR  Corporation 

508-689-9622 

Philip  WaJDker 

Naval  Postgraduate  School 

408-656-2674 

Zhengning  Wan 

University  of  California 

805-893-4541 

Jinxue  Wang 

PUGPOS 

617-377-4854 

Alan  Wetraore 

ARL 

505-678-5567 

Robert  C.  Wilson 

Boston  College 

617-552-8765 

Peter  P.  Wintersteincr 

ARCON 

617890-3330x231 

Rich  Woodford 

USAFETAC 

618-256-5211 

Kouichi  Yoshino 

CFA 

617-495-2796 

Douglas  Youmans 

WJ.  Schafer 

508-256-2070 

Daniel  Zhou 

SDUUSU 

301-953-5470 

818 


AUTHORS 


ABREO,  L.W. 
ACHARYA,  P. 
AHDERSON,  D.E. 
ANDERSON,  G.P. 
ANDERSON,  S.H. 
ADLER-GOLDEN,  S.H. 
ARMSTRONG,  P.S. 
ARMSTRONG,  R. 

BARZILAI,  B. 
BELAND,  R.R. 

BEN  SHALOM,  A. 
BERK,  A. 

BERNSTEIN,  L. 
BEVILAQUA,  R.H. 
BINGHAM,  G.E. 
BLDHBERG,  W.A.H. 
BOISHOT,  A. 
BONZAGNI,  H. 
BRANIMN,  N. 

BROWN,  J.H. 

CHEIFETZ,  H.G. 
CHETWYND,  J.H. 
CLARKE,  F. 

CLOUGH,  S.A. 
CONANT,  J. 
OORNETTE,  W.M. 
CRANSTON,  K.E. 
CROW,  K. 

CRUNCELTON,  J.P. 

DAVIS,  B. 

DE,  P. 

DEBRESTIAN,  D. 
DEHAJISTRE,  R. 
DESCHAHP,  J. 

DEVIR,  A.D. 

DEWAN,  E. 

DODD,  J.A. 

DUFF,  J.W. 

DURAND,  G. 

ESMOND,  J.R. 
ESPADA,  L. 

FANESLW,  J.L. 
FROMM,  H. 

GALLERY,  W.O. 
GATHHAN,  S.G. 
GLACCUH,  W.J. 
GRANT,  K.E. 
GREGOIRE,  G. 
GROSSBARD,  N. 
GROSSMAN,  A.S. 
GRUHNINGER,  J.H. 


Al,  A2,  1,  26 

A2,  A3,  A43,  26,  46,  739 

A12,  230 

Al,  A2,  A3,  A21,  A37,  A43,  A44,  1,  26,  46,  351,  629,  739,  764 

A19,  335 

A33,  557 

A36,  620 

A22 


A15 

A30,  507 
A15 

A2,  A43,  26,  739 

A2,  A33,  A43,  26,  557,  739 

A38,  649 

A37,  A44,  629,  764 

A36,  620 

A9,  162 

A26 

A15 

A28,  A29,  A33,  448,  468,  557 
A17,  320 

Al,  A2,  A43,  1,  26,  46,  739 
A4,  71 

Al,  A42,  1,  726 
26 

A3,  46 
A31,  536 
AlO,  187 
A31,  536 

A13,  247 
A28,  448 
A38,  649 
A12,  230 
162 
A15 

A25,  394 

A27,  A36,  431,  620 
A33,  557 
162 


A20,  346 
304 


A39,  673 
A38,  649 

A21,  A40,  351,  690 
A6,  101 
A38,  649 
A45,  776 
162 

A25,  A29,  394,  468 
A45,  776 

A28,  A33,  448,  557 


81^ 


GUZZI,  R. 

HALL,  L.A. 

HARADA,  L. 

HAYES,  P.B. 

HEALY,  R.J. 
HORNSTEIN,  J.S. 
HUHHEL,  J.R. 

lACONO,  H.J. 

ITO,  K. 

JAIH,  A. 

JONES,  J.R. 
KIMBALL,  L.H. 
KNEIZYS,  F.X. 
KRIGHAN,  S. 
KDHHER,  J.B. 
KDROCZ,  R.L. 
KDTEPOV,  A.A. 
KYROLA,  E. 

LARSEN,  N. 

LESLIE,  D. 

LIPSON,  S.J. 
LWIGTIN,  D.R. 
LOPEZ-POERTAS,  H. 
LOWELL,  J.R. 
LOHPE,  J. 

LYNCH,  D.K. 

HCCOHB,  J. 
MCKENZIE,  S. 
MALHERBE,  C. 
MARAVUDA,  S. 
MATHEW,  H.W. 
MATHEWS,  L.A. 
MATSDI,  T, 

HAZUK,  S. 

MEGILL,  L.R. 
MERGENTHALER,  J.L. 
HICHON,  P. 
MINSCHWANER,  K, 
MISKOLCZI,  F. 
MONCET,  J.L. 

NADILE,  R.H. 

NANDI,  S. 

NEER,  M. 

OGIBALOV,  V.P. 
OIKARINEH,  L. 

PARKER,  J. 
PARKINSON,  W.H. 
PATTEN,  K.O. 
PATTERSC®,  J.E. 
PHILBRia,  C.R. 
PICARD,  R.H. 


A26 

A2 

A14,  276 
A24,  373 
A33,  557 
A38,  649 
A17,  330 

A42,  726 
A20,  346 

A45,  776 
A17,  320 

A2,  A43,  26,  739 
A2,  26 
A38,  649 
A37,  629 
A18,  333 
A34,  580 
705 

All,  202 
A14,  276 
A36,  620 
A17,  320 
A35,  601 
A36,  620 
A38,  649 
A8,  143 

A39,  673 
A22 

A9,  162 
A23,  357 
A33,  557 
A7,  118 
A20,  346 
A8,  143 
A37,  629 
A37,  629 
162 
A2 
A26 

A21,  A40,  1,  351,  690 

A36,  A44,  620,  629,  764 
A39,  673 
AlO,  187 

A34,  580 
A41,  705 

A39,  673 
A20,  346 
A45,  776 
A39,  673 
A23,  357 

A27,  A35,  431,  601 


820 


RE2AI,  B.K. 

A37,  629 

RICKLIH,  J.C, 

A31,  536 

ROBERTSON,  D.C. 

A2,  A3,  A33,  A43,  26,  46,  557,  739 

ROCHE,  A.E. 

A37,  629 

SAVACOOL,  H. 

276 

SHARMA,  R.D. 

A33,  557 

SHETTLE,  E.P. 

A2,  A19,  A38,  26,  335,  649 

SHIRKEY,  R.C. 

A13,  A16,  247,  304 

SPED,  G.H. 

A34,  580 

SIPOLA,  E. 

705 

SIHONEAD,  P. 

A9,  162 

SMITH,  R.W. 

A31,536 

SMELL,  H.E. 

A24,  373 

SPARKES,  L. 

A39,  673 

STAHNES,  K. 

All,  202 

STEED,  A.J. 

A44,  764 

STEVENS,  T.D. 

A23,  357 

SDNDBERG,  R.L. 

A28,  A33,  448,  557 

TAHARESIS,  J.S. 

A45,  776 

TAHHINEN,  J. 

705 

TELFORD,  J.W. 

A32,  546 

THERIADLT,  J.H. 

1 

TOFSTED,  D.H. 

A13,  A16,  247,  304 

TSAY,  S. 

All,  202 

VAIL,  J.J. 

26 

VANZANT,  T. 

A25,  394 

WALKER,  P.L. 

A7,  118 

WANG,  J. 

Al,  1 

WETHORE,  A.E. 

A5,  A13,  A16,  83,  247,  304 

WILLIAMS,  J. 

83 

WINICK,  J. 

A27,  A35,  431,  601 

WINTERSTEINER,  P.P. 

A27,  A35,  431,  601 

WORDEN,  H. 

A39 

WUEBBLES,  D.J. 

A45,  776 

YANG,  Y. 

A44,  764 

YEH,  H. 

All,  202 

YODMANS,  D. 

276 

YOSHIHO,  K. 

A20,  346 

ZARDECKI,  A. 

A13,  A16,  247,  304 

ZHOD,  D.K. 

A37,  A44,  629,  764 

17th  Annual  Review  Conference  on 
Atmospheric  Radiation  Models 


ABSTRACTS 


7-8  June  1994 
PL/Geophysics  Directorate 
Science  Center 


822 


FASCODE:  AN  UPDATE  AND  RECENT  VALIDATIONS/APPLICATIONS 
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Phillips  Laboratory 

Geophysics  Directorate 

29  Randolph  Road 

Hanscom  AFB,  MA  01731-3010 

*NRC  Fellow  (on  leave  from  Univ.  of  Michigan) 


FASCODE,  the  DoD  standard  line-by-line  (LBL)  code,  has  undergone  only  minor  changes 
over  the  past  two  years,  since  the  release  of  FASCOD3P  in  March  1992.  Present/future 
capabilities  and  their  driving  requirements  will  be  reviewed.  In  addition,  the  results  of  two  recent 
comparisons  of  FASCODE  against  other  LBL  codes,  both  initiated  by  the  WMO  (World 
Meteorological  Organization)  will  be  presented:  Intercomparison  of  Transmittance  and  Radiance 
Algorithms  (ITRA)  in  the  microwave  region,  and  InterComparison  of  Radiation  Codes  used  in 
Climate  Models  (ICRCCM)  in  the  more  historic  infrared  spectral-  range. 


M0DTRAN3,  M0DTRAN4,  AND  BEYOND 


A.  BerkS  P.K.  AcharyaS  L.S.  BernsteinS  D.C.  RobertsonS  G.P  Anderson^  J.H.  Chetwynd^ 
F.X.  Kneizys^y  L.A.  Hall^,  L.M.  Kimball^,  E.P.  Shettle^,  L.W.  Abreu^,  and  K.  Minschwaner^ 


‘Spectral  Sciences,  Inc. 
99  S.  Bedford  St.  #7 
Burlington,  MA 
01803-5169 


^Phillips  Laboratory 
29  Randolph  Road 
Hanscom  AFB,  MA 
01731-3010 


^Naval  Res.  Lab 
Code  7227 
Washington,  DC 
20375-5351 


.  "‘Ontar  Corp. 

9  Village  Way 
No.  Andover,  MA 
01845 


^NCAR 
Box  7000 
Boulder,  CO 


The  MODTRAN  code  continues  to  be  an  invaluable  tool  for  predicting  and  modeling 
radiation  transport  in  the  lower  atmosphere  (<  60km).  The  newest  version,  MODTRAN3,  will 
be  released  shortly.  Its  enhancements  include:  (1)  arbitrary  layer  resolution  (2)  improved 
treatments  for  single  scatter  solar  and  multiple  scatter  radiances  (3)  CFC,  UV  NOj,  UV  SO2  and 
new  O3  Chappius  cross-sections  (4)  a  spectrally  varying  surface  albedo  input  option  (5)  HITRAN 
94  based  band  model  data.  Current  research  which  wUl  lead  to  upgrades  for  MODTRAN4  and 
later  include:  improved  cooling  rate,  vertical  flux,  and  weighting  function  calculations;  higher 
spectral  resolution;  an  updated  aerosol  database;  vector  optimization;  and  incorporation  of  the 
NRL  climatology  database  via  the  SAG  code. 
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MOSART:  AN  UPDATE  ON  THE  MODERATE  SPECTRAL  ATMOSPHERIC 
RADIANCE  AND  TRANSMITTANCE  CODE 


William  M.  Cornette  Prabhat  K.  Acharya 

David  Robertson 


Gail  P.  Anderson 


Photon  Research  Assoc.  Inc. 
10350  N.  Toney  Pines  Road 
Suite  300 

La  Jolla,  CA  92037-1020 


Spectral  Sciences,  Inc. 
99  South  Bedford  St.#7 
Burlington,  MA  01803 


Phillips  Laboratory 
Geophysics  Directorate 
29  Randolph  Road 
Hanscom  AFB,  MA 
01731-3010 


The  MO^derate  ^ectral  Atmospheric  Radiance  and  Transmittance  Code  is  currently 
imdergoing  beta  testing  and  is  scheduled  for  release  in  August  1994.  Since  the  overview  last 
year,  a  number  of  new  capabilities  have  been  added.  These  include  an  earthshine  and  skyshine 
model,  improved  classification  for  tenain  types,  and  upgrades  to  the  global  terrain  database  and 
aerosol  model.  General  verification  and  validation  has  been  performed  primarily  with 
comparisons  to  MODTRAN,  while  limited  direct  validation  has  been  performed  with  data 
comparisons. 


PLEXUS 


PHILLIPS  LABORATORY  EXPERT  UNIFIED  SOFTWARE 


Dr.  Frank  O.  Clark 


Phillips  Laboratory 
Geophysics  Directorate 
29  Randolph  Road 
Hanscom  AFB,  MA  01731-3010 


PLEXUS,  or  Phillips  Laboratory  Expert  Unified  Simulator,  places  PL  Codes  in  a  unified 
knowledge  based  interface.  PLEXUS  is  intended  for  the  intelligent  non-expert  as  well  as  expert 
users.  PLEXUS  contains;  the  standard  Atmospheric  Generator  (SAG),  MODTRAN  (& 
LOWTRAN),  a  validated  integrated  MODTRAN  and  SHARC  pre-  and  post-  processor 
combination,  CBSD  (Celestial  Background  Scene  Descriptor).  With  a  single  set  of  input 
parameters  from  the  user,  the  preprocessor  automatically  nms  MODTRAN  below  50  km  and 
SHARC  above  50  km  and  produces  a  single  output,  permitting. sight  paths  which  cross  the 
equilibrium/non-equilibrium  boundary.  PLEXUS  also  contains  the  new  Standard  Atmospheric 
Generator  which  is  based  on  NASA  MSIS-90  and  NRL  lower  atmosphere  phenomenology  data. 
This  seamless  model  implementation  is  essential  for  modem  space  based,  TMD,  and  GPALS 
applications.  The  single  geometry  interface  applies  to  all  codes. 


EOSAEL92  UPDATE 


Alan  E.  Wetmore 
U.S.  Army  Research  Laboratory 
White  Sands,  Missile  Range,  NM  88002-5501,  U.S.A. 


Since  the  distribution  of  EOSAEL92  several  of  the  modules  have  been  updated  and  an  acous¬ 
tics  module  has  been  selected  for  incorporation.  The  TARGAC  model  has  been  extended  and  an 
interface  to  run  the  module  under  Microsoft  Windows  has  been  developed.  The  UVTRAN  model 
also  has  a  version  that  runs  under  Microsoft  Windows.  Both  of  these  will  be  demonstrated  at  the 
conference.  The  millimeter  wave  module,  NMMW,  has  been  updated  with  a  new  gaseous  absorp¬ 
tion  algorithm.  The  Scanning  Fast  Field  Program  (SCAFFIP)  is  an  acoustics  model  that  predicts 
acoustic  propagation  based  on  spherical  spreading,  molecular  absorption,  refraction,  acoustically 
complex  ground  impedance,  and  diffraction  over  being  terrain.  We  are  in  the  process  of  making  all 
of  these  modules  and  documentation  available  through  the  Test  and  Evaluation  Community  NETwork 
(TECNET).  In  addition  the  Validation  and  Verification  of  the  laser  transmission  module,  LZTRAN, 
has  been  performed. 


THE  NAVY  OCEANIC  VERTICAL  AEROSOL  MODEL 
IN  A  FORTRAN  SUBROUTINE  FORMAT 


Stuart  G.  Gathman 


NCCOSC  RDTE  DIV  543 
53170  Woodward  Road 
San  Diego,  CA  92152 


A  FORTRAN  version  of  the  Navy  Oceanic  Vertical  Aerosol  Model,  NOV  AM,  has  been 
completed  and  is  ready  to  be  interfaced  with  other  types  of  transmission  models.  The  model  is 
available  in  the  form  of  a  NCCOSC  RDT&E  DIV  technical  report  #1634  which  contains:  (a)  a 
written  description  of  the  model,  (2)  the  printed  code  of  the  subroutine,  (3)  the  printed  code  of 
the  associated  subroutines  and  functions,  and  (4)  a  sample  driver  which  can  be  used  for  batch 
processing  of  data.  The  model  requires  meteorological  data  in  the  form  of  both  surface  and 
atmospheric  sounding  data  files.  NOVAM  predicts  the  vertical  distribution  of  aerosol  in  the  first 
6  km  above  the  ocean.  Outputs  of  the  model  include  optical  properties  in  the  wavelength  band 
from  0.2  to  40  micrometers  at  any  altitude.  These  properties  include  the  volume  extinction  and 
absorption  coefficient,  the  relative  humidity,  as  well  as  parameters  describing  the  aerosol  size 
distribution  in  terms  of  dN/dr  at  that  altitude.  This  paper  will  describe  the  model  and  some 
sample  applications. 


THE  NATURE  OF  DESERT  AEROSOLS:  SUMMARY  OF  CHINA  LAKE 

STUDIES  TO  DATE 


P.L.  Walker 

Physics  Department,  Naval  Postgraduate  School,  Monterey,  CA  93943-5117 

and 

L.A.  Mathews  retired  from 
Naval  Air  Warfare  Center,  China  Lake,  CA 


The  Naval  Air  Warfare  Center  has  been  conducting  desert  aerosol  characterization  measure¬ 
ments  at  China  Lake  at  ground  level  in  1987;  at  15,000  feet  in  connection  with  Long  Jump  in  1988; 
at  ground  level  in  the  western  Mojave  including  Edwards  AFB  and  China  Lake  in  1990;  and  will  be 
conducting  further  measurements  in  the  western  Mojave  in  the  summer  of  1994.  This  presentation 
will  present  an  over  view  of  the  analysis  of  data  to  date.  The  properties  of  desert  aerosol  that  are 
appearing  are  somewhat  different  than  those  proposed  by  Longtin  and  Shettle  in  that  the  accumula¬ 
tion  mode  is  dominated  by  the  presence  of  large  numbers  of  organic  carbon  particles  while  dust  is 
characterized  by  the  presence  of  clay  particles  in  the  dust  size  distribution  mode  with  few  quartz 
particles  being  present  and  no  quartz  dominated  blowing  sand  mode.  We  will  compare  particle  size 
and  composition  at  15,000  feet  and  ground  level,  further  micrographic  studies  of  the  composition 
of  high  altitude  aerosols  and  the  effect  of  wind  speed  on  aerosol  characteristics  at  ground  level. 


APPLICATION  OF  THE  MODTRAN2  CODE  TO  THE  MODELING  OF  SILICATE 

DUST  CLOUDS. 


S.  Mazuk  and  D.K.  Lynch 
The  Aerospace  Corporation 
P.O.Box  92957 
Mail  Code  M2/255 
Los  Angeles,  CA  90009-2957 


The  MODTRAN2  code  has  been  used  to  model  the  thermal  infrared  radiance  from  a  dust 
cloud  formed  during  tin  explosives  test  at  the  White  Sands  Missile  Range  in  1993.  High  signal-to- 
noise  infrared  spectra  of  the  dust  cloud  were  measured,  indicating  the  presence  of  silicate  particles. 
A  silicate  dust  cloud  was  defined  as  part  of  the  model  atmosphere,  and  MODTRAN  was  used  to 
calculate  radiance  spectra  between  6  and  16  microns.  A  description  of  the  analysis,  comparison 
with  experimental  data,  and  the  problems  encountered  will  be  presented. 
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MODELIZATION  OF  CLOUD  RADIATION  AND  COMPARISON  WITH 
AIRBORNE  MEASUREMENTS 


C.  Malherbe,  P.  Simoneau  and  A.  Boischot 

Office  National  d’Etudes  et  de  Recherches  Aerospatiales 

BP  72,  92322  Chatillon  Cedex 

France 


We  have  developed  a  computer  code  named  NUALUM,  dedicated  to  the  computadon  of 
the  atmospheric  radiance  in  presence  of  clouds. 

This  radiance  is  the  sum  of  three  terms: 

-  the  atmospheric  radiance  between  the  cloud  and  the  observator 

-  the  thermal  radiance  of  the  cloud 

-  the  scattering  by  the  cloud  of  the  solar,  the  terrestrial  and  the  atmospheric  fluxes. 

The  bidirectional  reflectivity  and  the  emissivity  of  the  cloud  are  computed  with  the 
DISORT  program.  The  solar  flux,  the  terrestrial  flux,  and  the  atmospheric  flux  are  calculated 
by  MODTRAN.  Numerical  results  are  compared  with  spectral  and  spatial  infiraied  measurements 
performed  by  airborne  devices. 
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OSIC  -  AN  ULTRAVIOLET  TRANSMISSION  AND  MULTIPLE  SCATTERING  MODEL 


Michael  Neer  and  Kathy  Crow 


SciTec,  Inc. 

100  Wall  Street 
Princeton,  NJ  08540 


A  novel  analytical  approach  is  presented  for  predicting  not  only  conventional  transmission 
associated  with  molecular  and  aerosol  scatter  and  absorption,  but  also  the  contribution  of  multiple 
scattering  including  angle  of  arrival.  The  approach  has  been  implemented  in  a  model  called 
MWR  and  sample  signature  versus  time  predictions  will  be  presented  along  with  comparisons 
to  experimental  data.  The  program  will  run  much  faster  than  Monte  Carlo  models  so  it  is  better 
suited  to  carrying  out  missile  signature  parametric  studies.  The  model  can  also  be  used  to  predict 
signatures  when  the  source  is  not  in  the  sensor  field  of  view  and  the  received  signal  comes  only 
from  scattered  radiation. 
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ATMOSPHERE  UV-VISIBLE  RADIATION  FIELD  PERTURBATIONS  PRODUCED  BY 
CLOUDS  AND  AEROSOLS:  COMPARISON  OF  RESULTS  USING  MONTE  CARLO, 
DISORT  AND  INTEGRAL  EQUATION  METHODS 


D.E.  Anderson  and  R.  DeMajistre 


The  Johns  Hopkins  University 
Applied  Physics  Laboratory 
Laurel,  MD  20723 


The  UV- Visible  stratosphere  radiation  can  be  strongly  perturbed  from  the  clear 
atmosphere  case  when  clouds  and/or  aerosols  are  present.  The  clear  atmosphere  is,  to  within 
-^2%  accuracy,  isotropic.  When  clouds  and/or  aerosols  are  present  the  radiation  field  can  become 
anisotropic  to  a  degree  that  isotropy  is  no  longer  a  useful  assumption,  an  inclusion  of  these 
effects  is  necessary  either  by  approximation  or  by  a  more  representative  numerical  description. 
We  compare  the  results  from  three  different  approaches  to  the  calculation  of  the  radiation  field 
when  clouds  and/or  aerosols  are  present.  Results  are  presented  for  the  280-800nm  wavelength 
region  and  for  a  range  of  solar  zenith  angles. 

We  find  that  the  integral  equation  solution  assuming  isotropic  scattering  is  accurate  to  for 
a  clear  atmosphere  and  when  similarity  transformations  are  included,  accurate  to-~5-20%  when 
aerosols  and/or  clouds  are  present,  the  variation  in  accuracy  depending  on  wavelength,  solar 
zenith  angle,  altitude  and  cloud  ot  aerosol  extinction  optical  depth.  The  Monte  Carlo  and 
DISORT  results  employ  either  the  actual  phase  function  calculated  from  a  Mie  scattering  code, 
or  a  Henyey-Gieenstein  representation  of  the  phase  function. 
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INCLUSION  OF  ACCURATE  MULTIPLE  SCATTERING  IN  MODTRAN 


K.  Stamnes,  S.  Tsay,  N.  Larsen  M.  Yeh 


Geophysical  Institute  Caelum  Research  Corp. 

University  of  Alaska  1 1229  Lockwood  Drive 

Fairbanks,  AK  99775-7320  Silver  spring,  MD  20901 


An  accurate  multiple  scattering  scheme  has  been  incorporated  into  MODTRAN  by 
merging  it  with  the  general-purpose  multiple  scattering  algorithm  DISORT.  The  multiple 
scattering  source  functions  for  solar  as  well  as  thermal  radiation  are  computed  in  DISORT  and 
fed  into  MODTRAN  where  they  are  combined  with  the  single  scattering  solar  and  thprmai  source 
terms.  The  complete  intensity  distribution  is  then  computed  in  MODTRAN  by  integrating  along 
the  line-of-sight.  The  importance  of  this  improved  multiple  scattering  computation  is 
demonstrated  by  comparing  results  from  the  original  MODTRAN  with  those  produced  by  the 
merged  MODIRAN/DISORT  package  for  radiation  in  the  solar,  thermal,  and  solar/thermal 
crossover  spectral  regions.  In  the  ultraviolet/visible  region  the  improvements  are  significant  even 
for  clear  sky  conditions,  while  multiple  scattering  effects  are  important  in  the  thermal  region  only 
in  the  presence  of  clouds  and  aerosols.  Comparisons  with  experimental  data  wiU  also  be  used 
to  demonstrate  the  superior  performance  of  the  merged  package. 
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AN  APPLICATION  OF  RADIATIVE  TRANSFER  THEORY 
TO  UNDERSTANDING  AEROSOL  MTF 


David  H,  Tofsted 
Alan  E.  Wetmore 
Richard  C.  Shirkey 
U.S.  Army  Research  Laboratory 
White  Sands  Missile  Range,  NM  88002-5501 

Andrew  2^rdecki 
Los  Alamos  National  Laboratory 
Los  Alamos,  NM  87544 

Brian  Davis 

Physical  Sciences  Laboratory 
Las  Cruces,  NM  88003 


Aerosol  forward  scattering  effects  can  significantly  blur  distant  objects  when  the  Mie  size 
parameter  becomes  large  (such  as  snow,  rain,  fog  in  the  visible).  Starting  with  previous  radiative 
transfer  forward  scattering  approximations,  we  have  developed  a  new  means  for  determining  these 
effects  on  the  aerosol  MTF.  This  derivation  allows  for  inclusion  of  horizontal  inhomogeneities  of 
aerosol  densities  and  species,  uses  the  aerosol  forward  scattering  approximation  of  radiative  trans¬ 
fer  theory,  and  models  the  aerosol  phase  function  forward  peak  as  a  sum  of  Gaussian  components. 
Using  Gaussian  phase  function  components  yields  an  exact  solution  to  the  radiative  transfer  func¬ 
tion  under  the  forward  scattering  approximation.  The  results  of  the  radiative  transfer  calculation 
have  been  translated  directly  into  terms  compatible  with  a  linear  filter  interpretation  of  incoherent 
imaging.  Thus  the  aerosol  MTF  derived  directly  from  the  radiative  transfer  analysis  and  does  not 
depend  on  system  characteristics. 


MOLECULAR  AND  AEROSOL  EFFECTS  ON  AIRBORNE  LASER  PROPAGATION 


y 


Larrene  Harada  and  Daniel  Leslie 


WJ.  Schafer  Associates,  Inc. 

1901  North  Fort  Myer  Drive,  Suite  800 
Arlington,  VA  22209 


Recent  missions  for  lasers  on  airborne  platforms  require  propagation  through  long, 
horizontal,  high  altitude  slant  paths.  Propagation  is  affected  by  upper  tropospheric  and 
stratospheric  molecules  and  aerosols.  Extinction  coefficients  and  tranmissions  for  specific  high 
energy  laser  and  laser  radar  wavelengths  have  been  calculated  using  FASCOD3P/HITRAN92. 
The  relative  importance  of  various  molecular  constituents  are  affected  because  the  influence  of 
water  vapor  is  significantly  less  than  it  is  near  the  surface.  Stratospheric  aerosols  can  cause 
considerable  attenuation  following  a  major  volcanic  eruption.  We  have  used  NASA  SAGE 
(Stratospheric  Aerosol  and  Gas  Experiment)  satellite  data  to  evaluate  the  effects  of  the  Mount 
Pinatubo  eruption  on  airborne  laser  propagation  at  certain  wavelengths.  We  have  also  determined 
probability  of  occurrences  of  cirrus  at  the  air-to-air  mission  path  altitudes  for  different  seasons. 
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GROUND-TO-SPACE  ATMOSPHERIC  TRANSMITTANCE  MEASUREMENTS  IN  THE 
3-5  AND  8-12)iin  SPECTRAL  REGIONS:  COMPARISON  WITH  LOWTRAN  7 


Adam  D.  Devir,  N.  Brandman,  B.  Barzilai  and  A.  Ben  Shalom 
Electro  Optics  R&  Division  (EORD) 

Technion  R&D  Foundation,  Haifa  32000,  ISRAEL 


The  computer  code  LOWTRAN  is  widely  and  extensively  used  for  the  prediction  of  the  propa¬ 
gation  of  IR  radiation  through  the  atmosphere.  The  latest  version  of  this  code,  the  LOWTRAN  7,  is 
assumed  to  be  the  most  elaborate  and  accurate  one.  Hence  it  was  decided  to  test  its  validity  by 
comparing  its  prediction  for  ground-to-space  slant  paths  with  actual  transmittance  measurements 
with  the  sun  as  a  blackbody  source.  A  good  agreement  between  the  theoretical  predictions  and  the 
experimental  results  was  obtained  in  the  8-12jx  m  spectral  region  for  all  zenith  angles  between  90° 
and  60°  (0°  to  30°  above  the  horizon).  In  the  3  -5|xm  spectral  region  some  spectral  discrepancy  was 
observed  though  the  value  of  the  integrated  measured  transmittance  agreed  well  with  the  predicted 
one.  The  reasons  for  this  can  be  assumed  to  be  in  the  new  band  model  absorption  parameters  (in  the 
3-4pm  region)  and  in  the  water  continuum  model  (in  the  4.4-5.2|im  region). 


WEATHER  AND  ATMOSPHERIC  VISUALIZATION  EFFECTS  (WAVES) 

FOR  SIMULATION 


X 


Richard  Shirkey 
David  Tofsted 
Alan  Wetmore 
Army  Research  Laboratory 
Battlefield  Environment  Directorate 
White  Sands  Missile  Range,  NM 


Andrew  Zardecki 
Los  Alamos  Consulting 
Los  Alamos,  NM 


The  WAVES  suite  of  codes  ascertains  and  visualizes  environmental  effects  due  to  natural 
clouds/haze  and  fog.  WAVES  determines  the  illumination  passing  through  multiple  inhomogeneous 
cloud  layers  and  the  resulting  radiance  field.  The  effects  of  aerosol  forward  scattering  and  optical 
turbulence  are  also  included. 


<T 
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SPARTA’S  LIDAR  SIMULATION  CODE,  BACKSCAT  VERSION  4.0 


D.R.  Longtin,  M.G.  Cheifetz,  J.R.  Jones,  J.R.  Hummel 


SPARTA,  Inc. 

24  Hartwell  Avenue 
Lexington,  MA  02173 


Since  1988,  SPARTA,  Inc.  has  been  developing  a  lidar  simulation  package  for  the  IBM  PC 
environment  called  BACKSCAT.  The  package  simulates  the  backscatter  return  from  lidar  systems 
of  different  designs,  viewing  geometries,  and  atmospheric  conditions.  This  paper  will  describe  the 
latest  version,  BACKSCAT  Version  4.0,  and  present  sample  results.  The  package  will  be  available 
for  distribution  within  a  few  months.  BACKSCAT  Version  4.0  includes  many  new  technical  fea¬ 
tures.  Foremost  is  the  ability  to  evaluate  signal-to-noise  performance  for  direct  detection  and  co¬ 
herent  Doppler  lidar  systems.  Signal-to-noise  ratios  are  then  used  to  estimate  the  accuracy  of  range 
and  wind  speed  measurements  for  the  chosen  lidar  configuration.  These  capabilities  transform 
BACKSCAT  into  a  more  general  purpose  lidar  simulation  tool  for  users  who  want  to  design  a  lidar 
system  to  explore  specific  aspects  of  the  atmosphere.  In  response  to  the  user  community,  the 
second  new  feature  is  the  ability  to  include  water  clouds  in  lidar  simulations.  Although  the  existing 
user-defined  aerosol  option  could  be  used  to  represent  water  clouds,  it  required  detailed  cloud 
particle  size  distribution  data  that  might  not  be  available  to  some  users.  Thus,  as  an  alternate  ap¬ 
proach,  water  clouds  now  can  be  “clicked  on”  automatically  as  a  built-in  cloud  model.  The  third 
new  feature  is  auxiliary  software  that  estimates  the  molecular  absorption  profile  for  a  particular 
lidar  wavelength.  Although  previous  versions  of  BACKSCAT  provided  the  means  to  include  mol¬ 
ecule  absorption  in  a  simulation,  the  code  did  not  provide  a  tool  to  calculate  it.  The  molecular 
absorption  package  has  been  nicknamed  MABS  and  generates  output  that  can  be  directly  used  in  a 
lidar  simulation.  It  must  be  noted  that  the  spectral  resolution  in  MABS  probably  exceeds  that  of 
most  lidars,  so  the  package  is  not  intended  to  be  a  complete  treatment  of  the  molecular  absorption 
problem. 


*  Work  Performed  Under  Contract  F19628-C-9 1-0093 


THE  SOLAR  IRRADIANCE  BY  COMPUTATION 


Robert  L.  Kurucz 


Harvard-Smithsonian  Center 
for  Astrophysics 
60  Garden  Street 
Cambridge,  MA  02138 


I  am  now  able  to  compute  a  purely  theoretical  model  photosphere  that  reproduces  the 
irradiance  measurements  of  Neckel  and  Labs  in  the  visible  for  bandpasses  of  approximately  2  nm. 
That  model  and  Avrett’s  empirical  quiet  sun  model  are  used  to  predict  the  irradiance  out  to  200 
microns.  At  200  microns  the  chromosphere  in  the  empirical  model  is  responsible  for  most  of  the 
flux.  Actually,  monochromatic  irradiance  is  required  to  determine  the  absorption  by  molecular 
lines  in  the  atmosphere.  I  am  producing  atlases  of  the  solar  flux,  central  intensity,  and  limb 
spectra  taken  by  James  Brault  at  Kitt  Peak.  One  atlas  "Solar  Flux  Atlas  from  294  to  1300  nm" 
by  Kurucz,  Furenlid,  Brault,  and  Testerman  (1984),  has  been  published  thus  far.  I  also  have  the 
"Photometric  Atlas  of  the  Solar  spectrum  from  1,850  to  10,(X)0  cm-l"  by  Delbouille,  Roland, 
Brault,  and  Testerman  (1991)  also  taken  at  Kitt  Peak.  In  addition  I  have  the  ATMOS  central 
intensity  atlas  from  650  to  48(X)  cm-l  taken  by  Farmer  and  Norton  (1989)  from  Spacelab  3. 
Parts  of  the  flux  atlases  directly  give  the  irradiance  spectrum  but  much  is  confused  or  obscmed 
by  terrestrial  lines.  The  intensity  atlases  can  be  used  as  guides  for  deconvolving  blends  because 
the  solar  lines  are  narrower.  One  path  I  am  pursuing  is  fitting  the  transmission  spectram 
empirically  and  simply  dividing  it  out  However,  the  regions  where  the  transmission  is  too  low, 
there  is  no  alternative  but  to  compute  a  purely  theoretical  solar  spectrum.  I  use  these  atlases  to 
test  the  pure  calculations  of  solar  spectra  and  of  transmission  spectra.  I  identify  problems  with 
the  line  data  and  I  try  to  make  generic  corrections  that  improve  hundreds  or  thousands  of  lines 
at  a  time.  If  the  spectrum  calculations  look  good  in  the  regions  that  have  high  transmission  are 
computed  accurately.  The  main  problem  has  been  continuum  placement  Ozone  and  01  "dimer" 
features  are  difficult  to  determine  because  the  atlases  are  each  made  up  of  a  number  of  sharply 
peaked  FTS  scans.  The  continuum  placement  affects  the  appearance  of  line  wings  and  the 
apparent  depth  of  weak  features.  To  get  a  feel  for  the  scope  of  the  monochromatic  irradiance 
problem  I  have  computed  the  spectrum  from  150  nm  to  200  microns  at  a  resolution  of  500000 
using  58  million  lines,  both  predicted  and  observed.  These  are  the  same  lines  that  went  into 
computing  the  opacity  for  the  model  photosphere.  If  this  spectrum  is  degraded  to  the  resolution 
of  the  model  it  looks  like  the  model.  At  any  given  monochromatic  point  the  spectrum  is  not 
reliable.  But  at  a  resolution  of  10(X)0,  say,  it  approaches  measurement  accuracy.  In  regions  of 
low  transmission  it  is  more  reliable  than  existing  measurements.  I  will  publish  tables  of  these 
irradiance  spectra  at  various  resolutions. 


NEW  VISIBLE  AND  NEAR  IR  OZONE  ABSORPTION  CROSS-SECTIONS 

FOR  MODTRAN 


Eric  P.  Shettle,  Code  7227,  Remote  Sensing  Div.,  NRL,  Washington,  DC  20375 

and 

Stuart  M.  Anderson,  Dept,  of  Physics,  Augshurg  College,  Minneapolis,  MN  55454 


Knowledge  of  the  ozone  absorption  cross-sections  is  important  for  modeling  the  propagation 
of  radiation  through  the  atmosphere,  and  is  critical  for  deriving  the  atmospheric  ozone  concentra¬ 
tion  from  measurements  of  the  spectral  character  of  transmitted  or  scattered  radiation.  While  the 
ultraviolet  cross-sections  for  ozone  were  completely  replaced  in  LOWTRAN  7,  the  visible  cross- 
sections  are  unchanged  since  the  original  version  of  LOWTRAN  eind  the  cross-sections  are  not 
given  at  all  for  wavelengths  longer  than  769  nm.  Recently  there  have  been  a  several  new  spectro¬ 
scopic  studies  of  the  ozone  Chappuis  and  Wulf  absorption  bands,  along  with  measurements  of  the 
absolute  cross-section  at  selected  wavelengths.  We  have  normalized  the  spectral  measurements  to 
the  absolute  cross-section  data,  and  smoothly  joined  the  different  sets  of  measurements,  to  provide 
the  ozone  cross-section  from  407  to  1089  nm.  The  temperature  dependence  is  also  given  for  407  to 
762  nm.  It  is  anticipated  that  these  cross-sections  will  be  included  in  MODTRAN  3. 
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ABSORPTION  CROSS  SECTION  MEASUREMENTS  OF  CARBON  DIOXIDE  IN  THE 
WAVELENGTH  REGION  118.7-175.5  nm  AND  THE  TEMPERATURE  DEPENDENCE 


K.  Yoshino,  J.R.  Esmond  and  W.H.  Parkinson  K.  Ito  and  T.  Matsui 


Harvard-Smithsonian  Center 
for  Astrophysics 
60  Garden  Street 
Cambridge,  MA  02138 


Photon  Factory,  KEK 
Tsukuba,  IbaraJd  305 
Japan 


In  absorption,  CO2  is  transparent  down  to  at  least  210  nm.  Photoabsorption  cross  section 
measurements  of  COj  in  the  UV  to  VUV  region  have  been  made  by  several  investigators. 
Because  of  lack  of  resolution,  none  of  the  studies  in  the  period  1950-70  resolved  the  structures 
in  the  COj  absorption  cross  sections.  After  investigating  photographic  plates  of  CO2  absorption 
spectra  taken  in  high  resolution,  we  confirmed  that  structures  in  the  CO2  are  broad  enough  to 
allow  the  use  of  a  medium  resolution  instrument.  The  3-m  spectrometer  on  the  BL-20A  beam 
line  at  the  Photon  factory,  Japan  was  used  in  the  first  order  of  a  1200  1/mm  grating  to  obtain 
absorption  cross  sections  of  CO2.  The  background  continuum  is  provided  by  the  synchrotron 
radiation  source  without  any  order  separation,  because  any  high  order  radiation  is  cut  off  by 
MgFj  windows.  The  absorption  cross  sections  of  CO2  were  obtained  at  two  temperatures,  195 
K  and  295  K,  in  the  wavelength  region  120-175  nm  for  every  0.005  mn  step.  Almost  no 
temperature  dependency  was  observed  in  the  peak  cross  sections  in  the  wavelength  region  125- 
143  nm,  however,  the  ratio  of  peak  cross  sections  are  down  to  0.7  around  175  nm.  This  work 
is  supported  by  NASA  Upper  Atmospheric  Research  Program  under  the  Grant  No.  NAGS-494 
to  the  Smithsoiuan  Astrophysical  Observatory. 
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FASCODE  ENVIRONMENT  (EASE) 


J.L.  Moncet  and  W.O.  Gallery 
Atmospheric  and  Environmental  Research,  Inc. 

840  Memorial  Drive,  Cambridge,  MA  02139 

and 

G.  P.  Anderson 

Phillips  Laboratory,  29  Randolph  Road,  Hanscom  AFB,  MA  01731-3010 


The  Optical  Physics  Division  of  the  Phillips  Lab  is  developing  a  state-of-the-art  line-by-line 
atmospheric  transmittance/radiance  model  as  the  successor  to  FASCODE.  Under  a  joint  program 
with  DoE’s  Atmospheric  Radiation  Measurement  (ARM)  project,  the  best  features  of  FASCODE 
and  the  DoE’s  standard  radiative  transfer  program  LBLRTM  will  be  merged  and  new  capabilities 
developed.  While  retaining  all  of  FASCODE’s  capabilities,  EASE  will  specifically  address  envi¬ 
ronmental  applications  such  as  atmospheric  heating  rates  and  temperature  and  species  inversion. 

The  main  features  of  EASE  will  be: 

-  Highly  modular  all  the  basic  function-l-b-1  optical  depth,  radiative  transfer,  atmospheric 
profile,  filtering,  scanning  and  plotting- will  be  broken  out  into  independent  programs 
These  program  can  be  run  together  through  a  driver  just  as  FASCODE  is  now,  run 
separately,  or  incorporated  into  customized  applications:  e.g.,  multiple  scattering  codes 
or  inversion  routines 

-  Core  module:  line-by-line  calculation  of  optical  depth:  incorporate  best  Physics: 
including  continuua,  line  coupling,  cross-sections  of  heavy  molecules. 

-  Optimized:  for  vectorized  and/or  parallel  processing 

-  Interface:  improved  user-friendliness 

-  Improved  packaging:  makefile  for  easy  installation  on  a  wide  variety  of  platforms. 

-  Configuration  controlled:  for  ease  of  maintenance 


SSGM:  THE  STRATEGIC  SCENE  GENERATION  MODEL 


S.  McKenzie  and  R.  Armstrong 


Mission  Research  Coip. 
One  Tara  Blvd.,  Suite  302 
Nashua,  NH  03062 


The  Strategic  Scene  Generation  Model  generates  viewer-perspective  pixel  apparent 
radiance  maps  and  point-source  intensity  data.  These  scenes  are  computed  from  the  output  of 
state-of-the-science  phenomenology  codes  and  their  authenticated  input  databases.  Scenes  consist 
of  quiescent  and  enhanced  natural  and  perturbed  backgrounds  with  embedded  targets  and  target 
related  events.  Backgrounds  include  terrain,  clouds,  earthlimb,  space,  aurora,  horizon,  and 
nuclear  events.  Targets  include  boosting  missile  plumes  and  hardbodies,  fuel  vents,  and  post¬ 
boost  objects  including  RV  and  decoys.  A  review  of  the  robust  Earthlimb  Capabilities  in  Release 
6.0  will  be  presented. 


LIDAR  MEASUREMENTS  OF  ATMOSPHERIC  OPTICAL  PROPERTIES 


p 


C.R.  Philbrick,  T.D.  Stevens,  and  S.  Manivada 


Penn  State  University 
Dept  of  Electrical  Engineering  and 
Applied  Research  Laboratory 
University  Park,  PA  16802 


The  LAMP  (Lidar  Atmospheric  Measurements  Program)  lidar  instrument  has  been  used 
to  develop  and  refine  techniques  for  measurements  of  atmospheric  properties.  The  initial 
measurements  primarily  focused  on  the  measurement  of  middle  atmsophere  dynamics  and 
structure.  More  recent  efforts  have  concentrated  on  measurements  of  tropospheric  water  vapor 
and  temperature.  The  vibrational  and  rotational  Raman  scatter  signals  used  for  these  lower 
atmospheric  profiles  can  also  be  used  to  determine  the  backscatter  and  extinction  profiles.  The 
optical  extinction  profiles  determined  from  the  Raman  scatter  signals  at  several  wavelengths 
provide  a  measure  of  the  wavelength  dependence  of  aerosol  and  cloud  layers  in  the  troposphere. 
Several  examples  of  the  optical  extinction  profiles  have  been  examined  and  comparisons  with 
inversion  techniques  from  the  backscatter  have  shown  limitations  of  the  inversion  analysis.  The 
inversion  of  backscatter  profiles  does  not  generally  yield  satisfactory  results.  These  investigations 
show  that  the  Raman  lidar  measurements  provide  a  most  useful  way  to  measure  the  atmsoheric 
optical  properties. 
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GROUND-BASED  MEASUREMENTS  OF  HE  AND  HCI 


Hilary  E.  Snell  and  Paul  B.  Hays 
University  of  Michigan 


This  poster  describes  remote  sensing  observations  of  HF  and  HCI  column  densities  over  Ann 
Arbor,  Michigan.  Spectra  were  collected  between  October  1992  and  July  1993  using  a  ground- 
based  Michelson  Interferometer  and  the  technique  of  solar  absorption  spectroscopy.  While  the 
column  densities  computed  from  these  spectra  agree  with  other  groups'  measurements  of  these 
species,  we  have  noticed  that  the  HCI  spectral  line  at  2925.897  cm'*  appears  strongly  asymmetric. 
It  is  our  belief  that  this  is  due  to  spectral  contamination  by  another  chemical  species  and  not  instru¬ 
mental  phase  errors  as  previously  reported  by  other  groups. 


AN  ATMOSPHERIC  MODEL  FOR  GRAVITY  WAVE  INDUCED  TURBULENT  LAYERS 
(BLINI)  BASES  ON  THE  SATURATED  CASCADE  MODEL 


E.M.  Dewan 

N.  Grossbard 

T.  VanZandt 

Geophysics  Directorate 
Phillips  Laboratory 

29  Randolph  road 
Hanscom  AFB,  MA 

Boston  College 

NOAA 

01731-3010 

Turbulence  in  the  atmosphere  above  the  boundary  layer  occurs  in  thin,  pancake-shaped 
horizontal  layer  called  "blini".  These  are  generated  by  gravity  waves  whenever  the  Richardson 
number  is  less  than  a  critical  value.  A  theory  which  seems  to  accurately  predict  atmospheric 
gravity  wave  spectra  was  recently  published  in  GRL  (Dewan,  May  1, 1994,  vol.  21,  p.817).  This 
theory  also  implies  wave  number-period  relationships  used  below.  On  the  basis  of  this  theory, 
simulations  were  made  of  gravity  wave  fields  by  means  of  a  superposition  of  random  (phase  and 
amplitude)  sine  waves  consistent  with  the  vertical  wave  number  spectrum  of  horizontal  wind 
fluctuations.  Whenever  shears  exceeded  critical  values,  the  computer  program  graphed  the 
locations  on  an  altitude-plot  As  a  result,  realistic  blini  appear  which,  as  was  anticipated,  "blink" 
on  and  off.  Subsequently  this  model  will  go  by  the  name  "The  blinking  blini  model  for 
atmospheric  turbulence".  It  will  be  used,  in  the  future,  together  with  laser  propagation  codes  to 
ascertain  the  effects  of  wave  induced  turbulence  on  ABL  systems. 


PRELIMINARY  RESULTS  FROM  A  RECENT  ICRCCM  INITIATIVE 


F.  Miskolczi,  M.  Bonzagni,  R.  Guzzi 
University  of  Maryland 


In  this  paper  high  resolution  measured  and  computed  atmospheric  downward  radiance  spectra 
are  compared.  The  spectral  measurements  were  taken  during  the  SPECTRE  experiment  (Coffeywille, 
Kansas  1991)  using  the  UWI AERI  fourier  spectrometer.  The  radiance  modeling  is  based  on  the 
High-Resolution  Atmospheric  Radiance-Transmittance  Code  (HARTCODE)  (F.  Miskolczi,  M. 
Bonzagni  andR.  Guzzi,  1990).  In  this  comparison  the  seven  major  absorbers,  Fll,  F12,  N2  and 
CCL4  were  included. 


Reference: 

F.  MISKOLCZI,  M.  BONZAGNI,  R.  GUZZI  (1990) 

High  Resolution  Atmospheric  Radiance-Transmittance  Code  (HARTCODE). 
Meteorology  and  Environmental  Sciences,  World  Scientific 
Pub.  Co.,  1990,  pp.  743-790 


STRUCTURE  IN  RADIATIVE  EXCITATION  AS  A  SOURCE  OF  HIGH  ALTITUDE 
RADIANCE  STRUCTURE;  CO(v=l)  RADIANCE 


J.R.  WinickS  R.H.  Picard*,  P.P.  Wintersteiner^,  and  J.A.  Dodd’ 


‘Phillips  Laboratory 
29  Randolph  Road 
Hanscom  AFB,  MA 
01731-3010 


^Arcon  Corp.  ’Stewart  Radiance  Laboratory 

Waltham,  MA  129  Great  Road 

Bedford,  MA  01730 


We  discuss  a  largely  overlooked  source  of  upper  atmosphere  stractured  radiance,  stracture 
in  the  earthshine  pumping  rate  caused  by  a  structured  lower  boundary.  The  4.7  pm  radiance  from 
CO(v=l)  is  a  good  example  of  a  non-LTE  emission  that  is  dominated  by  radiative  processes. 
The  analysis  of  limb  spectra  of  both  the  ^^C‘®0  and  ^’C‘®0  from  CIRRIS  lA  data  (Dodd  et  al., 
GRL,  ^  2683,  1993)  clearly  indicated  how  the  differing  opacity  of  the  isotopes  controlled  the 
effect  of  pumping  from  the  lower  boundary.  We  examine  some  limiting  cases  that  show  that 
nighttime  earthlimb  radiance  variations  of  factors  greater  than  2  can  be  caused  by  cold  clouds 
at  the  tropical  tropopause. 
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INFRARED  RADIANCE  FLUCTUATIONS  IN  THE  UPPER  ATMOSPHERE 


J.  GruningerS  R.L.  Sundberg^,  P.  De*,  and  J.  Brown^ 


‘Spectral  Sciences,  Inc.  ^Phillips  Laboratory 

99  S.  Bedford  St.  #7  29  Randolph  Road 

BurUngton,  MA  01803  Hanscom  AFB,  MA  01731-3010 


Simulation  of  infrared  radiance  fluctuations  in  the  atmosphere  depends  on  detailed 
descriptions  of  fluctuations  in  atmospheric  species  number  densities,  vibrational  state  populations, 
and  the  kinetic  temperatures  along  the  sensor  Une-of-sight.  In  the  upper  atmosphere,  where 
nonequilibrium  effects  (LTNE)  become  important,  a  small  fluctuation  in  the  kinetic  temperature 
can  result  in  correlated,  anti-correlated,  or  no  change  in  the  vibrational  temperature  of  states.  The 
relationship  between  kinetic  and  vibrational  temperature  fluctuations  depends  on  the  subtle 
interplay  between  changes  in  the  total  number  densities,  changes  in  the  temperature-dependent 
kinetic  rates,  and  the  relative  contribution  of  the  radiative  relaxation.  The  model  developed  in 
this  paper  predicts  the  two-dimensional  radiance  covariance  function  for  LTNE  conditions.  The 
radiance  statistics  are  non-stationary  and  are  explicitly  bandpass  and  sensor  FOV  dependent.  The 
SHARC  model  is  used  to  calculate  mean  LOS  radiance  values  and  radiance  derivatives  which 
are  necessary  to  determine  the  radiance  statistics.  Inputs  to  the  model  include  the  statistical 
parameters  of  a  non-stationary  atmospheric  ten^erature  fluctuation  model  and  an  atmospheric 
profile.  The  radiance  statistics  are  used  in  a  simple  model  for  synthesizing  images.  The  model 
has  been  applied  to  calculate  the  radiance  structure  for  the  OH(Av=l)  swir  band  and  the  C02(V3) 
MWIR  band  under  nighttime  conditions. 
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SYNTHETIC  3-D  ATMOSPHERIC  TEMPERATURE  STRUCTURE: 

A  MODEL  FOR  KNOWN  GEOPHYSICAL  POWER  SPECTRA  USING 
A  HYBRID  AUTOREGRESSION  AND  FOURIER  TECHNIQUE 


James  H.  Brown  Neil  Grossbard 

Phillips  Laborarory  Boston  College 

29  Randolph  Road 
Hanscom  AFB,  MA  01731-3010 


Within  a  defined  domain,  geophysical  phenomena  often  are  characterized  by  smooth  continuous 
power  spectral  densities  having  a  negative  power  law  slope  dependence.  Frequently,  Fourier 
transform  analysis  has  been  employed  to  generate  synthetic  scenes  from  pseudorandom  arrays 
by  passing  the  stochastic  data  through  a  Fourier  filter  having  a  desired  correlation  structure  and 
power  spectral  dependency.  This  report  examines  the  possibility  of  producing  synthetic  structure 
by  invoking  autoregression  analysis  in  conjunction  with  the  Fourier  method.  Since  computations 
that  apply  multidimensional  fast  Fourier  transforms  to  large  data  arrays  consume  enormous 
resources  and  time,  the  goal  of  this  study  is  to  seek  an  alternative  method  to  reduce  the 
computational  burden.  Future  editions  of  the  Phillips  Laboratory  Atmospheric  Radiance  Codes 
will  feature  an  ability  to  calculate  structured  radiance.  The  metfiods  explored  herein  provide  a 
process  that  can  complement  or  in  some  cases  supplement  methods  presently  being  used.  Line 
of  sight  computations  using  the  three-dimensional  database  will  be  presented. 


PLANE  WAVE  SCINTILLATION  IN  AN  ONION  SKIN  MODEL 


Robert  R.  Beland 


Phillips  Laboratory 
29  Randolph  Road 
Hanscom  AFB,  MA  01731-3010 


The  geometry  of  an  onion  skin  model  for  atmospheric  optical  turbulence  is  presented  for 
application  to  near  horizontal  propagation.  This  is  then  used  to  derive  analytic  expressions  for 
the  log  amplitude  variance  of  a  plane  wave  in  the  weak  turbulence  (Rytov)  framework.  The 
geometrical  effects  on  the  altitude  weighting  function  are  investigated  and  used  to  define  the  limit 
of  validity  of  the  flat  earth  model.  The  formalism  is  then  extended  to  include  the  effects  of 
aperture  averaging  and  non-zero  timer  scale.  It  is  foimd  that  Rytov  theory  should  be  applicable 
to  zenith  angles  of  less  than  approximately  85°.  It  is  also  established  that  inner  scale  and 
aperture  averaging  are  much  less  important  for  near  horizontal  than  for  vertical  propagation  in 
the  weak  turbulence  case.  Predictions  of  log  amplitude  variance  are  presented  for  calculations 
using  the  CLEAR  t  model  as  well  as  individual  Cn2  profiles  derived  from  balloon  measurements. 
This  formalism  is  readily  extended  to  provide  analytic  expressions  for  other  turbulence  effects. 


COMPARISON  OF  A  MODEL  DESCRIBING  PROPAGATION  THROUGH  OPTICAL 
TURBULENCE  (PROTURB)  WITH  FIELD  DATA 


Robert  W.  Smith 

U.S.  Army  Test  and  Evaluation  Command 
Ft.  Belvoir  Meteorological  Team 
Ft.  Belvoir,  Virginia  22060 

Jennifer  C.  Ricklin 
U.S.  Army  Research  Laboratory 
Battlefield  Environment  Directorate 
White  Sands  Missile  Range,  New  Mexico  88002 

Kevin  E.  Cranston  and  John  P.  Cruncleton 
Physical  Science  Laboratory 
New  Mexico  State  University 
Las  Cruces,  New  Mexico  88003-0002 


The  PROTURB  model,  developed  by  the  Battlefield  Environment  Directorate  of  Army  Re¬ 
search  Laboratory,  calculates  an  estimate  of  optical  turbulence  strength  and  its  effects  on  visible 
and  infrared  imaging  and  laser  system  performance.  The  PROTURB  model  was  compared  to  data 
supplied  by  the  Army  Test  and  Evaluation  Command  Ft.  Belvoir  Meteorological  Team.  Here  we 
present  a  comparison  using  data  from  three  locations  that  differ  greatly  in  meteorological  condi¬ 
tions  and  in  site  characteristics.  We  discuss  the  data  collection  and  reduction  effort,  and  show  how 
this  data  compares  with  PROTURB  predictions. 


THE  ROLE  OF  TURBULENCE  IN  CLOUD  DROPLET  FORMATION 

AND  OUTSIDE  CLOUDS 


4 


by  James  W.  Telford 
Atmospheric  Sciences  Center 
Desert  Research  Institute 
P.O.  Box  60220 
Reno  NV  89506 


Turbulence,  it  is  becoming  increasingly  clear,  contains  organized  structures  which  carry  fluid 
from  one  place  to  another  through  the  turbulent  medium.  Much  of  the  atmosphere  is  in  a  stable 
stratified  state.  In  these  conditions  two  dimensional  vortices  form,  and  smaller  swirls  can  coalesce 
to  feed  larger  scale  motion.  The  Red  Spot  on  Jupiter  appears  to  be  such  a  structure.  In  clouds, 
cooled  parcels  near  cloud  tops  generate  descending  spherical  vortices  which  carry  diluted  cloud 
parcels  deep  into  the  cloud.  This  phenomena  drives  the  microphysical  development  in  both  cumu¬ 
lus  and  stratus  clouds,  and  hence  controls  the  droplet  spectra  and  the  distribution  of  ice,  which  in 
turn  plays  a  large  role  in  the  radiative  characteristics  of  the  cloud.  It  seems  likely  that  the  mammatus 
protruding  into  the  clear  air  below  shelf  clouds  are  of  this  form  of  motion. 
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SHARC-3:  A  MODEL  FOR  INFRARED  ATMOSPHERIC  RADIANCE  AT  HIGH 
ALTITUDES 


R.D.  SharmaS  J.H.  Brown^  J.H.  Gruninger^  R.L.  Sundberg^,  J.W.  DuJf^,  L.S.  Bernstein^ 
M.W.  Matthew^,  S.M.  Adler-Golden^  D.C.  Robertson^  and  R  J.  Healey^ 


Thillips  Laboratory 
29  Randolph  Road 
Hanscom  AFB,  MA 


Spectral  Sciences  Inc. 
99  South  Bedford  St.,#7 
Burlington,  MA  01803 


Yap  Analytics,  Inc. 
594  Marrett  Road 
Lexington,  MA  02173 


SHARC-3  is  the  latest  version  of  Strategic  High  Altitude  Radiance  Code.  This  talk 
describes  new  SHARC  features  like  its  ability  to  simulate  changing  atmospheric  conditions  along 
the  line-of-sight  (LOS)  paths  being  calculated.  SHARC  has  been  developed  by  the  Air  Force  to 
meet  systems  requirements  for  a  simulation  code  that  is  both  suitable  for  systems  studies  and 
which  provides:  rapid  and  accurate  calculations  of  upper  atmospheric  IR  radiance  and 
transmittance  spectra;  has  a  spectral  resolution  of  better  than  1  cm'*;  spans  the  2  to  40  pm  (250 
to  5,000  cm'*)  wavelength  region;  and  accommodates  arbitrary  lines  of  sight  in  the  50-3(X)  km 
altitude  regime.  SHARC  accounts  for  the  production,  loss,  and  energy  transfer  processes  among 
the  molecular  vibrational  states  important  to  this  spectral  region.  Auroral  production  and 
excitation  of  COj,  NO,  and  NO^  are  included  in  addition  to  quiescent  atmospheric  processes. 
Calculated  vibrational  temperatures  are  found  to  be  similar  to  results  from  other  non-LTE  codes, 
and  SHARC’s  equivalent-width  spectral  algorithm  provides  very  good  agreement  with  much  more 
time-consuming  "exact"  line-by-line  methods.  Calculations  and  data  comparisons  illustrating  the 
features  of  SHARC-3  are  presented. 


NON-LTE  IN  CO2  IN  THE  MIDDLE  ATMOSPHERE 


A.A.  Kutepov 


V.P.  Ogibalov 
G.M.  Shved 


Institute  for  Astronomy  &  Astrophysics 

University  of  Munich 

Germany 


Dept,  of  Atmospheric  Physics 
University  of  Sl  Petersburg 
Russia 


A  new  line-by-line  ALI  (accelerated  lambda  iteration)  technique,  initially  designed  to 
solve  multilevel  non-LTE  line  formation  problem  in  stellar  atmospheres,  is  applied  to  smdy 
effects  of  the  non-LTE  on  the  populations,  source  functions  and  limb  radiances  in  various 
infrared  CO2  bands  in  Earth’s  middle  atmosphere.  Approximations  used  in  previous  studies  of 
non-LTE  in  CO2  are  analyzed.  Implications  for  remote  sensing  with  high  spectral  resolution 
instruments  are  discussed. 


COMPARISON  OF  LINE-BY-LINE  AND  MODIFIED  CURTIS  MATRIX  NARROW- 
BAND  MODEL  APPROACHES  TO  RADIATIVE  TRANSFER  IN  THE  CO2 15  BANDS 


P.P.  Wintersteiner*,  M.  Lopez-Puertas^  J.R.  Winick^,  and  R.H.  Picard^ 


‘Arcon  Corp. 
Waltham,  MA  02154 


^Institute  de  Astrofisica 
de  Andalucia 
Granada,  Spain 


^Phillips  Laboratory 
29  Randolph  Road 
Hansoom  AFB,  MA  01731-3010 


We  have  undertaken  a  detailed  comparison  of  two  different  models  for  calculating  non- 
LTE  populations  in  the  terrestrial  atmosphere,  for  the  purpose  of  studying  the  effects  of  the 
radiative  transfer  algorithms  that  are  embodied  in  them.  The  models  are  the  iterative  line-by-line 
RAD  model  that  forms  the  basis  of  the  PL  ARC  code  and  the  inverse  modified  Curtis  matrix 
(MCM)  model,  developed  by  Lopez-Puertas  and  co-workers.  They  have  been  used  to  calculate 
the  vibrational  temperatures  and  cooling  rates  in  the  COj  15  pm  bands.  The  comparison  covers 
the  altitude  range  40-160  km,  includes  strong  and  weak  bands,  uses  model  atmospheres  having 
extreme  temperature  structures,  and  also  tests  the  effect  of  changing  kinetic  parameters.  The 
models’  vibrational  temperatures  and  cooling  rates  generally  agree  to  within  2  K  and  0.5  K/day, 
respectively.  We  will  show  these  results,  discuss  those  aspects  of  the  radiative  transfer 
algorithms  responsible  for  the  small  differences  that  do  appear,  and  mention  the  implications  for 
non-LTE  code  validation. 


SUBTHERMAL  NOell^)  SPIN-ORBIT  DISTRIBUTIONS  IN  THE  THERMOSPHERE 


S.  J.  Lipson,  P.S.  Armstrong',  J.A.  Dodd',  J.R.  Lowell',  W.A.M.  Blumberg, 

and  R.M.  Nadile 


Optical  Environment  Division,  Geophysics  Directorate,  Phillips  Laboratory. 

Hanscom  AFB.  MA  01731 

^Stewart  Radiance  Laboratory,  139  Great  Road,  Bedford,  .MA  01730 


The  populations  in  the  two  NO(X^n,v=l,Q=l/2, 3/2)  spin-orbit  manifolds  in  the  Earth’s  ther¬ 
mosphere  have  been  found  to  depart  by  as  much  as  a  factor  of  two  from  the  ratio  expected  from 
thermal  equilibrium.  The  effective  temperature  describing  the  observed  population  distribution  is 
hundreds  of  degrees  Kelvin  lower  than  the  kinetic  or  rotational  temperatures.  Absolute  spin-orbit 
specific  colunm  densities  of  NO(v=l)  were  derived  from  high-resolution  (1  cm'')  IR  earthlimb 
spectra  in  the  100  to  200-km  region,  obtained  in  the  CIRRIS  1 A  Space  Shuttle  experiment.  Nonlin¬ 
ear  least-squares  synthetic  spectral  fitting  was  used  to  analyze  the  NO  Av=l  fundamental  band 
emissions  near  5.3|im.  The  spin-orbit  distribution  represents  a  third  degree  of  freedom,  along  with 
vibration  and  rotation,  that  is  not  in  equilibrium  with  the  local  kinetic  temperature.  These  observa¬ 
tions  may  significantly  impact  the  remote  sensing  of  NO  in  the  upper  atmosphere.  The  subthermal 
distribution  is  most  likely  produced  in  the  collisional  uppumping  of  NO(v=0)  by  O  atoms,  which  is 
the  major  source  of  NO(v=l)  in  the  thermosphere.  Thus,  spin-orbit  distribution  measurements  may 
provide  a  means  of  distinguishing  NO  produced  by  chemical  reaction  from  NO  excited  by  colh- 
sional  uppumping. 

This  work  was  supported  by  the  Air  Force  Office  of  Scientific  Research  and  the  Ballistic 
Missile  Defense  Organization. 
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COMPARISON  OF  GLOBAL  VARIATION  OF  THE  RADIANCE  OF  THE 
STRATOSPHERIC  O3  CH4  SPECTRA  AS  VIEWED  CIRRIS  lA  FORM  DISCO VERY-103 
WITH  CLAES  EXPERIMENT  ON  BOARD  OF  UARS 


B.K.  RezaP,  G.E.  Bingham^  L.R.  MegilP,  D.K.  Zhou^  M.  Mergenthalei^,  A.E.  Roche^, 
J.B.  KumaH,  and  G.P.  Anderson^ 


^Center  for  Atmospheric  &  Space  Sciences 
Utah  State  University 
Logan,  UT  84322-4405 


^Space  Dynamic  Laboratory 
Utah  State  University 
Logan,  UT  84322-9710 


^Lockheed  Palo  Alto  Research  Lab 
Palo  Alto,  CA  94304 


Phillips  Laboratory/GP/GPOS 
29  Randolph  Road 
Hanscom  AFB,  MA  01731-3010 


We  report  on  meaurements  of  these  species  by  the  Cryogenic  Infrared  Radiation 
Instrumentation  for  shuttle  (CIRRIS  lA)  which  flew  aboard  the  Discovery  OV-103,  STS-39  on 
April  28  to  May  6,  1991.  By  looking  in  the  limb,  it  was  possible  to  obtain  approximately  a 
thousand  interferograms.  This  made  possible  a  certain  identification  of  many  important 
stratospheric  gases,  including  O3,  CH4,  HNO3,  and  chlorine  source  species  CCLjFj,  and  CCL3F. 
The  acciuracy  of  the  CIRRIS  lA  measurements  is  therefore  of  interest  in  polar  ozone  depletion 
studies,  and  validation  is  a  major  consideration.  This  paper  will  review  the  CIRRIS  lA 
observations  of  O3,  CH4,  and  HNO3  spectra  and  of  vertical  profiles  concentration  and  discuss 
their  accuracy  and  precision.  We  will  compare  the  CIRRIS  1 A  measurements  with  the  Cryogenic 
Limb  Array  Etalon  Spectrometer  (CLAES)  on  board  the  Upper  Atmosphere  Research  Satellite. 
Specifically  we  will  compare  CIRRIS  1 A  and  CLAES  spectra  in  the  wavelength  regions  targeted 
for  O3,  CH4,  and  HNO3  concentration  retrievals,  and  we  will  compare  the  retrieved  quantities  for 
species  in  the  southern  and  northern  hemispheres.  Since  CIRRIS  1 A  data  were  collected  shortly 
before,  and  CLAES  data  shortly  after,  the  Pinatubo  eruption,  such  comparisons  are  of  significant 
value  both  scientifically  and  from  a  validation  standpoint. 
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THE  POAM  EXPERIMENT  AND  EARLY  MEASUREMENT  RESULTS 


E.P.Shettle,  R.M.Bevilacqua,  J.S.Hornstein  (Code  7227,  NRL,  Washington,  DC  20375) 

WJ.Glaccum  (ARC,  Landover,  MD) 

S.Krigman,  JXumpe,  M.Fronim,  and  D.Debrestian  (CPI,  Fairfax,  VA) 
and  the  POAM  II  Science  Team 


NRL’s  Polar  Ozone  and  Aerosol  Measurement  (POAM  II)  instrument  was  developed  to  mea¬ 
sure  the  distribution  of  atmospheric  aerosols,  PSC’s,  and  several  molecular  species  critical  for 
understanding  ozone  chemistry  in  the  polar  stratosphere.  It  was  launched  on  France’s  SPOT  3 
satellite  on  25  September  1993.  It  makes  solar  occultation  measurements  in  nine  narrow-band 
channels  from  the  visible  to  the  near  IR.  From  the  measured  transmissions,  it  is  possible  to  deter¬ 
mine  the  densities  of  aerosols  at  altitudes  of  10-40  km,  ozone  at  altitudes  of  10-60  km,  water  vapor 
at  altitudes  of  15  -40  km,  and  nitrogen  dioxide  at  altitudes  of  20-40  km.  The  SPOT  3  satellite  is  in 
a  98°  inclined  sun-synchronous  orbit  at  an  altitude  of  833  km.  This  results  in  all  the  POAM  n 
measurements  being  in  the  Polar  regions.  We  will  discuss  the  POAM  n  instrument  design,  the 
retrieval  algorithm,  and  the  preliminary  results  of  the  measurements,  including  early  data  valida¬ 
tion  studies.  POAM  II  is  sponsored  by  the  Innovative  Science  and  Technology  Office  of  BMDO 
(BMDO/DTI)  and  the  DoD  Space  Test  Program  (STP). 
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LARGE-SCALE  RETRIEVAL  OF  ATMOSPHERIC  PARAMETER  PROFILES 


L.  Sparks 
J.  L.  Fanselow 
J.  McComb 
S.  Nandi 
J.  Parker 
J.  E.  Patterson 
H.  Worden 

Jet  Propulsion  Laboratory 
Pasadena.  California 


We  have  developed  and  tested  SEASCRAPE  (Sequential  Evaluation  Algorithm  for  Simulta¬ 
neous  and  Concurrent  Retrieval  of  Atmospheric  Parameter  Estimates),  a  code  designed  to  automate 
and  expedite  the  retrieval  of  atmospheric  profiles  from  large  quantities  of  remote  sensing  data 
generated  by  nadir  and  limb  soundings.  Our  algorithm  uses  a  sequential  parameter  estimation  method 
known  as  the  Square  Root  Information  Filter,  a  numerical  technique  for  solving  least  squares  prob¬ 
lems  that  is  particularly  suited  for  handling  large  volumes  of  data.  The  algorithm  permits  simulta¬ 
neous  retrieval  of  many  atmospheric  profiles.  By  implementing  the  algorithm  in  a  paraUel  process¬ 
ing  environment,  concurrent  data  analysis  will  achieve  a  dramatic  increase  in  the  rate  at  which  the 
data  can  be  processed. 


GENERAL  INVERSION  PACKAGE  FOR  RETRIEVAL  OF  ATMOSPHERIC  SPECIES 
FROM  HIGH  RESOLUTION  SPECTRAL  MEASUREMENTS 


J.L.  Moncet  and  W.O.  Gallery 


Aonospheric  &  Environmental  Research,  Inc. 
840  Memorial  Drive 
Cambridge,  MA  02139 


A  general  inversion  algorithm  based  on  the  Rodgers  (1976,  1990)  optimal  inversion 
approach  is  being  developed.  The  algorithm  is  applicable  to  the  simultaneous  retrieval  of  any 
ensemble  of  selected  molecular  species  as  well  as  atmospheric  temperature  from  down-looking, 
up-looking  or  limb  measurements.  The  algorithm  has  been  designed  to  use  FASCODE  for  the 
molecular  optical  depths  calculations  but  can  be  modified  to  run  with  faster  transmittance  models 
(e.g.  MODTRAN).  Substantial  computational  gain  is  achieved  over  the  existing  "Path 
Characterization"  algorithm  (Isaacs  et  al.  1990),  by  performing  the  K-matrix  calculation 
analytically  and  by  making  the  code  fully  vectorizable.  The  algorithm  has  been  applied 
successfully  to  the  simultaneous  retrieval  of  temperature  and  water  vapor  from  the  ground  based 
interferometric  measurements  (Theriault  et  al.  1993).  The  algorithm  will  be  applied  to  CIRRIS 
data  and  to  MSX  data  when  available. 
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SIMULATION  OF  STELLAR  OCCULTATION  MEASUREMENTS 


L.  Oikarinen 


Finnish  Meteorological  Institute 
Dept  of  Geophysics 
P.O.  Box  503,  FIN-00101 
Helsinki,  Finland 


An  instrument  for  monitoring  stratospheric  ozone,  GOMOS  (Global  Ozone  Monitoring 
by  Occultation  of  Stars),  will  be  launched  by  ESA  on  the  ENVISAT-1  satellite  in  1998.  The 
instrument  measures  atmospheric  UV-visible  transmission  spectra  in  250  -  675  nm  at  a  resolution 
of  0.6  nm.  Absorption  spectra  of  Oj  and  HjO  are  detected  by  two  IR  channels  of  spectral 
resolution  0.12  nm  located  around  760  and  930  nm.  Transmission  along  slant  paths  of  tangent 
heights  down  to  15  km  is  measured.  Solar  light  scattered  into  the  instrument  FOV  forms  a 
background  signal.  Simulation  of  both  the  transmittance  and  the  scattered  background  are 
required  in  developing  the  instrument  and  data  analysis  methods.  Special  requirements  of 
occultation  measurements  are  the  modeling  of  scintillations  and  refractive  attenuation  of  light 
(dilution  of  the  flmc  density  due  to  the  decreasing  of  refractive  index  with  altimde).  In  the  case 
of  GOMOS  absorption  by  minor  species,  such  as  NOj,  NO3,  BrO,  and  OCIO,  should  also  be 
simulated. 
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LINE-BY-LINE  CALCULATIONS  OF  ATMOSPHERIC  FLUXES  AND  COOLING 
RATES  II:  APPLICATION  TO  CARBON  DIOXIDE,  OZONE,  METHANE,  NITROUS 

OXIDE,  AND  THE  HALOCARBONS  ^ 


M.  J.  lacono,  and  S.  A.  Clough 
Atmospheric  and  Environmental  Research,  Inc. 
840  Memorial  Drive 
Cambridge,  MA  02139 

Tel:  617-547-6207 
FAX:  617-661-6479 


Determining  the  spectral  contribution  of  atmospheric  molecular  components  to  the  total 
longwave  fluxes  and  cooling  rates  is  a  necessary  step  toward  understanding  the  effect  and  re¬ 
sponse  of  these  gases  to  global  climate  change.  Clear-sky  longwave  computations  using  a  line- 
by-line  radiative  transfer  model  (LBTRTM)  have  been  completed  for  several  atmospheric  pro¬ 
files.  We  have  previously  presented  results  for  both  water  vapor  [Clough  et  al. ,  19921]  and  the 
water  vapor  continuum  [Clough  et  .al.,  1989].  In  the  current  paper  we  present  the  effects  of 
carbon  dioxide,  ozone,  methane,  nitrous  oxide  and  the  most  important  halocarbons  in  addition  to 
water  vapor.  Both  the  spectral  and  spectrally  integrated  flux  and  cooling  rate  contributions  from 
each  of  these  significant  greenhouse  gases  are  evaluated.  Carbon  dioxide,  an  extremely  strong 
atmospheric  absorber,  is  generally  saturated  in  the  troposphere  with  the  exception  of  two  narrow 
spectral  bands  with  the  result  that  carbon  dioxide  plays  an  anomolously  small  role  in  this  altitude 
regime.  Additional  studies  include  the  effect  on  the  longwave  cooling  rates  and  fluxes  due  to  the 
anticipated  ten  year  changes  in  C02, 03,  CFC-11,  CFC-12,  GFC-22,  and  CCL4  from  1990  to  2000 
as  well  as  due  to  the  doubling  of  carbon  dioxide 


Clough,  S.A .,  F.X.  Kneizys,  and  R.  W.  Davies  Line  shape  and  the  water 
vapor  continuum,  Atmos.  Res.,  23,  229-241,  1989. 

Clough,  S.  A,  M.  J.  lacono,  and  J-L.  Moncet.  Line-by-line  Calculations 

or  Atmospheric  Fluxes  and  Cooling  Rates:  Application  to  water  Vapor.  ^ 

J  Geophys.  Res  .,  97,  15761-15785, 1992. 
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VERY  NARROW  BAND  MODEL  CALCULATIONS  OF 
ATMOSPHERIC  FLUXES  AND  COOLING  RATES 
USING  THE  MODTRAN  CODE 


L.  S.  Bernstein,  A.  Berk,  P.  K.  Acharya,  and  D.  C.  Robertson 
Spectral  Sciences,  Inc.,  Burlington,  Massachusetts 

and 

G.  P.  Anderson,  J.  H.  Chetwynd,  and  L.  M.  Kimball* 

Geophysics  Directorate,  Phillips  Laboratory 
Hanscom  AFB,  Bedford,  Massachusetts 


Because  of  the  large  computational  speed  advantage  afforded  by  band  models  over  line-by¬ 
line  (LBL)  models  the  development  and  application  of  band  models  to  the  calculation  of  atmo¬ 
spheric  radiative  fluxes  and  cooling  rates  continues  to  be  an  active  area  of  research.  A  wide  variety 
of  narrow  band  model  (NBM)  and  wide  band  model(WBM)  approaches  have  been  developed. 
Typically  the  calculation  of  spectral  resolution.  Aw,  of  the  NBMs  falls  within  the  5<Aw<20  cm  ' 
range,  while  for  the  WBMs  Aw>100  cm  '.  Direct  application  of  these  band  models  to  the  calcula¬ 
tion  of  atmospheric  fluxes  and  cooling  rates  usually  results  in  poor  to  only  fair  agreement  with  LBL 
calculations.  In  order  to  circumvent  this  problem  the  band  models  are  treated  as  parameterizations 
rather  than  as  first  principles  models.  These  tuned  band  models  can  yield  atmospheric  flux  and 
cooling  rate  predictions  of  LBL-like  quality.  This  paper  discusses  a  new  first  principles  based  very 
narrow  band  model  (VNBM)  approach  with  Aw=  1  cm''  which  yields  atmospheric  flux  and  cooling 
rate  predictions  of  near  LBL  quality.  The  VNBM  is  incorporated  in  the  Phillips  Laboratory 
MODTRAN  atmospheric  transmittance-radiance  code  which  enables  general  specification  of  the 
atmospheric  species,  concentration  and  thermodynanuc  profiles,  calculational  spectral  interval, 
and  path  geometry.  The  relatively  high  spectral  resolution  of  this  approach  combined  with  separa¬ 
tion  of  the  transmittance  algorithm  into  truncated  line  and  line  wing  components  treats  the  problem 
of  line  correlation  in  a  physically  realistic  manner.  In  order  to  evaluate  the  VNBM  approach  com¬ 
parison  is  made  to  LBL  benchmark  atmospheric  flux  and  cooling  rate  calculations. 


'  Mathematics  Department,  Worcester  Polytechnic  Institute 
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GREENHOUSE  GAS  CONCENTRATION  PROFILES  RETRIEVED  FROM  CIRRUS  lA 
MEASUREMENTS  IN  THE  7-13  ,um  WINDOW  DURING  STS-39 
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G.E.  Bingham,  D.K.  Zhou,  Y.  Yang,  and  AJ.  Steed  (All  at  Space  Dynamics  Lab., 
Utah  State  University,  Logan,  UT  84322:  801  770  0256:  e-mail:  cimis@cc.usu.edu) 
G.P.  Anderson  and  R.M.  Nadile  (Geophysics  Directorate,  Phillips  Laboratory, 
Hanscom  Air  Force  Base,  MA  01731) 


The  spectral  and  spatial  distributions  of  the  greenhouse  gas  infrared  emissions  are  mea¬ 
sured  by  the  CIRRUS  1 A  (Cryogenic  Infrared  Radiance  Instrumentation  for  Shuttle)  inferometer 
sensor  during  STS-39  flight  (28  Apr  -  6  May,  1991).  The  spectra  in  the  important  7- 1 1  |im  infrared 
window  were  recorded.  High  resolution  spectral  scans  (-1.5  cm'^  instrument  resolution)  and 
FASCOD-3  with  HITRAN92  database  and  NMC  temperature  and  pressure  profiles  are  used  in  an 
onion  peel  routine  to  retrieve  gas  concentrations  from  absolutely  calibrated  radiance  data.  Scan  to 
scan  and  profile  to  profile  variability  are  compared  to  show  the  measurement  capability  of  the 
cryogenic  infrared  interferometer.  The  local  aerosol  is  estimated  by  comparing  the  modelled  and 
measured  CO,  hot  band  (10.5  ,Lim)  spectra.  The  results  of  concentration  profiles  retrieved  from  the 
spectra  (42°  south  latitude  32°  west  longitude,  local  time  of  22.4  hr.)  are  shown.  Profiles  of  CFC- 
11,  CFC-12,  HNO-,  O3  and  CH^  developed  in  this  fashion  and  at  different  locations  are  compared 
with  ATMOS  (Atmospheric  Trace  Molecule  Spectroscopy),  LIMS  (Limb  Infrared  Monitor  of  the 
Stratosphere),  and  situ  mean  data  for  these  gases. 
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Changing  concentrations  of  atmospheric  methane  can  significantly  affect  the  concentration 
and  distribution  of  ozone  in  the  troposphere  and  stratosphere.  Changes  in  the  concentrations  of 
these  gases  can  thus  influence  the  radiative  forcing  on  climate.  In  this  study  we  determine  the 
indirect,  time-varying  tropospheric  radiative  forcing  and  Global  Warning  Potential  (GWP)  of  O3 
due  to  an  impulse  of  CH^  that  is  equivalent  to  36  percent  of  the  atmospheric  mass  of  CH^.  The  direct 
CH  radiative  forcing  and  GWP  are  also  calculated.  The  LLNL  2-D  radiation-chemistry-transport 
code  is  used  to  evaluate  the  resulting  changes  in  the  O3  and  CH^  atmospheric  profiles  as  a  function 
of  time.  A  correlated  k-distribution  radiative  transfer  model  is  used  to  calculate  the  radiative  forc¬ 
ing  at  the  tropopause  of  the  globally-averaged  atmosphere  profiles.  The  O3  indirect  GWPs  vary 
from  ~  25  after  20  yr  integration  to  ~3  after  500  years.  The  direct  CH^  GWP  vary  from  "48  at  20  yrs 
to  ~  6  at  500  yrs.  The  20  year  direct  CH^  GWP  is  -20%  larger  than  previous  estimates  due  to  an 
increased  CH^  lifetime  resulting  from  the  applied  impulse.  Estimates  are  made  of  latitude  effects  in 
the  radiative  forcing  calculations. 
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